﻿Editura Didactică și Pedagogică MINISTERUL EDUCAȚIE! Șl iNVATAMiNTULU! Prof dr ing VLAD IONESCU sisteme liniare Ш Q ț у Editura Didactică șl Pedagogică București torni N ч К PREFAȚĂ Luciarea prezintă, cu specificul și cerințele impuâe unui curs universitar, al teoriei sistemelor liniare, orientat cu precădere spre bazele Elaborarea textului s-a bazat în primul rînd pe încercarea autorului ele conținutul uzual al sistemelor automate și proiectarea sistemică a acestora" a satisface trei exigențe pe care le reclamă în prezent eficiența procesului instructiv: — asigurarea unui cadru definitoriu și metodologic capabil să faciliteze o pregătire în profil larg și adaptabilă la dinamica informațională actuală — prezentarea numai a acelor noțiuni și metode care eminamente stimulează creativitatea studentului, cu evitarea abuzului de informație inutilă' — orientarea către aplicații ce contribuie efectiv la dezvoltarea deprinderilor și a unei gîndiri care să permită viitorului specialist să rezolve complexele probleme impuse de practică în contextul permanentului progres științific și tehnologic Plecînd de la dezideratele menționate, conținutul cursului a urmat o dezvoltare care să cuprindă, cu amendamentele unei asimilări comode și cu o pondere adecvată, atît metodele inginerești devenite tradiționale în studiul sistemelor automate (cum sînt cele frecvenționale), cît și cele moderne, indispensabile proiectării asistate prin calculator (cum sînt abordările de tip geometric și polinomial și care acoperă în cvasitotalitate întreaga problematică a sistemelor liniare) Capitolul , evitînd o axiomatică aridă, introduce conceptul de sistam privit din punctul de vedere țelor sale Capitolul începe prin a actualiza sub raportul abordării sistemice implicit al rigurozității, o serie de noțiuni clasice din tehnica tradiționali , a reglării automate Totodată, au fost prezentate într-o manieră Wicata atit L pe parcursul Jrării, esențiala atenție ui"lnct“marcabS de fuziune intre sistemele continue și discrete l-a re-p JLtat analiza unificau îna®foVpmzeS unele și efectele stabilizante alo acesteia ț 'm itica CtfUM are o importanții ing nerease; țMggg ■ л ' distinctlee »**!'■**-»* și atributele acestuia plecmd de la noțiunea de proces, acesta din urmă fiind privit din punctul de vedere al funcționalității tehnologice și al perfonm cursului este cu precădere de reiecție (asimptotică) a perturbațiilor persistente Totodată cititorul este familiarizat cu noțiunile de reglare cu localizarea perturbațiilor și regxare robustă Sub raportul proiectării, conforme cu familia performanțelor convenționale, s-au evidențiat o serie de aspecte de universalitate ale metodelor in Capitolul reprezintă trecerea efectivă de la cadrul abordărilor tradi-ționale la punctul de vedere sistemic modern cu evidențierea expresă a proprietăților structurale remarcabile și a „geometrizării" acestor proprietăți Partea finală a capitolului introduce cititorul in metodele de abordare, geometrică metode care și-au dovedit eficiența procedurală și în cazul neliniar Capitolul , și ultimul, formulează și rezolvă cîteva din problemele semnificative ale sintezei sistemice cum sînt: problema acoperirilor dinamice alocate (care stă la baza sintezei estimatoarelor), problema stabilizării (prin compensarea cu extensie dinamică), problema reglării Referitor la această ultimă problemă s-a insistat mai mult pe problema reglării intern stabile și a metodelor structurale de rezolvare a acesteia Demonstrațiile matematice apelează strict la programele analitice ale cursurilor de matematici (analiză și matematici speciale), numeroase exemple ușurînd cititorului înțelegerea și subtilitățile diferitelor concepte și metode De asemenea, în alcătuirea materialului, în dezvoltările și structurarea acestuia, s-a avut permanent în vedere completarea și corelarea dintre textul scris și cursul verbal expus, dialogul obligatoriu care trebuie să existe în procesul de asimilare între cadrul didactic și student, dialog în mod firesc extins și în aria valențelor sale educative Așa cum s-a menționat, efortul metodic axat pe automatică, mai exact pe descrierea sistemică a conceptelor și în special a metodelor acesteia S-a insistat pe prezentarea, inclusiv algoritmică, a acelor metode care permit studentului, viitorului inginer automatist, să înțeleagă și să-și însușească ușor particularitățile proiectării asistate cu care va lua contact la disciplinele din anii terminali (Proiectarea structurilor de reglare Sisteme și echipamente pentru conducerea proceselor Roboți indus triali etc ) Această preponderentă axare a cursului pe automatică, evident încadrată în contextul general de pregătire al tuturor disciplinelor care definesc specialitatea Automatizări și calculatoare, rezultă ca un imperativ major firesc al Directivelor Congresului al XlII-lea al Partidului Comunist Român cu privire la dezvoltarea economico-socială a României, a sublinierii exprese a rolului ce revine automatizării, electronizării și robotizării în creșterea nivelului tehnic și calitativ al producției în acest sens, eficiența aportului învățămîntulm tehnic, în particular al celui de automatică și calculatoare, la traducerea în viața a obiectivelor menționate constă, și aceasta avînd o valoare axiomatică, în procesul de integrare al învățămîntului cu cercetarea -Și producția, cu practica socială, concept original, revoluționar și științific ■elaborat de tovarășul Nicolae Ceaușescu, secretarul general al partidului Autorul ține să mulțumească tovarășului șef lucrări dr ing Șerban Sever pentru permanentele discuții purtate pe marginea materialului și pentru ec ura atentă a manuscrisului cursului Aceleași mulțumiri sînt aduse, pentru concursul dat tovarășilor șef lucrări Florin Stratulat și Soare Călin, precum și tovarășului asistent I Corâci F cărei bfinS î Selrite * > * CUPRSNS Prefață CAP Sisteme dinamice Definiții Clasificări Proces și sistem Definirea noțiunii de sistem dinamic Clasificări Structuri particulare Axiomatica CAP Sisteme liniare invariante în timp, съі o intrare și o ieșire Considerații generale Comportarea intrare/ieșire Funcția de transfer Realizări Echivalență Stabilitatea sistemelor liniare cu o intrare și o ieșire Răspunsul sistemelor liniare la mărimi de intrare standard Discretizarea sistemelor netede Reprezentarea (în frecvență) a funcțiilor de transfer Corelația dintre modulul și argumentul funcției de transfer Determinarea răspunsului indicial din reprezentarea în pulsație a funcției de transfer Conexiunea sistemelor Reacția negativă Stabilitatea conexiunii în reacție * negativă Sisteme diferențiale și cu timp mort ' CAP Sisteme de reglare automată Formularea problemei reglării Considerații și rezultate cu caracter general Analiza sistemelor de reglare automată Analiza stabilității (asimptotice interne) Performanțele regimului dinamic Precizia sistemelor de reglare automată Performanțele regimului staționar (al erorii) Sinteza convențională a sistemelor de reglare automată Metode de sinteză exactă Metode de sinteză în frecvență G, (dezacumulare daca ct cu #сШНяХІГхН; SKczlSixlSA я> y, u, v păstrîndu-și semnificația de mai sus (v fig ) Pentru o reprezentare mai compactă mărimea de comandă și cea perturbatoare se reunesc și formează mărimea de intrare, iar cu mărimea măsurată și cea de calitate se procedează similar, rezultînd mărimea de ieșire In consecință ( ) devine £ = X, Ы) У = g(t, X) ( - ) unde cu «еЕГ , veR” s-au notat (în mod abuziv) mărimea de intrare și respectiv de ieșire Introducem acum: Д Definiția Numim sistem dinamic neted tripletul ăvînd semnificația matematică ( ) cu ^( •) e A'o, -« ( ’)) ( ) Mulțimea { \ avtad semnificația matematică J x(t -I- ) = f{ț, x(l), u(t)) y(l) = g{t, x(l)) ( ) cu t e Ж, unde i’) t%A{M( ) JE L”' —> JI LW continuă în raport cu x și и Adesea vom utiliza, în special pentru sistemele liniare, scrierea xr = Д/, x, u) (L ) У = s(t x) unde, în mod convențional, A Definiția Un sistem dinamic (neted sau discret) se numește liniar dacă \ A(f) e cu cu f este liniară în x și ut adică /(/, x, u} = A(t)x - B(t)u MXW, B(/) g JELwxw (continue dacă £eX L) g este liniară în x adică g(t, x) = C(t)x ) e JESLpxn (continuă dacă /eJR ) Cititorul este invitat să constate că dacă procedăm la desfacerea mărimii • » ; ’ : • ' ’ • ■ • AXIOMATICA Construcția axiomatică copiază proprietățile soluției ( ) ale ecuației diferențiale din ( ), Aceasta se face cu următoarele elemente constitutive: , mulțimea momentelor de timp, cu ordonarea naturală din Uzual T — ”*} — clasa funcțiilor da intrare n — „spațiul stărilor" ) |y(-): T-> У c } “* clasa ieșirilor / • ’ ’ - * ; : ■ e ’ • ’ • ч • " • ’ * * • • , • •* t X X x QZ —♦ A — funcția de tranziție a stărilor Y — funcția de ieșire, unde U și Y introduse prin și sînt eub mulțimi ' v Se consideră adevărate următoarele axiome: Al Netrivialitatea czie exprimă faptul că ( ) A Concatenaritatea Fie e T, t == , , cu tt (/; h> фй; zo> xo> *(•))> “(')) pentru V хее X și u(-) e A Cauzalitatea Fie «(•), u'(?) -У două intrări arbitrare pentru care kn[ u (’)) Vv e F Л [Zq> ( ) ( ) ( ) Avem astfel Definiția Se numește sistem dina ic octuplul S = (T, U, Y, у, X, Л%( V^Cj , De asemenea, prin substituții succesive se obține din prima ecuație a lui ( ): x(ț) = (Z; % , «(•)) — Alxe - P Al~T >u(t), Z e IV ( ) т^О Punînd ( - ) unde Ф(/) se numește matricea fundamentală, egalitățile ( ) și ( ) se re-scriu în forma rt %(/) == Ф(^) x l ф(£ — т) w(x)dT, t G JHL ( ) Jo și respectiv £ —* x{i) = ®(t)x - Ф(і — t — ) bu(x), t e M Тм ( ) Sa observăm că în cazul sistemelor netede matricea fundamentală este nesingulară oricare ar fi matricea A, iar Ф" ^) = Ф(—t), We Et* în cazul sistemelor discrete însă, Ф(/) este nesingulară dacă și numai dacă Л este nesingulară, în care caz, de asemenea, Ф’х(і') = Ф(—t), Wg Față de aceste constatări, obținem din ( ) Г* Xq zsx %( ) я= Ф( — t) x(t) — l Ф(—т) u(t) df «■o iar din ( , ) - x( ) ? ф(—t') %■(/) — ^ ) Ф( — v— ) u(v) ( ) ( ) - = Cum A este de tip companion, matricea modală este de tip Vandermonde și deci forma Jordan a matricei A este Atunci unde Pentru cazul discret avem ФШ A Și deci ' (— )‘ Ф( ) « > ф( ) > Ф(/) «=a f - •> Întrucît A este singulari, Ф(/) cu t ( ) și respectiv y(t) сФ(/)% -p ( ) n I w* hc (Щ ) И) ret ) ) ) nsi- ) •ll! И) ) ) Atunci, utihzînd produsul de convoluție*, relațiile ( ) și ( ) serescriu unificat în forma ’ MO = стФ(і) x + (Лс* u) (/), /> ( ) cu ZelRL sau te'Z,, după cum sistemul este neted sau discret, obținîndu-se o exprimare compactă a răspunsului sistemelor liniare Primul termen din membrul drept al lui ( ) definește componenta liberă a răspunsului, У’М = x , t^Q ( ) •• * * • asociată condiției inițiale x ; cel de-al doilea termen din membrul drept al lui ( ) definește componenta forțată a răspunsului (V й) w i > ( ) asociată comenzii (intrării) u(-) Se constată imediat din ( ), ( ) că jy(O) = Așadar, cu ( ) și ( ), ( ) se poate scrie ca MO =>»(*)+#(*>> ( - ) • ’ •* •*** a ’ V ’ • j ( * Ж Ш flb • >• > Жл I • * ‘ e * -* * - — * *■ U • • • evidențiind faptul că răspunsul unui sistem se explicitează prin suprapunerea componentelor liberă și forțată Răspunsul cauzal la impuls, care intervine expres în definirea răspunsului forțat al sistemului, este susceptibil de o anumită interpretare care-i justifică denumirea primită Astfel, în cazul sistemelor netede, dacă notăm cu Д distribuția regulată asociată funcției hc (evident local integrabilă), rezultă că are suportul în [ , -Ș- oo) și deci hc e ф'+ care este o algebră cu convoluție Atunci ( ) devine în Ф'+ de forma j = h ж и cu yf e Ф+ ( ) Punînd = , distribuția lui Dirac, obținem direct din ( ) ca = & Deci, răspunsul cauzal la impuls, interpretat ca distribuție (regulata), răspunsul forțai al sistemului descris în ф\ prin { ), la o intrare de trpu* distribuției Dirac (fig ) ■■ f - - J~ i Г ’ со Г ‘ Se reamintește ca (/ » S} W A ț ~ = ) ю ?(т)/(Z “ т) - «o "K °® semnalelor continue, iar (/ » «)(/) « £ JW CaZ“ “т“а' Л) lelor discrete о Fig în cazul discret, interpretarea este imediata, impulsul unitar într-adevăr considerînd ( ) t = o t e { } ( ) furnizează, cu ( ), yf(t) — hc(t), ^ (vezi fig ) Fie acum spectrul matricei A, c(A) = {Xi, X , •••, Xn} și fie «=max{Re X^l R și deci l(f) Ф( » z— transforniabilă Notînd cu fE transformarea’ corespunzătoare, obținem din ( ) ( ) Ж{ (/) Ф( } ■= I - - г’Л’ - ••• ” Л)' în consecință, din ( ) obținem cu ( ) rWl) - - = z^cWO) -I- г>ф( ) -I- Г*Ф( ) - ] b «= Gr(xZ - Л)- b г I > R ( ) ( ) Introducem acum A Definiția Funcția de variabila complexă H(s) A ^hc(t) in cazul sistemelor netede și ( ) in cazul sistemelor discrete, se numește funcția de transfer asociata sistemului (Л, o, c ), considerat ca sistem neted și respectiv ca sistem discret Din ( ) și ( ) rezultă că pentru un sistem liniar (A, b, c "), indiferent de accepțiunea în care acesta este considerat, funcția de transfer are aceeași exprimare Я(Х) = ст(Х/ - A)- b ( ) unde prin X se înțelege o notație unificată pentru variabilele complexe s și respectiv z Explicitînd ( ) obținem Я(Х) = су(ХІ - A)+ b det(Xl — A) ( ) unde (ХІ — Л)+ este matricea formată cu complemenții algebrici ai transpusei matricii rezolvente ХІ — A Rezultă imediat că cT(\I — Afo este un polinom de grad cel mult n — Punem u(X)A^(x/-Л)+ еЖХ], [u(X)] Le ce are avantajul ехрІІсІИЛ ireductibila ( ), sa observam ■, nn(ic U(X) este polmomul minimal ÎT-Ж■ (s) Яр = s" - «„ îs"" - - «is - «o ( ) r(s) Jir — “ ••• ” *Ь % ( ) Atunci, pentru otice distribuție и е ф\ care, este un original, rezultă direct din aplicarea transformării Laplace convoluției ( ) că y(s) — —-— r(s) uts) — H(s) u(s) ( ) Ms) unde s-a ținut cont de ( ), ( ) și ( ) Așadar, funcția de transfer se asociază cu imaginea operațională a tranziției distribuționale descrise prin relațiile convoluționale ( ) sau ( ) Explicitînd ( ) cu ( ) și ( ; obținem interpretarea tranziției intrare/ieșire prin ecuația diferențială în distribuții (») (Я— ) («-!) у “ an-i У “ ••• ” ai ' * ѵ-ъУ = Pb~i u ( - ) considerate ca elemente ale algebrei cu convoluție Ф’+ și prin urmare avînd suportul în [ , -со) Așadar, interpretînd ( ) la nivel de funcții (interpretare evident abuzi ă), această ecuație trebuie privită ca o relație (diferențială) între funcți' definite pe EZ —cu suportul în [ , со) și în consecință perfect specificate pentru Z + Pi* + Pq y(z) zn +| -p atf + a ( ) Reamintim că aceasta este o dacă S+ este mulțimea funcțiilor algebră cu convoluție Fie p, q^S+ definite prin cu suport pozitiv Z*( = ^oW ~b an-i^o(^ ~ ••• “ ao^o(^ r(Z) = — ) » - powo(^ — n) unde ^ (£) a fost introdusă prin ( ) Fie ecuația în ( ) ( ) ( ) Cum p este cu suport marginii, p~r este bine definit în S±(p * p = Uq) Obținem astfel exprimarea (unică) у = (p" *r) * и ( ) a tranziției intrare/ieșire în algebra cu convoluție S+ Fie Z>M = ^W)} = + + - “ = z= z~n(zn + + ••• + «o)- ( - ) Atunci, pentru orice ( ), prin aplicarea ueS± z -transformabil, obținem cu ( ), ( S) și transformării ^ relației ( ), ( ) y(z) = — r(z) p(z) cuirea lui și r din ( ) și ( ) obținem v(t\ «- V(t - n -I- + - w) ₽ " + •” + “ rt) ( ) CU esenval tr-o ecuație cu diferențe ftmte Е„„рШ RetaS» «Iplot di" Л' PM “ *"*? de transfer este H(X și deci o exprimare ireductibilă este H(X) = polii și zerourile fiind >[H(X)] = { , - } și respectiv [H(X)J - { } în cazul neted, răspunsul cauzal la impuls este Лс(і) = l(t) стФ(і) = l(t) [ - e-‘] * - * V-j > * LJД • — > к de unde Xn cazul discret Deci Atunci «W = -=H(X)x Ecuațiile care descriu tranziția intrare/ieșire sînt în cazul neted și în cazul discret j/(Z) y(t — Dacă, pentru sistemul neted considerăm и « atunci, așa cum s-a văzut, у ry£ ! unde ; ( ) =o ( ) unde я-au considerat relațiile ( ) Soluția este șl deci -Я— b I тш „г ■Lkl I СИ О, X”) э ?( (Z) Același rezultat w u r ' " r " ''' t > """ “"■« I cT(Xl — A)~^ b = H(\)\ A A • V/ intii ca Să arătăm mai ^(д) » Pentru aceasta fin ^аъіа, uc #(X) exprimată FI(X) =; n se mmest& ( , ) în forma Atunci X«- -l~ Ș|X " po ■ ■ ’*■ ■ ■ OC[X (Xy I o o » o xo O tr o ~~aj o ao ( ) este о realizate a lui ( ) Fie Ѳ(Х) (ХІ - Я) Ѳ(Х) = о о Lao X о â= [IX Х" ]т Atunci "О О е о Lx(x) "О • Х(Х) — »Х(Х) ( ) -Ле у(Х)Ах» + а^Х-* + - + аіХ + care eSte caracteristic al matricei A, aceasta fiind de tip companion Din ( ) rezulta = (X - A)- b x(x) ( ) Cum cT (X) = ₽ + ₽iX + •;• + ₽п-іХв = ѵ(Х)г rezultă prin înmulțirea la stînga a ambilor membri ai lui ( ) cu c că H(X) = v(X) x(x) = cT(XJ-Л)'Ч ( ) și deci ( ) este o realizare a lui Я(Х) fiind în consecință un element al lui se face introduetnd noțiunea dc echivalență c, cT -b'ri Să observăm cft c(Ă) e [X]k Așadar, oricare ar fi Щ\), conține o infi* A Definiția Doua sisteme (A, b, cT) și (Â, b, cT) se numesc echivalente, notîndu-se (A, b, cT) ~ (A, b, cT), dacă există o matrice nesingulară T astfel incit ( ) Se verifică direct că ~este o relație de echivalență, adică este reflexivă simetrică și tranzitivă Se constată imediat că două sisteme echivalente au aceeași dimensiune De asemenea, două sisteme echivalente au aceeași funcție de transfer într-adevăr, fie (A, b, cT) (Â,b, cT), avînd funcțiile de transfer Я(Х) și respectiv Я(Х) Atunci utilizînd ( ) avem Я(Х) = c^-kl-Ây^b = c^T-^TT- - TAT~ )~ Tb = == cT(xZ — A) b = H(X) ( ) во ) Cu noțiunile introduse mai sus — asociată unui Я(Х) fixat Așadar, dacă (A, b, cT}e^ toate sistemele echivalente cu ( (Х) pentru ca sens • Teorema Doua polinoame /(X) -= Po "ІГ Ріл l " ^X) = a -I" « X -И - + “• Iх” * a"X“ ( ) ( ) • Să nu albă divizori comuni nelnversabih lu КРІ- I \ Д - к неіМіб nule, cu an* , sînt со prime dacă și numai dacă, următoarea nx X n m»itnce: • • * an- «o «O « ОС * * • ( - ) и ♦ P numită eliminantul polinoamelor r(V) și ^>(Ă), este nesingulară Demonstrație • * » Fie ЦХ) A[ XX X ”' ] Atunci din ( ), ( - ) și ( ) obținem egalitatea " PW "I W) X“" />(X) r(X) * X/(X) ! Х”-Ѵ(Х) M£(X) ( ) Condiția este necesară Fie așadar /(X), y>(X) coprime și să presupunem că M este singulară Atunci există l e ]Е§ в, I , astfel încît lTM == ( ) Punem Z == [йп й| >q by ( ) șj «(X) A a + «xX ’l- + «„ iX*’ ( ) fc(X) A- > + йхХ-F -Ь б^Х '" ( ) Înmulțind la sțînga ambii mșmbri ai egalității ( ) cu l'J' și ținînd cont de ( ), ( ) și ( ) obținem й(Х) />(X) + (X) r(X) ~ ( ) Deoarece r(X), />(X) sînt neidentic nule rezultă imediat că a(X), Z>(X) sînt de Д Г?(,JH'd п^елЬс ииІе Presupun» ni în ( ) că u(X) și Z>(X) sînt coprime, în caz contrar se Mmpldică divizorul comun Obținem din ( , ): r(X) «= — (X) Ж * * • ? : • • | ♦ ( ) I I I • c fel І M SÎnt c°Primc, rezultă că Cum a[/>(x)j ~ «(«„ o) și a[ (x)] (x) în ( ) obținem r(X) divizoi comun neinvcrsabil al lui />(X) și r(X) ceea ce este absurd întrucît r(^/(*) sînt copnme „ v ( ondiția este Deoarece M este nesingulară, există l e R " de (X)| />(X) și deci/>(X) = (X) înlocuind = «(X) > i ■ ’,K : ■ ■ înmulțind la stînga ambii membri ai lui ( ) cu lT dat de! ( ) și ținînd cont de ( ), ( ) și ( ) obținem ! v ' ** ' ' I li X) p{\) (X) /(X) (X) sînt coprime ‘ I >» ic H(X) = = r” ' ^(X) a„X« + X, (v și ( ), ( )) Atunci există R în serie Laurent I ••♦ | ’ * \ * astfel incit avem dezvoltarea ІП unde coeficienții y , Л se obțin prin identificare Obținem astfel ‘ Л** r Ш % Pq “ “ “ ” - «oY «»Y»+ sînt sau nu sînt !- ' coeficienții Y (X) sînt coprime dacă și numai dacă *'* * » rangX,n = « ( ) adică Л пп este nesingulară • ■ ’ * ' / • - Demonstrație Conform cu ( ) și ( ) avem următoarea factorizare a matricei obținută dintr-o permutare completă a coloanelor eliminantului ( ) al polinoamelor r(X) și £(X): и Ot д * ч * i Ya a % % - o В ( ) membri ai lui ( ) obținem Xnl - a,? ( ) numai dacă det ^/ Luînd modulul determinanților în ambii | det Д/ | = | det FSreZU\tă din ( J ) că del M daCa Și numai dacă det Utilizmd atunci teorema , demonstrația este complet terminată ■ ie (X) o uncție rațională proprie Se numește ordinul lui H(X) gradul numitorului oricărei fracții raționale ireductibile care explicitează pe Я(Х), Avem atunci următorul rezultat fundamental exprimat prin * ^eorema un șir de numere reale dat prin mulțimea Jfl din ( ) in care există cel puținjun i^ \ astfel incit у*/ Atunci există o funcție rațională pi opi ie ai cărei coeficienți Markov sînt exact elementele șirului Л dacă și numai dacă există un întreg n^ , astfel incit n = rang X+ rezultă automat din condițiile ( ) Condiția este suficientă Fie așadar șirul Л și condiția ( ) îndeplinită Atunci primele n relații din ( ) se scriu în forma ( ) Punem în ( ) aw (vezi ( ) cu i = j = ) a , an (monic) />(Х) «o + ,ai^ "I" " = (a„^ ) Atunci, deoarece este nesingulară a, л sînt unic determinați, deci polinomul W—A , X’*- Xм este unic determinat Mai mYlt datontăcondițieț(Y ), n coeficienții = determinați pnn(—r, /? qa\ (iP(hic unic coefi- • (X) sînt coprime, după cum re-sînt integral satisfăcute Mm mult r(X ч/И? i r zultă din ( ) și din teorema ordinul lui Я(Х) este deci tocmai « De precizat Că dacă л sa" {^гсЧап|'Ь)^гІтаШ prin; Teorema are o consecință iiducsaiH i XI t „ /•, A CU JC) ~~ ti- Corolarul I'uncția ]■ К „ dacă f numai daca condiția ?i are transformata z o funcție, raționala ({> P l и ( ) este îndeplinită^ ( ) satisfac relațiile ( ) pentru O n Л ' -a •«! / ^\ \ Z' Odată a,, Xе- fiind astfel după cum re Deci, r(X)//>(X) este ireductibilă, ia determinată transformabilă Fie acum ( , b, cT) un sistem Atunci avem (X/ — / ) = X Ч “ X I X | > R ( ) unde R Amax {| XJ | R ( ) Din ( ) rezultă că Yo=° = с ’Л ( ) sînt exact coeficienții Markov ai lui H(X) Din acest motiv, numerele , Л * , se numesc coeficienții Markov ai sistemului (A, b, cT) Re- i zultă că toate sistemele cu aceeași funcție de transfer au aceiași coeficienți Markov în consecință, revenind la familia a funcției raționale strict proprii H(X), rezultă că toate sistemele din CfHțX} au aceiași coeficienți Markov Fie acești coeficienți y , Y > ••• obținuți din Щ\) și unde sînt adevărate T • •• obținuți din H(\) și unde sînt adevărate egalitățile ( ), oricare ar fi (A, b, c ') Fie atunci (A, b, (?) arbitrar aleasă Avem, cu ( * > • \ " Yi Y ••• Y>- «« * •* A « , cT) Vom spune că realizare (Л» â, { a e ț/i, o, c ) e^ff|X) este minimală dacă orice altă / ' ?, r\ : - ' i* v ale dimensiunea mai mare sau reală cu a lui (Л, , cT) (adica cu dimensiunea matricei Л) S melfr din W deoarece dimensiunea tuturor siste- prin numere întregi și pozitive ?i»wus|«nca oricărei realizăm minimale din фп{ este egală cu ordinul var, fie H(X) reprezentată printr-o ireductibilă Dacă ordinul lui H(X) este n, evident [/>(X)] = n Fie de exemplu (Л £ Z realizarea controlabilă a fracției r(X)//»(X) și avînd evident dimensiunea Aceasta realizare este minimală Intr-adevăr dacă ar exista (A', b’, c'T) e CȚ-cu dimensiunea n' (X) nu este ireductibilă, ceea ce este absurd Rezultatul fundamental este dat de • Teorema Toate realizările minimale din Cfr Demonstrație n, n, atunci numitorul funcției de transfer al acestui sistem? n și deci r(X)/ Я(Л) sînt echivalente Demonstrație ■ Fie n ordinul lui H(X) Atunci toate realizările minimale au dimensiunea n Fie (Л, b, cT) și (Â, b, cT) dQuă realizări minimale Conform teoremei nesingulară Avem atunci din ( ) și ( matricea Hankel JLnn f " ‘ € * Г'ІТГТ rrtj\ f nn rezultînd că matricele sînt nesingulare Mai mult, din ( ) rezultă și deci Scriem ( ) în ^д(+ = ОЛ, ,, Dar, din, ( ) rezultă că / Д = [& | ~ l' Introducînd ( ) în ( ) obținem Șl ;(vi/ '(X) Vom proceda astfel: a Se calculează coeficienții Markov Yr, Y > •—> Y »-> unde d[/>'(X)] (evident y = deoarece (X) este strict proprie) Calculul se face utilizînd recurent relațiile ( ) și ( ) b Se formează matricea Y • • • Y»- Уз • • • Yn'+i Yn'+l • ’ ’ Y n'- și i se determină rangul Fie rang JCn'n^ Evident n^n' c, Se construiește matricea JCnn care obligator este nesingulară și se rezolvă sistemul ( ) cu an — Este astfel determinat numitorul mcmic /(X) =â- Xя + + - + a al ireductibilei d, Utilizînd ( ) se determină coeficienții ( , , ( n l ai numărătorului ireductibilei r(X) =^= РЯг Хп -f- -f- pxX -|- ( Din r(X)/y>(X) se construiește imediat o realizare minimală (vezi punctul de mai sus) d' Realizarea minimală se poate construi și direct fără a mai efectua punctul d Ea este (vezi ( )) f • • • ♦ • * t • o ~« — a ♦ ♦ — a cî> — [Yi Ya ••• Y»] ( ) o o sau dual A (X) = X + X - d Utilizînd ( ) obținem £ =l, pr = Deci r(X) — și Realizarea controlabilă este atunci d' Fără a mai utiliza ( ) obținem • Г - din ( ) direct г = [o ij L o Să considerăm acum un șir y — , Yi> Ѵг> unde exista cel puțin un Ol astfel îneît Y , daca și numai daca condiția ( ) este îndeplinită) Or, condiția ( ) este m general greu de probat întrucît dacă șirul Y / O nu are o regulă de generare vizibila, este greu de probat că, condiția/ ) are loc Vt,;> a presupunem însă că există un astfel îneît rang « Anlicînrl пячіі a b c si d de la punctul , cu referire la у o avîn/ordinul» și m cărei coetaenf Markov comod “ ( ,»;«') “ M /'(X) vom W uuuri Yi = ( ) , z л I, / f} este o realiz(lTf‘ Рат(М a șirului Y Vom spune că sistemul /Л M acestui sistem stnf egali cu termenii Y astfel incit Ф(/) ( ) asimptotic stabil dacă г ф(/) —► CU t ' - Sc 'constlW ^ ^/ - " eS } ( ) atunci condiția b de stabilitate asimptotică din teorema de mai sus este echivalentă cu unde а(Л) este spectrul matricei A • în cazul sistemelor discrete avem din ( ) = Mpjp,- + ^p,Y₽,-i(xr, + ^p,-iYpr- (Xr, Z)X* + + — unde t max {p | astfel încît w>o ( S) b strict stabil (extern sau strict stabil în sens intrare măr^inităHeșire măr-ginită) dacă ( ) în cazul sistemelor netede și, ( ) in cazul sistemelor discrete Fie Н(Х)Асг(ХІ — Л)' b funcția de transfer a, sistemului pe care o explicităm printr-o fracție rațională ireductibilă //(X) r(X)//>(X) cu y>(X) monic Punem / ' • PW — (Ъ ~~ Рі}Яг Pa) '' Pi) ' ( ) unde pi> p , >-> Pi sînt rădăcinile distinct^ ale polindinuliii p(),j Datorită faptului că r(X)/'/>(X) este ireductibilă polii funcției raționale //(X) sini \ ,{Ръ Pii ••••» ,/’i> •••■» Pi> Ph "Ч Pi}‘ ( ) ' Чу 't In consecință avem descompunerea U #(X) = (X — pi)' ' c, ", | a'i ■— imr ■ ■ ■■■ " ■■ » ■■ ii, '■ -(x ^piYx** ( ) W >иэО*Д t A ь\ь ич Аж чж Obținem atunci; • Teorema , Un sistem neted (Я, b, c' ') este,', a extern stabil dacă ■>/ numai dacă toți polii /itncției dc transfer ll(s) au partea reală w [О, оо), rezultă că -» г Reciproc Dacă Re O, i O cu t —► oo și deci obligator Re pi /) -л Ѵ// ( ) rezultînd din ( ) că z - *• * = l> Atunci, cu o tehnică de demonstrație similară celei din teorema obținem din ( ) Teorema Un sistem discret (A b, er) esle: ' a extern stabil dacă și nwnm daca Щ' b extern strict stabil dacă și numai ducă ta(i palia famțaea d tram hr rtV ) au modulul astfel încît oricare ar fi comanda uf} cu u{t) = p&târu t , răspunsul forțat corespunzător yf ) esre rrdz-ginit, mai precis, | yff) | Demonstrație Abordăm cazul sistemelor netede urmînd ca cititorul, într-o manieră similară, să procedeze la demonstrarea teoremei în cazul discret Condiția este necesară Din stabilitatea strictă rezultă că • A există un T> astfel Z strict stabil Avem atunci ( ) N există T> și comanda ( ) astfel îneît А contradicția dintre ( ) și ( ) încheind demonstiația ( ) ( ) ( ) Din dezvoltările de mai sus rezultă în final că, problema stabilității este o (sau caracteristic) al matricei A, în căzu stabll%t^^eă ftl etia de'trar^er ( ), ( )) este înțr-o mutțime dedusă ^dăcinilor înC~ Ui( ) în cazul sistemelor discrete) Cu leteim Cm z m; nncesară si suficientă ca rădăcinile poli- Criteriul lui Hurwitz Condiția necesaia , nomului de la numitorul unei i> H(X), în ca și în nomului PW = aMXM * a»-i^ ' + + a ' ( ) să fie toate plasate în ■ este ca minorii principali ai wxn matricei I t - ' c i ' t x - • V «» o o o “n o o Q * ( ) i ) I К f ГіЧГП f I I ГГr ~ i > Referitor la locația rădăcinilor unui polinom în U^O) se utilizează de regulă o transformare conformă a lui U^O) în [Я(Х)] în particular, implicația inversă implicației ( ) are loc în cazul cînd forma primară a funcției de transfer este ireductibilă, ia care caz incluziunea ( ) devine egalitate De remarcat că dacă ^[Я(Х)] c: ’ u »‘ recunoscînd imediat că sistemul este reălizareă observabilă acestei funcții de transfer, Să considerăm sistemul (Atb,c^) neted Atunci aplicînd criteriul lui Hurwitz polinoir ului v(X) obținem АіР і o cu minorii principali Д, = До — Д тиШТ/Ь ti іоілі ig: îaqhatz зяаятиі за ’ Sistemul nu este deci (intern) asimptotic stabil^; gj эппС{ ШО f£Tg£ £Cj * rezultă ponform evului ІФ* iul doi de la numitor sini Cum toți coeficienții trinomului do g-a și deci sistemul este asimptotic stabil Rezultă direct Că sistemul este și strict RĂSPUNSUL DE SISTEMELOR LINIARE LA MĂRIMI INTRARE STANDARD în acest paragraf vom pune în evidență unele particularități ale răspun-sului forțat al sistemelor liniare (v ( )) la așa-numitele mărimi de intrare standard Vom înțelege prin intrări standard două categorii de intrări: — intrări polinomiale de tipul u(t) = !(/) Zu- în cazul sistemelor netede și, u(t) = - ( , /), u> ( ) ( ) u> in cazul sistemelor discrete; — intrări armonice de tipul гірі) = ejto/ (/) (m ) ( ) cu t g]E£ pentru sistemele netede și t e pentru sistemele discrete în cazul intrărilor polinomiale avem atunci și respectiv I IV • Pentru sistemele netede avem ( ) ( ) ( ) >z(x) = //(s) Facem ipoteza că s nu este un pol al lui H(&j Atunci и—l Z/(s) • = / ( ) -L -h | , p W- + f/(s) ( ) s ! su~ (и — )! s unde fî\s) este o rațională strict proprie originalul lui ( ) este /и— ^(f) = H( ) cu polii identici cu ai lui H(s) Atunci H'( ) ( ) unde уМ А \ Я( ) ( ) și obținem л( — Ур^) + Уі( ) • Pentru sistemele discrete avem ( ) yĂz) = H (*) ( ) • «• v Cu ipoteza că z — nu este un pol al lui H(z) obținem Я'( ) • - = H( ) - u (г — H(l) + Я(г) ( ) Щ$)- Punem + y,(t), unde H(z) este o rațională strict proprie cu polii identici cu ai lui Я( ) diu ( ) avem yt{t) = со , iar consta- сц menținci^■■ipotezei de stabilitate strictă, răspunsul forțat al sistemului se mai numește și răspuns indicial, iar componenta permanentă este ’ ’>„’(/) ^=H( j, /> , ( ) în cazul sistemelor netede și, №(/)==//(l),v/> ( ) ■ în cazul sistemelor discrete Întrucît aceasta componentă este constantă în timp, ea se mai numește și componenta staționară notată adecvat y* у r(tf Așadar, dacă sistemul este extern strict stabil H( ) și H(l) pot fi interpretați ca factori de amplificare, primul asociat sistemului neted, iar cel de-al doilea sistepiului discret, ț , л în cazul intrărilor armonice avem > • ț L- ( ) x ' Î s — j O) și respectiv - -ă ■ ■ - ; v ej • ‘ fU? i î З'/О = Ж V/(z) г Я(е ix К V ( ^ ) și ( ) capătă forma ( ) In ipoteza stabilității (ex-terne) stricte ,a sistemului ‘“* cu / oo devenind o componentă tranzitorie \ (a răspunsului forțat) în acest context y^(t) definită prin ( ) devine componenta permanentă (a răspunsului forțat) Atunci din ( ), cît și din ( ), rezultă că în regim permanent, intrarea armonică provoacă o ieșire , de asemenea armonică, modificată însă ca amplitudine și defazaj £ , î D v» Ca o concluzie generală a celor prezentate mai sus se evidențiează faptul că? dacă sistemul este (extern) strict stabil, atunci în cazul intrărilor polinomiale sau al celor armonice, răspunstil forțat al sistemului se desface în două corn- ^н(о) de același lip cu intrarea H(s) = r O | ~ t\ numit și element de ordinul Să studiein comportai — Pentru intrarea în treapta unitară avem к «A fK o л ' г—' I ’ r а acestui element la intrările standard //(()) i I И I J Se constată imediat că yp(l) Uryft)^- -Ко este componenta Ігаи/ДІопе > l b > (), este componenta permanentă \ v De asemenea ( ) T ( ) interpretîndu-se mărimea Г ca — Pentru intrarea în rampă o „durată a regimului tranzitoriu", unitară (u(t) — !(/)/) avem ( = z ( ) Punem tga=A=T ( ) partea reală a răspunsului ( ) corespunzînd intrării ( ) cos со/, iar partea imaginară corespunzînd intrării l(/)sin ( ) Avem de asemenea — » V l — - - / - > ț dl Din ( ) ) rezultă că extremele faliei » momentele înlocuind ( ) în ( ) obținem >z(/t) = ( ) cf* ZT « И ;л° Fig £ ІпшЬіо эЬ ідо/пэіэ ц Нспяп • Îîtqn eij rn “ (o cr; rk'! sau ѴЖа w Cu ( ) și ( ) obținem pentru adică apare un regim oscilant ayînd pulsația egală cu cea naturală; — pentru £ =* , așa~n urni tul regim aperiodic critict care de asemenea rezultă direct din ( , ) ca fiind yz(/) = — e Ambele regimuri limită sînt ilustrate în fig, ( ") Din locul (v ( ) și fig ) rezultă pentru regi mul tranzitoriu, conform cu ( ), că geometric al extremelor ч\ > I ; ( ) " r ; »( cu cea a elementu-I - — ’S’ , ?- Fig ț \ numit și element de ordinul cu un zero Avem atunci n n Punem tg a — C și obținem уАѴ = unde p Un calcul imediat arată că d^(Q d/ /=o Polii y>j și zeroul — z sînt ilustrați în fig , unde -p ( ) unde К = — ѵій! + я Și deci H(l) = - Pentru treapta unitară l(t) răs- sin w/ Punem cos w A tj Atunci ( ) devine COS * + J â ( l - v)«)/Vl obținem din sin (м/-F a), /> sin a У/Ѵ) Introducînd în continuare ctg a У/W » ctniiască o paralelă cu cazul elementului Cititorul este invitat si consternase» neted de ordinul I $ ' Г\ ri Un exemplu tipic în H ț Un exemplu tipic în care nu are loc descompunerea răspunsului forțat în regim tranzitoriu și permanent este cazul așa-numitelor elemente cu caracter / л л / ) л лл Л rf/і/ Г\ЛГ (multi) integrat or i ■ ! ж V s ( , ) j sau Я(~) \ G XI г ,n , ((ШЛ) + l> ]E£ în funcția pentru kh^t O^iar ГІС“ CAN efectuează conversia funcției g: Й prin g(k) = "^(kh), k'G’Z ’ Sistemul ( ) se numește discretizatul cu pasul - h(-h > ) al sistemului neted ( ) semnificația operației de discretizare fiind ilustrată In figura , explicitează modalitatea de fe) oW un sistem, care natural tranzi- - - —~ J! tranziteze -funcții definite pe I ' * în funcția g: care v tează funcții definite pe ЗЕ^ , să I / I Exemplul \ t ) oc | л I A г o b = - j T гпиэ hnttJA Să calculăm discretizațu) cu pasul h Avem (exercițiu pcntiu cititor) ~ •) ' ” ■ ' ™ —n - /J) Obținem cu eJ ( 'ЛІ r\ —I ( ) - Bl prin fd(k) = =f(kh), keZ și h> arbitrar ales, numită discretizata cu pasul h a funcției f Cum /( \^Meat cu obținem ІЖ I = l/(^) I К* I > e** și deci ( ) definește transformata z a funcției //£)’, notată Fd(z) Obținem astfel din ( ) ВД =/( ) +>) + /( Ă) г~ Dar c> a ( ) unde F(s) = J>f(t)- înlocuind ( ) în ( ) rezultă ( ) Atunci cum Res = c > a, luînd | z | > есЛ > > °^ rezultă \esh/z | a, | z |> eeft ( ) Formula ( ) explicitează legătura dintre transformata Laplace a funcției original/și transformata z a discre tiza tei cu pasul h Vom prescurta ( ) prin notația ( - iZ[F(s)-] ( ) ech revine la a pune z = eph și a lua Hep> c (v figura ) Dacă F(s) este strict proprie, atunci integrala ( ) pe arcele de rază infinită С\ sau C este nulă în consecință, aplicînd teorema reziduurilor, cu referire la conturul care conține polii lui F(s) obținem din ( ) sh i=l L * formulă deosebit de utilă în calculele efective Exemplul Fie f(i) = l(^) e i{(; J/ib'^nto !| I J '• Іп'шйюэ я! Atunci r \ w este p ri(X) = (X — ех>л)рі • (X /г ex»z,)p« (X — e^*)* Aceeași proprietate de conservare o arc și stabilitatea externă Minimalitatea nu este însă o proprietate care se conservă întotdeauna, din simplul motiv că H(s) și nu se exprimă în mod obligator prin ireductibile de același ordin Mai exact, dacă (A, b, c ') este un sistem minimal, adică forma primară a funcției sale de transfer v(ș)/x(s) este ireductibilă, iar (cî) discreti- zatul său cu pasul h (v ( )), atunci ireductibilitatea formei primare a funcției de transfer discrcțizate Яй(я, Л) — cq\ș I — A^h)” ) ba(h) == vd(z, h)/ l^(z, h) va depinde de ht adică de acele valori pozitive ale lui h care anulează determinantul matricei climinanț, a polinoamelor u(f(^, Л) și ^rt(^, Л)* Or* acest determinant este evident o funcție analitică de h și în consecință st? va anula pentru o mulțime cel mult număr abilă de valori izolate ale lui h Așadar, dacă (zi, b, c ) este un sistem minimal, discretizatul va fi de asemenea minimal, pentru aproape toate valorile lui h Deci, conservarea mini-malității prin discrctizare este o proprietate generică REPREZENTAREA (ÎN FRECVENȚĂ) A FUNCȚIILOR DE TRANSFER Așa cum s-a văzut în § , răspunsul permanent al unui istern la o i ntrare armonică cu care sistemul nu rezonează, este у fit) = HI fi,) în cazul neted și у fit) = e^ în cazul discret Se numește reprezentare (în frecvență sau mai corect în pulsație, a funcției de transfer, o reprezentare grafică a dependenței numărului complex H(fii) și respectiv Я(е’“) de ) = | H(j ) A arg Я(» ( ) care definesc dependența modulului și respectiv fazei (argumentului de pulsația со Cu ( ) și ( ), ( ) se rescrie prescurtat H(jco) = Я(со)е^) ( ) De asemenea, punînd U(a) AReE(je) ( F(o) AlmH(jG)) ( evidențiem partea reala și partea imaginară a funcției de translez Așadar H(jco) = Щи) + j (*>) ( ) Să observăm că deoarece Я($) este o rațională peste ® avem “Я(У) = #($) ( ) unde prin suprabarare am notat conjugatul Atunci din ( ) cu ( -xj rezultă ZZ(joj) == £ (g>) + j= AT(jco) = “■ Л (“’ / ( ) Așadar, pentru weJES , unde Z/(jw) este definita, avem £ (g>) == D(—со) ( »оэ) și ^(co) —V(— со) ( ) rezultînd că partea reală este o funcție pară de w, iar partea imaginară este o funcție impară de o> De observat că dacă « = este un punct oe continuitate pentru P(o>), obținem din ( ) ca ( ) = ( ) Cum Z («) = [(/ ( )] / rezultă că Я( are exprimarea în dB, (X)dB = log x, unde log este logaritmul zecimal * •v ш * • • * * * ** • • * ы • > t • Caracteristica amplitudine-fază este reprezentarea în planul (Я, ) a dependenței parametrice H = și = ( ), unde abscisa este gradată în radiani (grade), iar ordonata H în dB Avem așadar următoarele planuri de reprezentare ilustrate în fig a, } Cj d e în cazul sistemelor discrete să observăm mai întîi periodicitatea funcției Я(е]'“) cu perioada îc în rest se introduc exact aceleași tipuri de reprezentări ca și în cazul sistemelor netede, utilizîndu-se notații identice Astfel Я(е^) = Я(со)е^“) = Я(со) + jF(co) ( ) unde Я(со) = IЯ(е]ы) |, ( ) =A= arg Я(е]’°), Я(со) = Re Я(е^), P(co) = = ІтЯ(е^) Relațiile ( ) și ( ) fiind valabile, intervalul efectiv de variație al pulsației va fi со e [ , тс] (similar cu со € [ , oo) în cazul neted) Ilustrăm în continuare, pentru cîteva exemple remarcabile, aspectul caracteristicilor menționate mai sus Precizăm însă că, eu toate că funcția de transfer Я($) sau H(z) a unui sistem neted sau respectiv discret esteorațio- v ddk H(cS) b nală strict proprie, în cele ce urmează vom ignora din punct de vedere al reprezentărilor această restricție și vom înțelege prin (s) sau Я(г) funcții “J “ de exemplu H(s) = s«, q>\ nu are nici’o raționale arbitrare, chiar dacă semnificație sistemică Fio clementul multiintcgrator К -r, — > iL > Cum H(jco) = • к и Pentru caracteristicile logaritmice avem H(co) = si deci [H(co)]dB = [Щав • q log co radată în scară logaritmică, ( ) reprezintă o dreapta care Cum abscisa este gradată în scară logaritmică, ( ) reprezintă o dreaptă care trece ?nn punctul ( , №) și are panta de - q dB/decadă deoarece de la O obținem somidroapta din figura /, b, care inversată față do cercul unitate se transformă în semicercul din fig, , a Pentru trasarea caracteristicii amplitudino-pulsație obținem din ( ) *= ~ W + ȚW) ( ) A/om atunci °> Ю Ig = ( ) ?i — =s> [ /( s = - Igîb ( ) care este o dreaptă care taie abscisa în = — și are panta de — dB/decadă Pentru II Din ( ), ( ) și ( ) rezultă următoarea alură a caracteristicii amplitudine-pul-sație (v fig ) Caracteristica [ - j(co)]dB — ( ) Avem FT(jco) ( ) unde л' = со/соп se numește pulsația normată Obținem din ( ) V(«) = ІП— -И(ж) = ( ) ( -л- ) + ^« ( - ж ) + £ * care sînt ecuațiile parametrice ale locului de transfer ilustrat în figura cu evidențierea dependenței de factorul de amortizare £ De remarcat că pentru , n) => [Н(*)ЬВ = - Ю lg = ( ) x > (w > O n) => E (*)]dB = x ( ) x=l (оз = wn) => [H( v)jdB = — lg £ ( ) De asemenea, extremele se obțin din anularea derivatei funcției ( ) care are loc pentru # = și pentru x = sj — £a iar valoarea (maximă) corespunzătoare lui ( ) este - - ig W-W care este valoarea amplitudinii de rezonanță Din ( ) și ( ) se deduce imediat caracteristica asimptotică [ялМ](ів ~ , — Ig x, ( ) fiind formată din două semidrepte: prima confundată cu axa absciselor pînă la x = și a doua plecînd din punctul x~ de pe axa absciselor și avînd panta de — dB/decadă Din ( ) rezultă că rezonanța apare pentru —r = , Din V ( ) se pot deduce valorile caracte risticii exacte la pulsația de frîngere x = (со = con) a caracteristicii asim ptotice Cele menționate mai sus se ilustrează grafic în fig în care sînt trasate caracteristicile amplitu-dine-pulsație cu evidențierea dependenței de factorul de amortizare Pentru trasarea caracteristicii fază-pulsație obținem din ( ) Z-OJ Гй; ср(лг) = — arctg ( ) — x care este ilustrată în figura De observat că la pulsația de frîngere a caracteristicii asimptotice x — l(co = con) defazajul introdus este — î/ , iar cel total este — тс în cazul rezonanței, £= , are loc „saltul de Fig fază" egal cu — Caracteristica amplitudine-fază parametrizată cu factorul de amortizare rkOdBjdec este redată în fig Pentru elementul „invers" H" »^) = - * + j toate caracteristicile se obțin prin simetrizarea caracteristicilor* din figurile , și Trecem acum la sistemele discrete și considerăm - J • ' л , ' • Elementul (multi) integrator discret: H(z) = -*> ( ) ( - )« Rezultă / x л ( со i / со I cos q — — sin I —• Ч | К = K l ț J ( ) Я(е]“) = " a jesin pH ‘ ( ) — i ( ) Introducem următoarele noțiuni: — excesul polilor față de zerouri e A d[p(s)] - [r(s)] > *>; ( ) — factorul de amplificare al sistemului este constanta К (vezi ); — tipul funcției de transfer este întregul (pozitiv sau negativ) q care intervine în ( ) Datorită condiției ( ) avem f autorizarea în factori ireductibili Hr - nV r(s} - П VA J(s) = U ( , ) unde r£(s) și pf^s) sînt polinoame de gradul întîi ayînd rădăcini-reale, sau polinoame de gradul avînd rădăcini complex conjugate în consecință, /г-($) și pfs) se reprezintă sub forma unor /actori liniari unde T,it Tpi sînt de semn arbitrar și au dimensiunea unui timp nummdu-se și constante de timp, sau sub forma unor factori pătratici' Г ■*- npi \ tri, numindu-se factori de amortizare, iar Tnri, TnJ)i constante de timp naturale Cititorul va recunoaște în ( ) forma standard a funcției de transfer inversate a elementului de oi dinul , constanta de timp avînd însă arbitrar De asemenea, introducînd pulsațiile naturale un semn ПРІ rp ■t npi nri x nri ( ) ( ) se rescrie în forma piwnpi npi ( ) = f I ) * ha^) = k °-+ ^ ^ + L(jco)« ] (s) respectiv Se vede atunci că, comportarea este dată (pentru p —> ) de către ( ) rezultînd că modulul tinde la zero (fapt evident deoarece H(s) este strict proprie), iar faza tinde către ( ) (v- fig , a și &) Detaliile intermediare ale hodografului se obțin prin reprezentare parametrică cu ajutorul, de exemplu, al unui tabel de variație I со- Exemplul , Fie dc unde sau este asimptotă — СО Pentru avem este dezvoltarea ■ rea ln jurul lui Н(іы) = C” tabeIuI (șimplificatj Comportarea peBtru Kczultă locul din figur e\ Exemplul p;e este: "(s) = = ? ) jco +- jw) o- — pi Vom presupune că în reprezentarea ( ) sînt îndeplinite condițiile uzuale: К, Tri> T Qpi -, —— > -; ri Pi nri npi de exemplu în cazul ( ) rl npk -* r nrl jpl rezultă pentru caracteristicile asimptotice corespunzătoare fiecărui termen din suma ( ) dispunerea din figura Caracteristica asimptotică globală, ilustrată în figura , va fi suma caracteristicilor fiecărui termen, rezultînd o linie poligonală cu variații de pantă la fiecare pulsație de frîngere Prima dreaptă, cu panta de — ț? dB/dec = — dB/dec (s-a ales spre exemplificare q — ) se numește asimptotă de joasă frecvență, iar ultima dreaptă cu panta — e dJB/dec, unde e este excesul poli/zerouri, se numește asimptotă de înaltă frecvență Pulsația co, se numește pulsație de tăiere kdB ■ [H(u dB Fig Pentru, trasarea caracteristicii exacte va trebui să se ia în considerație faptul că un termen liniar din suma ( ) introduce la pulsația de frîngere o abatere față de caracteristica asimptotică de — dB dacă este negativ și de + dB dacă este pozitiv, adică după modul cum este situat Ia numitorul și respectiv la numărătorul funcției de transfer (v fig ) Abaterea introdusă de un termen pătratic depinde de factorul de amortizare £ (v fig ) care condiționează abateri atît pozitive cît și negative față de caracteristica asimptotică Cu aceste precizări, calitativ, caracteristica amplitudine-pulsație va avea alura din figura / De precizat că abaterile de ± dB puse în evidență în figura sînt ipo~ Я tetice întrucît termenii se interinfluențcază (grafic) conducînd la variații ale ж abaterilor standard de la caracteristica asimptotică ale fiecărui termen con- siderat ca atare Evident, cu cit pulsațiile de frîngere sînt mai depărtate între r ele (cel puțin mai mari ca o octavă), abaterile caracteristicii exacte de Ia caracteristica asimptotică vor fi mai apropiate de abaterile individuale ale j termenilor care compun suma ( ) Din cele de mai sus rezultă următoarele reguli dc trasare a caracteristicii ( amplitudine-pulsație în condițiile uzuale precizate la început Se aduce H(s) la forma standard ( ) Se calculează Ig К și se figurează punctul de coordonate ’ ( , [K]dB) Pe axa со, gradată logaritmic, se dispun în ordine crescătoare pulsațiile de frîngere /Tn, /Трг, /Tnrt, \!Tnpi, ridieîndu-se verticale în aceste puncte Prin punctul de coordonate ( , Г?К]йв) se trasează asimptota de joasă frecvență de panta — q dB/dec La intersecția asimptotei de joasă frecvență cu prima verticală se modifică panta acesteia cu ± dB/dec sau ± dB/dec, după cum verticala corespunde pulsației de frîngere a unui factor liniar sau pătratic (plasat la numărător sau la numitor) La intersecția dreptei nou obținute cu următoarea verticală se va proceda similar Pentru trasarea calitativă a caracteristicii exacte se introduc în punctele de frîngere ale caracteristicii asimptotice corecțiile corespunzătoare factorilor liniari sau pătratici Evident, un calcul în puncte intermediare mărește precizia de trasare în exemplele care urmează vom parcurge etapele de mai sus: Avem: „ > + s = — • : • $ ( + IGs) ( - $) [ ]dB = g = dB — = , ) = - — ■ cu reprezenta rea din figura în fine, pentru exemplul , caracteristica fazâ-pulsație este trasată în figura Sa presupunem acum că, infirmate, adică К ilustrată în figura Fig Exemplul H(s) = - s — , s - , sa Se vede imediat сП IC ® — > £ = — , , precum și Tw — , și q Atunci î?(s) = - , s - , sa pentru care corespunde caracteristica amplitudine-pulsație din figura Caracteristica fază-pulsație rezultă din figura CORELAȚIA DINTRE MODULUL Șl ARGUMENTUL FUNCȚIEI DE TRANSFER > Dezvoltările care urmează se referă cu precădere la sistemele netede, dar pot fi adaptate și la sistemele discrete, operație lăsată în sarcina cititorului Stabilirea corelației dintre modulul și argumentul unei funcții de transfer se încadrează în așa-numita teorie a lui Bode, teorie încadrată în analizași sinteza circuitelor electrice Fie o funcție de variabilă complexă care satisface următoarele condiții; , ~ H(s) Este analitică în Re s > ?), Este analitică în jurul oricărui punct de pe axa imaginară, exceptînd un număr finit de puncte jefy, pentru care însă ($ — jo)J H(s) —► cu s jco i = , / ( ) lim —~ =» uniform cu | args | -JP> iw + ÎPJU и U [j° те Jo « — ? ( ) se rescrie cu ( ) în forma o f“ C (to) - C (w) , (ă>) = v p -\ -Ц -=^-dco те Jo со — со- care pune în evidență o exprimare integrală a părții imaginare în partea reală a funcției -fl(s) Facem acum schimbarea de variabilă P , satisfac condițiile , , și menționate mai sus Astfel de funcții de transfer se numesc de fază minimă Este ușor de constatat că o funcție de transfer este de fază minimă clacă și numai dacă îndeplinește condițiile uzuale de trasare a caracteristicilor semilogaritmice în cele ce urmează în acest paragraf ne vom referi exclusiv la funcțiile de transf er de fază minimă Atunci scriind In H(jcc) == п + j?(°) ( ) (unde Я(со) ^ |Н(] бз)|, ф(й)) =A= arg formulele ( ) și ( ) devin ) j do> J Obținem astfel un rezultat remarcabil exprimat de: ( ) ® Teorema (prima teoremă a Iui Bode ) Pentru orice funcție a: transfer de fază minimă, argumentul acesteia are exprimarea integrală г-u junc-Lie de modulul aceleiași funcții, dată prin formulele (echivalente) ( ) și ( ) Ca o aplicație imediată a formulei ( ), să evaluăm comportarea fazei unui filtru trece jos în banda de trecere a acestuia Cum Я(со) ~ct în banda de trecere [ , rezultă direct că о(оз) o> TC in g(o>)№) CO — CtF ( ) [ , coj Cum >co obținem în continuare In Я(б>)/Я(сЪ) ( ) Punem ) = — și ( , ) devine ф(б)) Ci ( ') • Y'ZU “Г ’>J’ а, pentru o» e ) d ) ( ) adică defazajul introdus este un multiplu de n/ , egal cu panta caracteristicii Considerăm acum caracteristica ( I din figura unde c=~ * (O Avem atunci (v fig ) a, pentru cog[ , coj, I j fco tp( j) = — \ (—?) ia cth— d = Г = ~#(U ( ) unde ( ) b pentru cog[coc, ], E Ș* deci ф(со) == —q—+ ( ) ( ) Rezumînd, obținem din ( ) și ( ) I—#(£c)> « e “J fee = ln ~ > °) \ / ■ ' • v • —+?/(-U» «e[oc) +oc) L=ln^ fii" + U + — t (— ? + ?) ln d£ = — + Л L + —(°° (—qr + q) In cth —( (—? + q) ln cth-=- dE, = * J- , тс Jț = -? — + (-qi + + (-? + ( ) ^^ИИІИМмИИИВММИИ Facem ipoteza că со și f) А? ( ) Din formula ( ) rezultă că dacă & q = atunci o(c^) = — • Așadar, daca pulsațiile de frîngere sînt simetrice, iar panta intermediara este mea ia aritmetica a pantelor extreme, atunci defazajul la pulsația de taiere cot este impus exclusiv de panta intermediara, fiind egal cit —q Din interpretarea formulei ( ) precum și din exemplele analizate mai sus rezultă următoarele concluzii generalizatoare referitoare la faza funcțiilor Fig de transfer de fază minima: IEI Cum funcția ln cth —- este „rapid descrescătoare" (v fig ) se poate aprecia că local faza la o pul saț if o este „în principal" determinata de panta caracteristicii amplitudiiie-pulsa-ție la acea pulsație (v primul termen din ) Defazajul total (adică pentru ulsat™ esU exclusiv determinată de panta asociată acesteia Pentru funcțiile de transfer de fază minimă sînt interesante două probleme: niodulul Я(со), să se determine H(s) de fază minimă astfel încît I Я(]со) I == Я( ) Я(—jco) = Я(]со) Я(]со) = Я ( ) astfel încît arg H (jw) — (w) COS ) — se efectuează factorizarea (s) — H(s) rțs) j>(-s) H(—s) p(s) r(—s) ( ) unde H(s) nu are poli și zerouri în к л« " ~ J j Й к ’fэ > Exemplul Fie tg Cp(w) = —со Atunci l-jb>s + (j ) > , Vw e K ), unde Л(і = , : , , , deoarece TlJs) este o funcție de ș Cum Hu bi) Iw =* ii ** ’(s) I» -« l)' ■’ ( ) ( ) ( ) * V de exemplu G li C a r t i a u u Анаіім f infam circuitelor electrice Editura Di- dactică și Pedagogică București, J ■ condițiile ( ) revin la (v ( ), ( ), ( ), ( )) гл = c» = = Со» = care prin identificare din ( ) implică On = ^ = — ^ n-l ( ) Cu ( ), ( Э ) devine unde b > C de normare Alegînd din considerente I Ѵ я i — \ / și inbodurfnd vâri bfia normată s/*>a, cu păstrarea abuzivă a notației, obținem din ( ) ( ) Factorizînd ( ) în raport cu axa imaginară vom obține din ( ) pe IIa(s) Polii lui Ha(s) sînt rădăcinile ecuației + (— )'* $‘;i = adică n e Л л * ( ) a căror dispunere pe cercul unitate formează așa-numita repartiție de tip Bvtterworth, Cum Ha& este strict stabi]ă, ea va fi perfect determinată de polii din O" ai repartiției Buttenvorth, care sînt cu reindexarea lui ( ) ( k -m n n în definitiv cu pk furnizați de ( ) și unde în final s se înlocuiește cu s/ob Exemplul Pentru n = avem și ceci sau după înlocuirea Iui cu O fiind un element de ordinul cu £ = = , ș’ ~ соь S Metoda expusă mai sus se mai numește firesc și criteriul modulului Există și o metodă duală numită criteriul argumentului; de asemenea, în locul aproximației de tip Taylor se poate utiliza o aproximație de tip Cebîșev,* DETERMINAREA RĂSPUNSULUI INDICIAL DIN REPREZENTAREA ÎN PULSAȚIE A FUNCȚIEI DE TRANSFER Ffg, , Vom fi interesați în stabilirea răspunsului indiciul al sistemelor neted' și a proprietăților acestuia utilizînd diferite modalități de reprezentare în pulsație ale funcției de transfer Este studiat de asemenea răspunsul indicial al sistemelor discrete prin „reductibilitatc" la cazul sistemelor netede» Fie în acest scop Я(з), ofuncțiede transfer strict stabilă, adică Atunci, așa cum am văzut în § , răspunsul indicial se obține pentru u(i) = (/) și decij(s) = H(s)/s Utilizînd atunci formula de inversiune obținem • y(i) =—ț+) ®-c’(ds, /> , Vc>o ( ) п] s (egalitatea ( ) fiind adevărată și în t = întrucît //( s)/ s este o rațională cu excesul e> ) Deoarece H(s)/s este olomorfă în semiplanul Re s> , integrala ( ) are aceeași valoare și pentru conturul din figura unde Y este un semicerc de rază p Avem atunci * ,,(i) L rez ^ + LC+'“(«ui S=SQ ] J— * unde piin(+ °° s-a notat integrala pe(—joo, j°o)\|—jp, jpj Obținem în continuare din ( ) Făcînd rezultă din ( ) y(t) = уЯ(°) + V-P- Punem reale și H(jo) = H(co) е^(“> și obținem din ( ) prin a celei imaginare ( ) explici tarea părții f + “ //(«) Sin ( /+ ф( ( ) ~ Jo w hi expresie care expliciLoază răspunsul indicial în funcție de modulul și argumentul funcției de transfer Dacă utilizăm explicitarea H(j fw)) L'(a>) sin ct>i , У( , obținem ГУ(со) sin ost = — Щ ) — v p do |- со Г ( ( ) Scăzînd membru cu membru ( ) dm ( ) obținem >( = v p E/W) sin л ■ -deo, Z > ~ Jo ( ) Expresiile ( ) și ( ) explicitează răspunsul indicial în funcție de partea reală și respectiv partea imaginară a funcției de transfer Datorită expresiei sale compacte, relația ( ) este remarcabilă, dezvoltările care urmează fiind în continuare consecințe ale prelucrării acestei relații Vom fi interesați în stabilirea unei forme tipice a graficului răspunsului indicial precum și a unor parametri care să caracterizeze sintetic această formă Pentru aceasta reamintim că (vezi § ) y( ) = o precum și (v ( )) = = H( ) = ( ) ultima egalitate rezultînd din ( ) ( ) ( ) Graficul pleacă în consecință din zero și tinde asimptotic către componenta staționară £ ( ) Evident, informații suplimentare asupra formei răspunsului indicial reclamă informații suplimentare asupra formei caracteristicii părții reale de frecvență Pentru început, vom considera că ( (co) este pozitivă și monoton descrescătoare (evident [ ( cu со —со), adică ( ) C (co) > , dco (v fig ) unde datorită faptului că t/(cb) este pară, Să încercăm să obținem informații asupra lui y(t), și ( ) Fixăm un t e [ , со) arbitrar Atunci, din sin / « = o dco rezultate din ( f sin u>t = rezultă că , , k = , , , ( ) și deci se poate scrie direct din ( ) j i/у = i%q — ^ # Г •” ( ) unde cu ( ( > ^ > со ( ) după cum rezultă direct și ( ) obținem din condițiile ( ) și ( ) Atunci j^)> , ) are o alură oarecare, atunci este posibilă totdeauna o descompunere de forma = i ~ Fig ( ) unde Lrf(°) Ș* (со) satisfac condițiile ( ) și ( ) (v fig , я și , ) Conform descompunerii ( ), răspunsul indicial va fi exprimai ca o sumă de forma unde Я, Я» 'W = £ W) - £ (*) i=l i=l * ( ) ( ) satisfac proprietățile ( ), ( ) ) Fig a b * Răspunsul indicial este У(^) = У+® - y-(t) ) prezintă «я singur maxim pozitiv și un singur minim negativ (y fig ) Din figura ( ) rezultă următoarea explicitare imediată a lui ( y(t) = ypp) - J (i!) ( ) unde dco I I C (co) sin ut ( ) Întrucît pe intervalul [ , co ] graficul lui U(y>) are alura din figura , a rezultă din ( ) cu ( ) ~ niax i i Să determinăm acum un majorant pentru | >’â(/) | Pentru aceasta, rescricm ( ) în forma C t (w) sin coi’ dw u / ( ( ) u nde Pentru funcția U(o>) descompunerea ( ) este І (w) — U (w) - ~ îl (w) după cum rezultă direct din figurile și ( ) ( Obținem atunci unde dw Utili zind o reprezentare de tipul ( ), cu considerarea lui ( ) și ( ) rezultă că " »o *- Utilizînd acum ( ), avem din ( ) majorarea o pentru în mod similar obținem și ininorarea №(/)>- , Umi„ Atunci din ( ) și ( ) rezultă ( Л О) I Ліѵ/ i mm» • v care introdus în ( ) furnizează majorarea max y(l) « , UmttJC -f- , Umi тін în vederea unor evaluări cantitative mai detaliate, caracteristice din figurile și se aproximează „trapezoidal" aproximație impusă chiar de forma acestor caracteristici Pentru început, să considerăm caracteristica trapezului din figura Se wdc din figura că C( ), (»= ( ( )-^— wd ~ — pulsația relativă de frîngere o precum și valoarea relativă ' ( ) Cu ( ), ( ) și ( ), ( ) se rescrie ( ) si care este răspunsul indicial normat corespunzător unui trapez unitar l ( = , £,= , cu pulsația de frîngere x e [ , ] Diferitele forme ale răspunsul indicial pentru cîteva valori ale parametrului x smt ilustrate m figura m Cunoscînd răspunsul indicial normat, parametrizat cu parametrul x se poate da o explicitate (aproximativă) a răspunsului indicialutilizmd o aproximare prin trapeze a caracteristicii părții reale a funcției de transter (v * Conform detaliilor din figura , și ținînd cont de ( ) ( > ) și ( ? \ obținem expresia răspunsului indicial de forma y(t) - ,( ) ^>(« , + IZ (O) y*^) + l a( ) ,y**( — indicele de oscilație T fiind variația relativă a amplitudinilor a două d - co^’ »-ționată de pulsația o a trapezului reprezentativ al acestei caracteristici Să vedem acum care este aspectul dezvoltărilor de mai sus în cazul sistemelor discrete, în ideea de altfel amintită a unei reductibilități la caz neted Fie așadar (?) o funcție de transfer strict stabilă, adică cu cz с ( x( ) și fie u(t) = !(/), /gK Atunci y(z) = H(z)z/(z — ) și deci: y(t) = — ( y(z) = A C AL s h ( ) Z IM — / Ы Introducem explici tarea //(jt)) - Я(^) ей{/,) f/(vj) l №) ( ) unde Й(т)) л I &(]■*)) Ь Atunci, cu ( ), obținem din ( ) H( ) i( со z —>— (deoarece |Я(— ) | )/(»-Л cefe- valentul neted este II (w) - - ) ’ unde ( )) ₽ - ( ) și pentru care corespund reprezentările din figura De remarcat că deși caracteristica părții reale a lui ft(w) este > și monoton descrescătoare, comportarea răspunsului indicial este total diferită de cazul neted (v figurile ), supraur-mărirea maximă fiind de % pentru caracteristica p— — din figura Remarcile făcute pe cazul particular de mai sus se confirmă și în general Intr-adevăr, considerînd ( ) și ( ) adevărate cu referire la С (т)) și observînd că sin (t • • arctg tj) = pentru t = tg —, t Fig k e ( ) se poate rescrie ca t v(^) = — ^ л=о ( ) ( ( ) Г ° t sin t arctg тк o C ( ) fK/ P) curbele de o și tt = ct Fig Fig ! i i • * £ ’ CONEXIUNEA SISTEMELOR REACȚIA NLGATlVĂ STABILITATEA CONEXIUNII ÎN REACȚIE NEGATIVĂ MV V К I f *• , (netede sau discrete) cu Fio două sisteme ( Ț, cî), v( intrătile nt și ieșirile u, i , ' Vom spune că cele două sisteme sînt conectate m — seric ducă «a = >'i ( ) definite ca o nouă intrare și respectiv o nouă ieșire; — paralel dacă | == '//- U ( ) iai ca nouă ieșire У У -I- У a — reacție (cu sistemul al doilea plasat pe calea de reacție) dacă ( ) ( ) ( ) ( ) cu w și у ca o nouă comandă și respectiv o nouă ieșire Semnul + în definește conexiunea în reacție pozitiva, iar semnul - în aceeași relație doinește reacția negativă (v fig а, b și c) Conexiune Jn reacție ț H ‘ w ?tan- * ace ipoteza eâ cele două sist,me fa ce n„ afeclea ă cu ]e • Conexiunea serie Din ( ) și ( ) ( ) ( ) f + ^ obținem un nou sis-Ьи, у = cc% pentru care Introducmd starea x А гxTxr tem numit sistemul rezultant, £ °L L& C A \ t^tdimSnsiune z(a) sînt fracții raționale ireductibile dm ( ' ) obțmem ' ? Я(Х) = -r^ r Л • г»/ л xu ( respectiv Atunci, dac?, rx(X) simplifică rădăcinile din a(L, , tenorul lui СГі(О) ale lui ^(X), polii lut Я(Х) se vor găsi m în ’ ( ) și sistemul rezultant va ti extern stabil netede Exemplul * Fio sistemele considerate i о(ЛА) •*» {— , ($ H ) ( — ) q(Av) t (s h ) (■- "• Atunci conform cu ( ) o( ) «« {— » » — , } sistemul rezultant nefiind intern asimptotic stabil, in schimb conform cu ( ) obținem //(>)= (s “ ) (s “h ) sistemul rezultant fiind insâ extern strict stabil și exemplul , evidențiază (aptul ca mod necesar minimalitatea sistemelor com- Discuțiile de mai sus, precum conexiunea în sene nu conserva țn ponente (dacă acestea sînt minimale) Conexiunea în serie poate fi generalizata prin inducție la un număr dc sisteme mai mare ca , relațiile ( ) și ( ), devenind n(zl) ■ ^МЛ) si respectiv H(X) = (~]Я lu % — d j v d“ ( ) у + cfza Rezultă similar ca în cazul precedent un sistem rezultant h nâ dimensional caracterizat de tripletul ( ) Din inspectarea tripletului ( ) rezultă că relația ( ) continuă să fie adevărată, iar echivalentul, în acest caz al relației ( ) «sv Я(Х) ( ) Hx(X) + Ha(X) în consecință, concluziile referitoare la stabilitatea кпіі^п^ № mimmalitale, stabilite în cazul conexiunii serie au aceeași valabilitate și in ca conexiunii paralel Prin inducție, relațiile ( ) și ( ) se genei*» i'/A 'i ' la cazul conexiunii în paralel a mai mult de două sisteme, în forma — а(Л ) și respectiv П(Х) <> • Conexiunea în reacție Utilizînd ( ), ( ), ( ) obținem și ( ) ( ) sistemul rezultant fiind descris prin tripletul Ц'гЧ J și avînd dimensiunea nr + и Să evaluăm spectrul lui A Avem: O == det O = det = det(XÎ — Лі) det (hi — A Ț "zPtvA ^i) °ісг/ = det(XÎ — Ді) det (ХІ — ri ) det(/Ț (X/ — Л ) b c^\I — Utilizînd formele primare de exprimare ale funcției de transfer H/X) -= v (X)/y (X) i = , , ( ) se sene m definitiv in forma ‘ xjx) = xjx) T vi(X) v (X) și care este expresia generală a polinomului caracteristic al sistemului rezultant în cazul conexiunii inverse лк+ п₽тп Hîn In ipoteza că ли»« v,(x)/zjx)xjă) mducMa, obținem : ( ) ( ) а(Л) = ( T formulă remarcabilă a conexiunii inee^^ * fa +; Inșpectind relația ( ), Ю modi[ica ftră ы іГа Tstelufoi UL k) «W rezultant? Kăspunsul la această întrebare este dat de: utilizat proprietatea cunoscută det(J + Л/N) det(/ + NM) Teorema Dtndu-se sistemul ( b c'{') de dimensiune n, presupus minimal (forma primară a funcției de transfer ѵх(Х)/%л(Х) este ireductibilă) si un polinom monic arbitrar x(X), avind gradul m, exista un sistem (A care conectat *n reacție negativă (pozitivă) cu primul (v fig, , , c) furnizează un sistem rezultant al cărui spectru coincide cu rădăcinile poUvomad/UA X(X), adică (v ) хл(х)s хл/х) хл/х) + vi(x)vs(x) ~ x'u Demonstrație Fie Z | (X) = » + a^X + + + Xя Ъ(Х) = Й ’ + ^’X + - P^iX- + Xя ( , ( ( ) și Ai eliminantul celor două polinoame presupuse coprime (vezi ) care este o ln X и matrice nesingulară (v teorema ) Punem X (x) = ao ) + aî ’x + - + a^iX»' + X® v (X) = p(“> + p^x + + e^iX”- + Xя ( ) ( , Fie de asemenea polinomul monic arbitrar, dar specificat, ( ) Cu ( ), ( ) și ( ) ecuația ( ) se scrie în forma / "xjx) " Х -лЯ [a a • ■ • Й j] ^■'Хл^) V (X) XV! (X) X'-iv^X) ( ) = Yo + YiX -p ••• "h (V» — “i) ') х'* + ”■ "I* (Yan-i ' a»-i) x" Utilizînd reprezentarea ( ), ( ) se rescrie în forma и> № № ••• Ш W) = = [Yo ■ Yt Y»-^ ’ ••• Y »-i~a!>-il y(x) ( ) Soluția (unică) a acestei ecuații este l«!r> a(f> ot!^! Й, ' = [Yo Yi ••• Yw— ’ ••• Yait-i—an- ] ( ) Odată aț \ & ), i = , , , n — determinați, sistemul căutat este orice realizare a fracției ѵ (А)/хл> (X), în particular, de exemplu, realizarea controlabilă Teorema este astfel complet demonstrată Consecința operantă a acestei teoreme se poate formula abreviativ și astfel: spectrul unui sistem minimal, poate fi arbitrar modificat printr-o conexiune în reacție negativă (pozitivă) cu un sistem convenabil ales în particular, dacă sistemul minimal este intern instabil, el poate fi făcut intern asimptotic stabil cu o locație a spectrului în C“ respectiv Uffp, corespunzător aleasă Această operație poartă numele de stabilizare (prin alocarea spectrului) Exemplul Fie sistemul neted, instabil Г avînd = { , } Avem direct v,(X) = Eliminantul celor două polinoame este M = o o fiind o matrice nesingulară Alegem de exemplu a(A) = {— , — , - ,- } adică yp-) — = (X - ) = i | X + X + л + >A Atunci ( ) devine: Rezultă [a P b , c?), astfel încît sistemul discret rezultant să aibă spectrul identic cu această locație; — se implementează schema din figura CNA UW neW U" (AtJ - rfd ГЛ Д,с/) u (tiist discret) Fig Exemplul Fie sistemul neted avînd Хл М — + x®> ѵг(Х) = Aplicînd formulele ( ) obținem (exercițiu pentru cititor) —sin h cos h Rei ’ltă XJjd(X) = —' kcos/r+ X , ѵ й (X) =( - cos h) ( -f- X) de unde se vede că, pentru h = ( Л + l)rc» k = , , , polinoamele ХЛ І(Х) Vui(X) — pcli-nomul ț- X ca divizor comun și în consecința operația de stabilizare nu se рэатз г elr rai nan tul polinoamelor ѵіл(Х) fiind o matrice singulară Fie așadar un ’ F\- — - e altfel arbitrar ales, de exemplu h = — - x vJd(X) — -p — X Alegem v(X) — X deci toate valorile proprii sînt locate în oiigine Aplicînd formula ( ) obținem I |X, ) aj ) PJa) PP>] = sistemul discret stabilizator fiind do exemplu Л Odată configurația de stabilizare din figura implementată (v spre ilustrare exemplul de mai sus) vom avea evident = Ад(Лй) — ► cu Ni э h — ► oo ( ) pentru orice inițializare a sistemului rezultant acesta avînd intrarea и (v ( ) și fig , c) identic nulă și unde ж, este starea sistemului (Av bt, cț), iar ;rw starea discretizatului acestuia Se pune însă întrebarea dacă și limita x^l) —> cu э t -> oo ( ) este adevărată, adică dacă operația de stabilizare are loc și la nivelul sistemului neted, care este sistemul original? Pentru aceasta fie un l> arbitrar ales Există atunci un astfel îneît t e [kh, [k )Л] Se poate scrie în consecință xx(/) — xfjth) + ( еЛі^“т) (т) eu h Luăm шип tu uuxl Intențional un Л w( /j p ) rv, /«eJS[și fUncțkv de transfer devine după simplificare Кщ(і)» —-— £ * l АрМчШЦ ptoedura din teorema » aceasta mit) acte uit ratei pentru orice realizare mini-ны?Л âv ex) a acsstei Jmcfii transfer pentru au ro avem dupl cum imediat rezultă vxW =■-• Xî W -a + l- «i “ Să punem» de exemplu : Хл(а) *= a* și vom obține a plicind formula ( > ) [aja) pj nj ® — l sistemul de pe reacție fiind deci Ла = , Ьа = ~» că = A«iest sistem va stabiliza numai realizarea minimală (Лх, Ьх,с'() și nu și sistemul original d й» c^)» în mod intenționat» să plasăm în reacția negativă a discretiza tului (ri^, (j, cfj sistemul ultim obținut (Лаі Cq ) Sistemul rezultant va avea conform cu ( ) (în care se ia semnul —), matricea Л (і Mi o i ’ ~ Matricea fundamentală va fi atunci Ф( = S-^jrZ - A)- = pentru - l pentru /»!, pentru (exercițiu pentru cititor) Dacă se consideră o inițializare do forma vom avea (pentru u(/)= ) д |(*) A«(/) = f = , obținem din ( ) explicitarea Я (Х) = -s AM ( ) A(X)У> (*) T n(X) r (x) PW Din ( ) rezultă imediat că dacă h(X)^(X)/^i(X)^) es^e ireductibilă* atunci și r(X)/j>(X) este ireductibilă în acest caz și numai în acest caz avem ^[Я (Х)] = зМЮ Ш T = ( ) = ЗВ) =F Ca și în cazul stabilității interne se poate da o formulare practic identica a teoremei pentru cazul stabilizării externe Deoarece, însă cele două tipuri de stabilitate nu sînt totdeauna echivalente, stabilitatea internă impli-cînd-o pe cea externă, problematica stabilizării, discutată mai sus se referă în exclusivitate la stabilitatea internă De aceea, operația de stabilizare descrisa, conform teoremei , în termenii funcțiilor de transfer, se efectuează cu formele primare ale acestora presupuse în mod obligator ireductibile și deci m fond operația de stabilizare se efectuează- la nivelul realizărilor minimale Altfel spus, operația do stabilizare este o problema sislemică internă și nu o problemă de tip i^rătej ieșire Considerînd acum conexiuni mat complexe, acestea se descriu in d descompuse, tot cu ajutorul celor trei tipuri de conexiuni dictai JJai precis ti polo ia unei conexiuni complexe se poate desene, (atenție ) n m •n termeni intrare/ieșire, printr-un graf de fluența n ra Л-Graful de fluență ascciat conexiunii din fig este redat in figura HM HM н,М н&) Transmitanțele cailor simple care leagă nodul и și nodul у sînt C = H H , C =H transmitanțele bu? ? o s - = —x? = ^ lui Mason este Нз iar Fig Formula и - Bt) = unde * elimină adiacența Avem atunci а >♦ i / * -h HaH - care, etic funcția de transfer rezultanta Să abordăm acum Analiza In fuMI« gativă) a sistemelor, istoricește, aceasta a * reacție fpe-constituit primul mod de abordate i ii li ' , cT) i = , conectate în reacție negativa (v, fig , , c, unde la sumator se ia semnul —) și fie H (s) = v (s)/x (s) formele primare ale celor două funcții de transfer Se face ipoteza că ‘ „fracția rațională t(s) Ms) eS|e irecjuctibilă“ Хл!(«) ( , ) Aceasta înseamnă, conform dezvoltărilor de mai sus, că atît fiecare sistem este minimal, cît și sistemul rezultat din serierea celor doua sisteme este, de asemenea, minimal în aceste condiții, din ( ) și din ( ) obținem egalitatea a(X) = ^[H (s)] = [ + H(s)] ( ) re unde (A, b, c ) este sistemul rezultant ( ) și unde în blocul ( , ) al ma-tricii / s-a luat semnul Problema analizei în pulsație a stabilității conexiunii în reacție negativă se formulează astfel: considerînd ipoteza ( ) adevărata, să se formuleze în termenii reprezentării H(jco) condiții necesare și suficiente de stabilitate asimptotică internă, adică de locație în a spectrului o(A) al sistemului rezultant, о(Л)с A arg ( + Я($)) |ser = ~Np ( ) (sensul trigonometric continuînd să fie sens direct) Să facem operantă dubla implicație ( ) Pentru aceasta, să observăm mai întîi că + я($)] = еГ- Fie No numărul polilor distincți de pe axa imaginară ai lui + Я(г), care sînt conform cu ( ) poli cu aceleași ordine de multiplicitate ai lui Я(х) Avem atunci A arg ( Я(з)) ser — A arg ( + Я(] Punem prescurtat As A arg ( + Я(]со))|иб( > OT)XJr ( ) и Se vede deci că A definit prin ( ) este variația argumentului vectorului cu origina în punctul critic și vîrful pe ramurile continue ale locului de transfer, atunci cînd parcurge în sens crescător intervalul ( , oo) (exceptînd punctele din M'\ Fie acum jco un pol de pe axa imaginară al lui H{s) avînd ordinul de multiplicitate r{ Avem atunci, dezvoltînd în serie Laurent în jurul lui jco : H{$) = - jof de unde Punînd No О Г- , г ‘ г adică suma ordinelor de multiplicitate a polilor lui H(s) plasați pe axa ima ginară, obținem din ( ) și ( ) că A arg ( “ Я($)) |seYi — Уо'* i=l ( ) in Ui ti în fine, cu Я(со) = , rezultă că A arg ( + H(s)) |sec = Д arS ~ = înlocuind în ( ) explicitările ( ), ( ), ( ) și A arg ( + H(s)) |ser = A — Cu ( ), dubla implicație ( ) se va scrie a (A) c A = (^o + ^p) у ( ) ( ) obțiiiem ( ) ( ) și care se citește sub forma enunțului din teorema care urmeaza, teorema cunoscută și sub numele de criteriul lui Nyquist Teorema (Criteriul lui Nyquist) Presupunem ipoteza ( )adevărată Atunci sistemul rezultat din conectarea in reacție negativa aS;mb- ( ») («■«Ptlnd punctele de discontinuitate), este egală cu (q + Nr) unde q și NP este numărul polilor de pe axa imaginară și respectiv dinG' a, lui H(s), Mmînd și ordinul de mulliplicitale vom considera * *•* «Wat шаі US “ "• '■“* * •— ь ь оИ МІЯ H(s) = intervalul ■Locul de lăsată ca Se vede transfer este ilustrat in figura t exercițiu cititorului, ’ rasarea lui fiind că unghiul a variază de la — (О, об) Deci: De asemenea, ?c = , AiP = ! și deci la - cind o, parcurg* oonstatîndu-se că este verificat criteriul : este invitat să probeze această afirmațieiеі^^гГ aSimptotic Cititori Exemplu! Să reIuăm £ ' lui Xyquist și deci Я( ) = CU hodograful ilustrat în figura , -Se vede că a variază de la n la și deci Pe care îl reluăm în figura > > Fig De asemenea, ? = și Л>= șj deci în bLcU ini eS , IVp = » (unde q este tipul sistemului, v § ) ■Cu ( ), ( ) devine и Deoarece !I(s) nu are decît poli în origine, locul de transfer va avea o singură ramură continuă pentru «e( , oo) Cu aceste precizări, tipologia formelor locului de transfer care respectă ( ) este ilustrată în figurile Vom nunii indice de traversare al locului de transfer, diferența dintre numărul de traT versări al intervalului (— oo, — ) al axei reale din cadranul II în cadranul III și numărul de traversări ale aceluiași interval, din cadranul III în cadranul II Fig Cu aceste elemente obținem: • Criteriul practic al lui Nyquist Presupunem îndeplinite ipotezele ( ) și ( ) Atunci sistemul rezultant în bucla închisă este intern asimptotic stabil dacă și numai dacă, indicele de traversare al hodografului este nul atunci cînd pentru co = + locul pleacă de pe axa reală sau din cadranele III sau IV, și este + , atunci cînd pentru со = + locul pleacă din cadranele I sau II (v și fig ) Este interesantă transpunerea criteriului practic al lui Nyquist în diagramele logaritmice, acestea construindu-se unificat în așa-numita diagramă Bode Diagrama Bode constă din suprapunerea celor două planuri de reprezentare a caracteristicilor logaritmice în sensul că, axa pulsațiilor este aceeași, dar ordonatele au originile decalate, la dB corespunzînd —к radiani (v fig ) Este evident că pulsațiile la care locul de transfer traversează axa reală în intervalul ( — oo, — ), corespund la cp(co) = — к și (co) > Atunci indicele de traversare în diagrama Bode corespunde la diferența dintre numărul de traversări al axei со, de către caracteristica fază-pulsaț-ie, de jos în sus și numărul de traversări de sus în jos, traversări contorizate numai în domeniul de pulsații j /(co)]dB > * Obținem astfel Criteriul lui Bode Presupunem ipotezele ( ) și ( ) îndeplinite Atunci, sistemul în buclă închisă este intern asimptotic stabil dacă indicele de traversare (calculat în diagrama Bode) este ~ pentru , iai pomoamee (s) și p(s) sînt hurwitziene O funcție de transfer cu astfel de proprietăți o vom numi de fază minimă în sens strict Obținem astfel: •Z * • • Criteriul practic al lui Bode Presupunem că funcția de transfer H(s) în buclă deschisă este de fază minimă în sens strict, iar caracteristica asimptotică amplitudine-pulsație posedă următoarele proprietăți: este monoton descrescătoare cu pulsația și are panta asimptotei de joasă pulsație egală cu sau — dB/decadă; este simetrică (sau cvasisimetrică) față de pulsația de tăiere, panta la această pulsație fiind de — dB/decadă Atunci, cu ipoteza ( ) îndeplinită, proprietățile și sînt condiții suficiente care asigură stabilitatea internă asimptotică a sistemului în buclă închisă într-adevăr, utilizînd ( ) cu A# = și q — , obținem iar indicele de traversai e este nul, ceea ce atestă stabilitatea internă asimptotică a sistemului în buclă închisă Precizare: Caracteristica asimptotică poate sa fie, cum dealtfel s-a menționat în paranteză în enunțul criteriului practic, cvasi-simetrică fața de pulsația de tăiere, cum este în cazul pantei nule a asimptotei de joasa pulsație, in acest caz, simetria exactă ar implica ca și panta asimptotei de înaltă pulsație să fie tot nulă, ceea ce este imposibil întrucît funcția de transfer este stuct proprie Cvasimetria revine în fond la faptul că într-un interval „suficient de mare" în jurul pulsației de tăiere, panta să fie de — dB/decadă (v tig ) Să transpunem rezultatele analizei de mai sus la sistemele discrete, fără, a dezvolta o teorie specifică, reducînd cazul sistemelor discrete la căzu; sistemelor netede ?|а, , Astfel, daca ipoteza ( ) este îndeplinită, stabiliiatea asimptotică internă a sistemului în buclă închisă este echivalentă cu t oii că a sistemului (discret) rezultant în buclă încinsă se poate proba, la nivelul condițiilor necesare și suficiente sau numai suficiente în cazul criteriului practic al lui Bode, aplicînd criteriile din cazul neted echivalentului neted de ordinul zero H(w) al funcției de transfer în buclă deschisă H(z) Exemplul Fie Echivalentul de ordinul zero este H(w) = Locul de transfer este ilustrat în figura , iar diagrama Bode în figura 'rac! к kdB Oh , - -w I A Din aplicarea criteriilor Nyquist si Bode, rezultă imediat stabilitatea asimptotică internă a sistemului în buclă închisă (exercițiu pentru cititor) Cititorul este invitat să constate efectul asupra locului de transfer și a caracteristicilor logaritmice, a faptului că H(w) este proprie SISTEME DIFERENȚIALE Șl CU TIMP MORT Am văzut în paragraful că dacă H(X) este o funcție rațională sfrid proprie, atunci familia (f a realizărilor acesteia, adică a tiipletelor (А, b, c ) pentru care c^’Chl — X)'"^ b = /Г(Х), este nevida în acest paragraf ne rom ocupa de extensiunea noțiunii de realizare lâ cazul cînd Я(Х) nu este strict proprie în acest caz, avem reprezentarea unica H(X) « rf(X) unde rf(X) este un polinom, numit cît sau parte mtieaga, й(Х) = d$ -p dfK -p “ rf,X (^г=/= ), ( ) г Я(Х) ( ) iar f?(X) este o rațională strict proprie Atunci uin ( ) rezulta imediată că avem Pe de altă parte,-noțiunea de funcție de transfer este asociată comportării intrare/ieșire, adică regimului forțat (v ( ), ( )) și în consecința veni scrie j(X) = Я(Х) ЦХ) == Я(Х) «(X) + d(f) u(f) = л(Х) + УгМ ( - ) unde h(X) А Я(Х) й'Х) ( ) și J (X) Aiî(X) «(X) ' ( ) Din ( ) rezultă că numai părții strict proprii Я(Х) a lui Я(Х) i se poate asocia noțiunea de realizare în sensul definiției și deci numai tranziția intrare/ieșire ( ) este susceptibilă de o interpretare sistemică 'ncadrabilă în axiomatica expusă în capitolul Așadar, ( ) rămîne să aibă exclusiv numai o interpretare intrare/ieșire căutînd să punem totuși în evidența dacă este posibil, o interpretare convoluti^ă • Vom începe cu cazul neted, în sensul că vom da o interpretare expri mării prin transformarea Laplace a relației ( ) pentru X == s, adică yfs) — d(s) u(s) ( ’ ) Fie atunci, în algebra cu convoluție a distribuțiilor cu suport pozitiv „ (r) A Wo§ + dt + + dr ) * и, V « e ф'+ ( ) sau explicitînd, (r) j/g = dQu , rezultă ca o tranziție pur anticipativă și deci necauzală Spre exemplificare, dacă d{z) = z, obținem y(i) = t ( ) unde t se numește timp mort sau timp de propagare (v figura ) Dacă «( ) este o funcție original atunci, aplicînd transformarea Laplace obținem r j/(s) == e~*Xs) ( ) și deci funcția de transfer (în sensul raportului jyțs)/u(s')), asociată tranziției de tipul ( ),'este: = e-f* ( ) • O primă problemă care se pune în ► jegâtur£ сц procesul de prOpagar e ( ) Г’? este problema realizabilită Pe scurt spus, dezvoltarea în jurul originii a raționalei ( ) pînă la puterea k a lui s, coincide cu dezvoltarea în jurul originii pînă la aceeași putere a lui e“ts Un astfel de tip de aproximare se numește aproximație de tip Pade a lui e~rs (de ordinul n + ^)- Spre exemplificare, să considerăm aproximațiile Pade pentru n = Relațiile maximale ( ) devin: —T + ri — T — + ă = Atunci, — pentru Л = : я, = « = obținem din pri: T^ d = t, b = — și deci ( ) pentru k = : « = , obținem din primele trei egalități — — t adică ( ) — pentru k O = , adică obținem din s ( ) T\/ T S + S > Așadar, aproximațiile Pade de ordinul ( + ), ( + ) și ( ), ( ) și ( ) Orice realizare a raționalelor ( ) ( ) și ( ), o vom numi realizarea finit dimensională de tip Pade de ordinul \n + k) a lui e- * Cititorul este invitat să probeze acum că e lizare (exactă) finit dimensională © O a doua problemă care se pune în legătură cu procesul de propagare ( ), este influența prezenței acestuia, într-o conexiune în reacție negațivă> asupra stabilității sistemului rezultant Pentru exemplificare, să considerăm conexiunea în reacție negativă din figura , unde, pe calea directă, funcția de transfer este asociată cu o realizare minimală a sa (de exemplu x =Ku, у = %) Funcția de transfer în buclă deschisă este atunci nu are o rea- H{s) — — e jco COS СОТ Д I «I HI COT KlS ( ) locul de transfer rezultînd din • • K (cos COT — j sin сот) sin COT sin X COS Л" ■ ți fiind ilustrat în figura pentru Д* l, împreună cu locul de transfer corespunzător absenței timpului mort (ț * ), ceea ce corespunde la k* ~ Să observăm că timpul mort, din punct de vedere al alurei, spitalizează în jurul originii» locul de transfer Mai mult, dacă pentru t и , sistemul în buclă închisă este, conform criteriului lui Nyquist, intern asimptotic stabil, oricare ar fi valoarea pozitivă a factorului de amplificare К (v locul = ), în cazul prezenței timpului mort (t > ), după cum Fig rezultă, conform aceluiași criteriu, din inspectarea locului de transfer k* = , sistemul rezultant își condiționează stabilitatea asimptotică de valoarea lui k* Mai-exact, după cum rezultă din figura (locul punctat), sistemul rezul tant este intern asimptotic stabil pentru • ■ • ■ *- ' ' Л * i / y{z) = L ф) ( ) adică ’ rezultînd că funcția de transfer, pentru т > , este o funcție raționala strici proprie O realizare imediată a lui ( ) este de exemplu cu Л gR* Așadar, în cazul discret, tranzițiile de tip propagare au realizări sistemice finit dimensionale, dimensiunea fiind egală cu numărul de pași asortați timpului (discret) de propagare ( ) SISTEME DE REGLARE AUTOMATĂ Capitolul de față, sub raport teoretic, este beneficiarul principal al noțiunilor introduse și al rezultatelor obținute în capitolul anterior; din punct • de vedere istoric, însă, problematica analizei și sintezei sistemelor de reglare automată este aceea care a generat și cristalizat conceptul de sistem dinamic, în particular de sistem liniar, conducînd la închegarea de sine stătătoare și totodată cu caracter generalizator, a construcțiilor teoretice unitare din capitolul Prin însăși obiectivele intrinseci și structura lor, sistemele de reglare automată au impus o semnificație nouă, teoretică și aplicativă, conceptului de conexiune în reacție sau inversă (feedback) Astfel, dacă în cadrul abordărilor circuitale, conexiunea în reacție (negativă) avea o destinație eminamente stabilizatoare (deci interesantă sub aspectul perturbațiilor trecătoare) sistemele de reglare automată au dat o utilizare calitativ nouă conexiunii inverseA utilizare care vizează comportamentul sistemic față de mărimile exogene cu caracter persistent Dacă sub aspect pur stabilizator, conexiunea inversă sau cum se mai spune „bucla de reacție" este o entitate globală și „simetrică", luată ca atare, din punctul de vedere al sistemelor de reglare automată, „bucla" prezintă o serie de puncte preferențiale impuse de localizarea mărimilor exogene și de evidențierea expresă a așa-numitei mărimi de eroare,, mărime avînd o destinație remarcabilă și anume una decizională Așadar, fără o pretenție definitorie riguroasă, vom spune că, sistemele de reglare automată sînt sisteme cu conexiune inversă, care își decid comportamentul față de mărimile externe pe baza mărimii de eroare automat generate, cu scopiri expres al anulării acesteia', elaborarea deciziei presupune o anumită cunoaștere apriorică a modelului mărimilor externe, ceea ce se reflectă prin prezența în sistem a unei copii a acestui model, copie numită și model intern Așadar, pe scurt, reglare automată înseamnă un proces sistemic decizional (de comandă), desfășurat în mod automat pe baza mărimii de eroare pentru anularea acesteia în contextul considerațiilor de mai sus, noțiunea de conexiune inversă nu mai prezintă nici un fel de filiație circuitală, prezența sistemică a acesteia fiind actualmente ridicată la rang de principiu, constituind așa-numitul principiu al modelului intern și conexiunii inverse Se ridică imediat întrebarea: care este raportul între noțiunile de reglare automată și respectiv de conducere automată, ultima noțiune avînd în prezent o semnificație unanim acceptată ca exhaustivă în general, așa cum am menționat mai sus, sistemele de reglare automată își explicitează funcționalitatea prin anularea mărimii de eroare; sistemele de conducere incumbă însă funcțiuni multiple și mai complexe, incluzînd și funcția de reglare (adică de anulare a erorii) Această înglobare a funcției de reglare se face de regulă în cadrul unei anumite „ierarhizări" funcționale, ierarhizare care definește prin nivelele sale cadrul și obiectivele- operației de conduceie într-o atare ierarhizare, funcția de reglare este de obicei plasată la primul nivel, nivelele superioare, fiind alocate funcțiunilor de identificare, optimizare, prelucrării datelor (în scopul predicției) etc , discutarea acestor funcțiuni ieșind din cadrul prezentului capitol De precizat însă că în orice structură de conducere, prezența funcției de reglare este obligatorie, fiind o necesitate sistemică intrinsecă Revenind la sistemele de reglare automată, ele se clasifică convențional după obiectivul final al funcției de reglare în două mari categorii: — sisteme de urmărire (sau servo sisteme), în care funcția de reglare (adică de anulare a erorii) are ca efect final urmărirea cît mai fidelă de către mărimea de calitate a sistemului (reglat) a unei mărimi externe, numită mărime de referința', — sisteme de rejecție a perturbațiilor (sau cu referință (consemn) fixă), în care funcția de reglare (adică de anulare a erorii) are ca efect final menținerea constantă la o valoare prescrisă, independent de „contextul perturbator**, a mărimii de calitate a sistemului (reglat) Să menționăm că îndeplinirea funcției de reglare nu se face exclusiv numai pe baza mărimii de eroare, ci, atunci cînd este posibil, în special în rejecția perturbațiilor, decizia (comanda) în sistem se poate lua și pe baza măsurării directe a mărimilor perturbatoare, cu condiția ca aceste mărimi să poată fi efectiv măsurate; este cazul reglării cu acțiune directă (feed-forward) ,v în fine, din punctul, de vedere al terminologiei nu vom utiliza termenul de sistem automat, dealtfel, din păcate încetățenit, mtrucit acesta nu expli-citează clar la ce se refera adjectivul de automat, cremd adesea confuzii sau generînd noțiuni ce au cu totul altă semnificație (v de exemplu noțiunea de automat (finit sau infinit)) FORMULAREA PROBLEMEI REGLĂRII CONSIDERAȚII SI REZULTATE CU CARACTER GENERAL J Fie un sistem liniar cu o intrare și o ieșire scris în forma sa completa z — dTx gi este mărimea de comandă, v e IR este mă- — • ♦ К ♦ " este starea, и e unde este starea, и e лиг l w , rimea perturbatoare, jye® este mărimea masurata (sau de leș re), e Sg este mărimea de calitate л J t Vom considera de la început că mărimea de calitate este direct măsurabila, adică: și vom numi mărimea у mărime de ieșire sau mărime reglată, b uncția de ti ansfer de la comanda и la ieșirea у este : y(s) tr/ x = ( ) № Ж ( ) (unde /%($)•• (unde r(s) Av(s), ^s) A z ,(«)) - Yț )s + - + y£-i^-‘ + sr* ( -l ) unde flr(s) este un polinom arbitrar de grad strict mai mic decît ;ir(s) Deci: ( = Ij'r'J'rW = б- —— > Var(s) e cu ^Ia(s)] QO £->oo ceea ce este echivalent cu ipoteza Polinoamele p (s) și ^(s) au toate rădăcinile în * • f I г ь г aceasta fiind о rațională proprie ( [rc(s)] л A(s) w(s) = rc(s) e(s) (З И) *^• J w % * * • • • Sistemul rezultant ( ), ( , ), ( ) și ( ) este: i = Ax у = cTx ( ) l sau eliminînd u, e și jy obținem : X = (A - hgccr) х + ЬДхс + bgcyr + cv Xc= -ЬссЧ + '-Ъсуг (З І £ = — СТЯ - Ут Introducînd стА[-^ ] ecuațiile ( ') se scriu cu ( ) în forma compactă: ijj — АдХд + bryT -j- b^v E = ( ) și reprezintă ccvcițiUc sistemului Tczuliuni in buclă închisa, cum se mai spune Cele de mai sus se explicitează în configurația din figurile , a și b Cu ajutorul reprezentării ( ) condițiile (S) și ( € ($) y(s) = r(s) u(s) -j- >(SMS) e(s) = ЗФ) - y(s) pc(s)w(s) — rc(s)E(s) care prin eliminarea lui ф) se scriu în forma matriceală r^>(s) —r(s)"l Fj'țs) p a>(S) Hs) /Ф) JU(s)J L° >'(«)' Se constată (cititorul fiind invitat să probeze) ca , x д л t , (v- “s- 'J’ “ Și Ѣ» figura a (utilizind eventual u„ grai de ftaențl asociat (v fig » obținem cu ( ), ( ) și ( ) л/ЛЛ/л /A fvds) - K«țs)t>(s)] «= —- înlocuind în ( ) expresiile lui S) ?’ obținem: , Kv л / \ / Al Л| ^(S) H,(s) a[w(s)D- ( - ) Fig Atunci, condițiile (S) și (A) capătă în limbaj operațional (sau polinomial cum se mai obișnuiește să se spună) următoarea explicitare echivalentă: — (S ) Polinomul vr(s) să fie hurwitzian, după cum rezultă din ( ); — utilizînd teorema valorii finale cu referire la ( ), în vederea explicității condiției ( ?), obținem de asemenea: !(T?'l) Я$ф)]сС“ ( ) Var(s), ^( )е]ВД, (а[бф)] * dr г», е [«,( )]} ( ) Рѵ(г) exprimare ce evidențiază în fond versiunea discretă a ipotezelor X- și o n ceea ce privește ipoteza , de persistență a mărimilor exogene (discre t sau discretizate), aceasta este automat îndeplinită întrucit prin operația de diștre-tizare, așa cum de altfel este cunoscut din capitolul anteuoi valori proprii din CT (CJ ) se transferă prin discietizaic ui Іц ) (С/ i( )/* Trebuie făcută însă precizarea de principiu ca dacă foima de undă a uiaii mii u, care comandă efectiv sistemul neted, este o foima de unda m scara, așa cum este ea furnizată efectiv de convertorul numenc ana ogic, plma; undă a perturbației v, avînd același tip de alura, este numai ipo e ica și încadrează în „erorile" sistematice tăcute prin operația de ciscre izare (v fig , a și b) Operînd acum în figurile , a și b, modificarea formală a literei s cu litera z, obținem varianta discretă a structurilor din cazul neted, structuri provenite în fond din configurațiile detaliate în figurile , a și b Cu precizările de mai sus se poate foi mula, pentru, cazul sistemelor discrete • Problema reglării Considerînd ipotezele II, , și , adevărate în versiunea lor discretă, să se determine un sistem (discret) compensator: xc(t + ) = Acxc(t) + bcz(t) u(t) + &e( ( ) astfel încît, sistemul rezultant ( ), ( ) și ( ) adică explicit : + == + Ш) + W) e(/) = c^xR{t} ( ) cu xR> AR) br, bv și cR furnizate de ( ) să posede simultan următoarele două proprietăți: (S) а(Лл) oo oricare ar fi și yțeCUr explicitate prin ( ) și respectiv ( ) Tot în varianta discretă avem și ' ’ J ’ • Problema r eglării (rcfor mulată in limbaj polinomial) Considerlnd ipotezele II, , I și adevărate în versiunea lor discretă, să se determine, pentru sistemul d iscret în reprezentarea polinomială de stare parțială: г (г)ф) ( ) p(z) У(г) “ z(z) u(z) un sistem (discret), compensator />с(г) u(z) = rc(z) s(z), e(z) ~л(г) ~ y(z) ( ) cu (X) jy(X) = r(X)î/(X) + ^(X)fl(X) e(X) = ^г(Х) — y(X) ( ) j >c(X)^(X) = ^e(X)e(X), unde )> = s sau X — z, după cum problema reglării este formulată pentru cazul neted sau respectiv cel discret în consecință, condiția (S) sau (Ș) din problema de reglare, condiție care se referă la stabilitatea asimptotică internă, trebuie să reflecte fidel în limbajul polinomial, o condiție de locație a spectrului matricei Лл(ст(Лл)) c (х)Л(х) + r(X)re(X) Хл(Х) «= det У>(Х) —r() jcW А(х) Pe de altă parte, cititorul poate constata direct sau referindu-se la paragraful , că dacă „fracția rațională r(X) rc(X)//>(X) Д (Х) este inductibilă" ( ) atunci зЬИ)] [i и- н(х)ад ’ k / ) ( ) unde, după cum imediat rezultă: - Н(х)Яс(х) - ■- ^ , '^ ,,^i ■ r(X)rc(X) exprimare ce atestă clar importanța condiției ( ) Din inspectarea lui ( ) se vede că în general; ț Л ) /* i \ Tf , (,/■ P ' [Хл( Й^З[ + (X) fc(X)], Г I * \ • / > / А incluziune care devine strictă dacă ( ) este infirmată, adică polinoamele r(X)rc(X) și ^(Ă) ^>C(X) nu sînt coprime Totuși, dacă „polinoamele r(X)/c(X), jâ(X)/>c(X) au c m m d c cu rădăcinile plasate în \ (S“) [ + H(X) H (X)] c (x)/>С(Х)/[АГ(Х), tcc(X)«“A(X) re(Ă)/pB(X) Cum condiția (S) este îndeplinită, adică з[хя(Х)]с cu X -> со dacă ЗІХіА)] со dacă зЕХл(Х)] c ^i( )- t - ) 'v A \ • - / ■ , , Așadar și condiția (R' ) este îndeplinită ’ Condițiile ( ), ( ) sînt necesare Intr-adevăr, din arbitraritatea pou-noamelor «r(X) și aB(X) (v ( ) și )) avem din ( ) ' ^(Х)Ж Лг(х) = i, «г,(х) = о xA) нА) s(X) =J sau ( ) —, ar(X) = , uB(X) = к xA) нА) Cum зЕхАМнА)] = зЕхА)! rізЕ нА)] = (V ( ) rezultă că pentru satisfacerea condiției (R ) iesPe? tl trebuie să dispară complet prin simplificare în consecință, condiții e (o ( ) se impun în mod obligatoriu drept condiții necesare Să punem în evidență prin cîteva exemple efectele teoremei Exemplul Fie sistemul de reglare din figura pentru care consideram (*r(s) = s și (lUs) = - Aceasta înseamnă că /r = ІУгІ>'г(/) -J «p® cori l (h eln lt V - *u h l ndhă luărlndb do rcfnrlnțft atnt do tip Iar perturbațille tfînt dc tip armonia (de pnlsațfe Sn vede «a sltdonud OHto іівногін prin /»(«) •* (« + ) (* f*w(s) » Ăa "I- l rompcn^dovul oato doaorto prin ’ Ker ultă atunci că Xb(«) "!>(*) /’«(«) + r o asemenea, |rr(s) st |Xy( ) и ? + | pc(^v(s) = (sQ - )/>C(S) fiind satisfăcute în consecință» conform teoremei și condiția (H) condițiile ( )» ( ) este îndeplinită pr(z) => p v(z) •=> ч - a lică mărimile exogene s!nt de tip treaptă unitară l(^) te У'(^) « Д , , t(r) , r (v) />„(r) = г — , f’ot't) “ ~~ ’ „ Avem atunci din figura obținem: polinom care aro tonte глИЛсіііПе рШЙ^Ѳ îu condiția ( ) c stt îndeplinită, 'menea (XrW ■» HvW “•■г ~ ^(*) ” s ~ condițiile ( ), ( ) fiind îndeplinite în consecința și condiția de reglaie (Ii) este satisf Să trecem acum la rezolvarea problemei reglării* Fie, în acest scop introduse următoarele polinoame, succesiv definite: ^(X) Ac m nulc {p r(X), />(^)}- ( - ) ■ ( ' ) ^(X) c m m d c {р^(Х)л ( ) Ш A ( * к • Cu aceste precizări, se vede imediat că problema reglării este echivalentă determinarea polinoame,lor zc(X) și jt>c(X), astfel incit să avem simultan Și unde (v ( )) a^(x)] iar M este o WiX »i matrice Procedăm în continuare astfel: I Se alege un polinom monic Xr(X) de grad «t — ale cărui satisfac incluziunea ( ) Fie acest polinom Xs(*) =₽o +₽iX + + Р В,- Х и>- + X "*- Se construiește eliminantul M al polinoamelor />х(Х) și r(X) Se calculează (prin rezolvarea adecvată a ecuației) l‘^l) ( ••• Ля ; /’о Ьі ] = (Pij pi р и,— ] zerouri și punem АД) = ao + «iX - + aMi t X"*- ( ) rcM =fy) + i’iX - + Z^-iX’i- ( ) de unde cu ( ) rezultă />C(X) Cititorul este invitat să arate că ( ), ( ), reprezintă soluția ecuației polinomiale Ă(X) A(X) - re(X)r(X) = Xr(X) ( ) (v paragraful ) Din ( ) avem ? (X)A(X) + , (X)r(X) t хл(х) xJĂ) Dar [r(X)] (X)] și deci cu atît mai mult d[/(X)] și deci £ (X) X xAOO/Xr(M t adică obligator d[^c(X)] = — Așadar ak(x)] ($) ч= r(s) «a r^s) l tc(s) — sa -h t ’ И • />x(s) »A p(s) V(s) Я= tt d" S> И » І aplicăm algoritmul din teorema , Se alege un polinom хл(«) do grad n — l » cu rădăcinile în ( h de exemplu d~ Ь sa d- ® -ț- d' s I I U у ^ ! £w;' ii Sc construiește ellmir mti '* cos i' ~ * de unde Д \lege m AW A’faW I ■* FixîUn de exemplu h я/З și obținem Г(Я) « Я /> W И г EKndnantul polinoamelov px(s) și г(я) esto o o rc/ o о о о о о о о о о о de ultima linie din M~l cere С > Ьс Datorit'! structurii polinomultii %K(^), ( ) este furnizată б/я /я KezultvL iar Obi пеш astfel schema din figura Să reluăm acum problema reglării dinfr-un punct de vedere mai larg și anume cel al stabilității structurale a soluțiilor acesteia Pentru a da un cadre riguros noțiunii menționate vom introduce cîteva elemente pregătitoare Fie în acest sens sistemul S = (A, b, cT, e) considerat în versiunea sa netedă ( ) sau în cea discretă ( ) Specificarea acestui sistem, ca dată inițială pentru problema reglării, se face prin elementele matricelor componente Л, b> c și r Aceste elemente pot fi asimilate cu un punct avînd N A» + n coordonate, , astfel îneît pentru orice ф' e ]HLX cu j| ф' — ф | ' Din definiție se vede imediat că o soluție structural stabila a problemei reglării în punctul ф este o soluție obișnuită în acel punct Fie așadar Se ° soluție a problemei reglării în punctul ф Scriem, atît sistemul v asociat punctului ф, cît și compensatorul Sc, în reprezentarea polinomială de stare parțială și obținem ecuațiile ( ) Faptul că Se este o soluție a problemei reglării în punctul ф revine la satisfacerea simultană a condițiilor ( ), cu ХЯ(Х), explicitat prin ( ) și, a condițiilor ( ), ( ) din teorema , condiții pe care resenem împreuna: ( ) cu ХЛ(Х) ^PW pe№ -I- r(X) rc(X) Ur(x) I p(\) A(X) H (*) I PcW r„(X) ( ) ( - ) ( Se constată direct că modificarea arbitrară a punctului p într-o vecinătate oricît de mică a acestuia, atrage după sine modificarea tuturor coeficienților polinoamelor p^K), /(X) și /Ѵ(Х) (cititorul este invitat, utilizînd teoria realizărilor făcută în paragraful să demonstreze riguros această afirmație), în consecință, pentru a menține în continuare adevărate condițiile de divi-zibilitate ( ) și ( ) este obligatoriu ca M*) |&(X) ( ) £ ИД) IPM ( ) Punînd p(X) Ac m in m c {p r(X), |х,(Х)} ( ) și denumind p(X) polinomul asociat mărimilor exogene, condițiile ( ) și ( ) sînt ambele echivalente cu condiția de divizibilitate P(X) IA(X) ' ( ) Pe de altă parte, datorită continuității rădăcinilor unui polinom în raport C ' rezulta că există un e > , astfel îneît pentru orice ф' e ' cu j ф —ф \ (X) și r(X), incluziunea ( ) continuă să fie adevarata Utilizînd atunci teorema , obținem, cu discuția de mai sus, • Teorema Dacă £c este o soluție a problemei reglării în punctul fp, atunci acesta este o soluție structural, stabilă în acest punct, dacă și nu nai dacă Teorema Dacă £c este o soluție a problemei reglării în punctul ($), f'(s)), specificată mai sus Avem Pi(s) — Eliminantul polinoamelor p^s) și r(s) este о о о о о о о м = = м- о о о Alegem $ obținem Mică, pc(s) == lx(s) P M ~ bținem astfel schema de reglare din figura , Уг + s fie Din dezvoltările făcute imediat mai rus, cît și din acest exemplu se des-prind cîteva observații remarcabile- v » , v L Sinteza unei soluții structural stabile este independenta de specificarea matricei e și respectiv a coeficienților polinomului rfs) în consecința specih-caisa acestor elemente nu este necesară în construcția unei soluții structura, stabile a problemei reglării* Se mai obișnuiește să se spună ca o soluție structural stabilă a problemei reglării este robustă în raport cu matricea e sau cu coeficienții polinomului r€(s) Compensatorul £c, scris în reprezentarea polinomiala de stare parțiala, л ($} k(s) — r (s) e(s), se vede că înglobează în mod obligator o copie a modelului mărimilor exogene, în sensul că p(s) este un divizor al lui pfs), adică pfs) — ^)?c(s) Atunci sistemul Sc mai poate fi scris și în forma; p(s) £($) = e(s) ( - ) pfs) u(s) = £(s) ( - ) „Subsistemul" ( ) poartă numele de model intern, prezența acestuia fiind o condiție necesară pentru existența oricărei soluții structural stabile a problemei reglării De observat că k(ă) | p(X) (v ( - )), și deci orice compensator care este o soluție a problemei reglării conține „o parte" din modelul mărimilor exogene, „restul" modelului fiind conținut de către proces, adică înglobat sub forma unor divizori ai polinoamelor p(s) și rfs) (v rel ( ) — ( )) Mai exact: modelul referințelor, sau o parte din acesta, se localizează ca un divizor al lui p(s), pe cînd modelul perturbațiilor, sau o parte din acesta se localizează ca un divizor al lui rD(s) Acest mod de localizare este obligatoriu și orice permutare a acestuia conduce la rezultate profund eronate Discuția făcută impune în fond distincția dintre o soluție obișnuită a problemei reglării și una structural stabilă, în sensul că: un compensator pfs) u(s) = rc(s) e(s) care este o soluție obișnuită a problemei reglării înglobează acea parte din modelul mărimilor exogene care nu este conținută de către proces, adică z s) \pfsj, pe cînd un compensator care este o soluție structural stabilă a ргсЫстеі reglării înglobează, sub forma modelului intern, adică p (s) I în totalitate modelul mărimilor exogene, chiar dacă acesta mai este conținut în parte sau complet de către proces Spre exemplificare, cititorul poate constata că în exemplul (în cazul neted și în cazul discret), avem h($) =/(? + compensatorul înglobînd numai o parte din modelul p (s) al mărimilor exogene și anume k(s) — s - în exemplul însă, unde x(s) == s, se vede că, cu toate că procesul înglobează modelul intern, acesta se mai regăsește încă o dată și în compensator, fapt ce conferă soluției calitatea de stabilitate structurală Faptul că un compensator, care este o soluție structural stabilă a problemei reglări), înglobează modelul intern, copie a modelului mărimilor exogene, iar acest model este activat de eroare (v ( )), eroare produsă cu ajutorul conexiunii inverse (v ( )), exprimă ceea ce actualmente este cunoscut sub denumirea de principiul modelului intern și al conexiunii inverse ANALIZA SISTEMELOR DE REGLARE AUTOMATĂ h «l«tE^r^rafulixantCiior au, fost eXPuse principalele rezultate privind Ti ' саоЫя ѵЬ«??і’ sistemelor de reglare automată, problematică ce a avut raior'^A w ’ț’s nu!':i unui compensator care este simultan stabili' Mor și regulator, adica conferă sistemului rezultant (în buclă închisă) pro- prietățile ( ) și (A‘) Soluționarea problemei s-a făcut într-un cadru general din punctul de vedere al arbitrarității mărimilor exogene, adică arbitrari-tății ,’olinoamelor p , (X) și |л„(Х) (sau |i(X)) Cu toate acestea, în aplicațiile concrete se impun sistemului de reglare automată proprietăți suplimentare celor fundamentale, adică proprietăților (S) și (R) Aceste proprietăți suplimentare explici tează ceea ce se numește calitatea procesului de reglare automată, calitate ce este descrisă convențional printr-o clasă de indici sintetici, care caracterizează, cu o exprimare unanim acceptată, performanțele sistemelo' de reglare automată Definirea indicilor sintetici de calitate și deci implicit a performanțelor, precum și aprecierea (într-o manieră inginerească) a eces/ora, fac obiectul analizei sistemelor de reglare automată Metodologiile de analiză „de tip convențional", care de altfel fac obiectul prezentului capitol, au în vedere în principal, comportamentul intrare/ieș’ e al sistemelor de reglare automată și, în consecință, vor opera în cvasitotali-tatea lor cu funcțiile de transfer și cu modurile de reprezentare al acestora Așadar, structura de bază pe care se vor dezvolta raționamentele din acest capitol este cea din figurile , a și , a, structură pe care o reluăm, pentru fixarea ideilor, în figura unde, așa cum am mai procedat, X — s corespunde sistemelor netede, iar X = z, celor discrete Graful asociat structurii din figura este ilustrat în figura (v și figura ) Se introduc următoarele funcții de transfer reprezentative asociate structurii din figura Funcția de transfer în circuit deschis, ca fiind funcția de transfer de la eroarea z la ieșirea у Д,(Х) ~ ~~ K- f^(X) WC(X) e(X) ( ) și care este egală cu transmitanța cu semn figura , B(X) = — (X) / „(X) Funcția de transfer în circuit închis, ca fiind funcția de transfer de la referința yr la ieșirea у л v(M„ к-гт**»****», fr-urf У,№ ( ) •i /iW rezultă imediat din absorbția buclei grafului din figura , schimbat a buclei* grafului din * Ceea ce juslificft Indicele „Ь" « fiinoțifl ч • к " perturbațiile fiind de același tip ȘP adică dentru sistemele netede și cu menținerea egalității ( ), pentru sistemele discrete Așadar yr(t) — (£) - p = * , , ( ) cu v(t) de același tip, pentru sistemele netede și, JO t proprietate detaliată prin — performanțele regimului staționar (al erorii), descrise de indicii sintetici corespunzători РеС^ а*,е a analizei sistemelor de reglare automată este aceea că aprecierea indicilor de calitate care descriu sintetic performanțele acestora, adică comportamentul în circuit închis, se face din datele și specificațiile sistemului în circuit deschis, adică ale funcției de transfer M^s) șt ale diferitelor tipuri de reprezentări în pulsație ale acestora Din acest ultim punct de vedere, analiza sistemelor de reglare automată este în fond o analiză prin pulsație (sau frecvență) ANALIZA STABILITĂȚII (ASIMPTOTICE INTERNE) Analiza urmărește în acest caz, verificarea proprietății (S), adică a inclu iunii trț ljj) c<> în cazul sistemelor netede și, а(Л,) cZ i( ) în cazul sistemelor discrete, incluziune explicitată polinomial prin locația zerourilor polinomului Xr(X) = />(X) />C(X) + r(X) ,rc(X) (v ( )) Aprecierea stabilității se poate face direct, prin aplicarea, de exemplu a criteriului lui Hurwitz polinomului /Л(Х), în cazul sistemelor netede, sau polinomului /л —!—— ) în cazul sistemelor discrete \ - Н/ Se preferă însă, ca fiind o procedură inginerească comodă și eficientă, să se utilizeze rezultatele analizei prin pulsație (frecvență) din paragraful / , referitoare la aprecierea stabilității asimtotice interne a sistemului rezultat din conexiunea în reacție negativă a două sisteme Mai exact, considerînd valabile ipotezele ( ) (mai dură) sau ( ) (mai relaxată), stabilitatea asimptotici internă, se apreciază, în cazul sistemelor de reglare netede (discrete), л \ direct din funcția de transfer în circuit deschis Hfs), I Hb(w) = Hb(z) I z = - I ț — wf prin aplicarea criteriilor Nyquist sau Bode, aplicare detaliată în dezvoltările făcuie în paragraful și asupra cărora nu vom mai insista Vom mai insista însă asupra a două considerente, care în diverse modalități de formulare au mai fost expuse Cu toate că criteriile analizei prin pulsație referitoare ia aprecierea stabilității (Nyquist, Bode) și operînd cu funcția de transfer în circuit deschis probează în fond incluziunea ( ), adică condiția (S")t care este o condiție necesară și suficientă de stabilitate strictă externă, aceste criterii sînt utilizate de fapt pentru a furniza informații asupra stabilității asimptotice interne (y ipotezele ( ), ( )) Maniera de utilizare a analizei prin pulsație în investigarea unor proprietăți interne ale sistemelor de reglare automată (în speță stabilitatea asimptotică internă) pune în evidență diferența calitativă dintre punctul de vedere ■bternic și punctul de vedere circuitul referitor la studierea configurațiilor în reacție negativă Aceasta diferență de esență marchează în fond trecerea la analiza de tip intrare/iesiro (de tip circuitul) la analiza făcută cu ajuto-rul reprezentărilor de stare (de tip sistențic) , PERFORMANȚELE REGIMULUI DINAMIC Acestea sînt descrise prin indicii sintetici ce caracterizează răspunsul ^dicial ai sistemului cu funcția de transfer ^o(s)(^o(^)L ud mă al tomponen-, ci forțate a mărimii de ieșire y, atunci cînd mărimea de leieiință yr aie o tip treaptă imitai a Acești indici sintetici au fost introduși și (-finiți în paragraful și sînt * suprailmiărirea care ш acest context capătă denumirea de suprareglaj, indicele, de, oscilație % durata regimului tranzitoriu t , precum și perioada oscilațiilor T (aproximativ definită) și numărul de oscilații N (v și fig ) Evaluarea indicilor menționați s-a făcut în paragraful și nu o reluăm Această evaluare are la bază reprezentarea în pulsație a funcției de transfer în circuit închis a sistemului de reglare automată • Astfel, considerînd pentru început cazul sistemelor netede vom pune H (jto) — H (w) ei‘₽’( ) -F/H/w) ( ) unde Я (б>) A | H (jto) I, ) și respectiv F(co) introduse (generic) prin ( ) și ( )) Așadar, evaluările din paragraful sînt identice, cu referire însă, în acest context, la funcțiile ( (m) și E ( ), Ub(u) A Re Hft(ju), E (ca) = Im#ft(jo>) Dezideratul exprimat mai sus obligă în mod natural la explicitarea următoarelor dependențe de tip „circuit deschis/circuit închis": — determinarea lui ио(ы) (У (со)) din ад Și r ( ); — determinarea lui Н (ы) ( ), adică din locul de transfer Hb(jo>); — determinarea lui (ci>) (F ( ) și v Vb = ct" Determinarea lui Г/ (ы) (*zo( ))) caracteristicile logaritmice caracteristica amj)litudine-fa%ă Considerînd relația de bază ( ) obținem TT z Hb cos l tg To ( ) sau din -I- llb cos - ]Hb sin cpt, ( ) de unde, prin identificare, Щ -I- Hb cos și celălalt în zona UQ , dependența ilustrată prin curba din figura și riguro s adevărată pentru elementul de ordinul doi capătă o cvasivalabilitate generală Așadar, T se determină utilizînd curba din figura prin cunoașterea valorii de rezonanță HQrnax Aceasta se poate determina în două moduri: — prin utilizarea locului de transfer și a cercurilor de Ho = ct ilustrate în figura (v fig ); — prin utilizarea caracteristicii amplitudine-fază (dedusă direct din caracteristicile semilogaritmice) și a diagramei Nichols din figura (v fig ) Este vizibil că odată cu valoarea de rezonanță se poate cui imediat, pe locul de transfer sau pe caracteristica amplitudine-fază, și pulsația d? rezonanță or(li^rnax — Й ( , ( ) litiliHnd ( ) pentru explici tarea lui C în funcție de Ям “ obținem din ( ) și făcînd ЯОжіі„ = k / —' max~b~ x/ffi) max J «ОвевЬШшАв cor ' ^ max “ ^^ max " ' ( ) unde , » Aprecierea numărului de oscilații N Acest indice arc sens să fie introdus numai atunci cînd toate rădăcinile polinomului caracteristic /K(s) sînt reale (și evident negative) In acest саг N este dat de jumătate din numărul trecerilor prin zero ale componentei tranzitorii în consecință, N va fi cel mult egal cu grad/id polfy&m/ulwi caracteristic • Să vedem acum cum se transferă rezultatele de mai sus în cazul sistemelor discrete de reglare automată în acest sens, din relația ( ), care devine pentru X == x h&) = (ЗЛ ) + ИьѴ~) obținem punînd x = + w/\ — w Яо(«>)^ ( ) + Hfw) Așadar ntre echivalenții netezi de ordinul zero relația care definește corespondența de tip „circuit închis/circuit deschis" este identică cu aceea din cazul sistemelor netede în consecință punînd ^(h) ~ ^(t)) + P№) ( ) ej o(’n) ( Л dependențele elementelor care definesc relațiile ( ) de cele care definesc relațiile ( , ) sînt grafic explicitate prin diagramele construite mai sus Răspunsul indicial este însă descris cu ajutorul echivalentului neted de ordinul I- //e(w) al lui // (z), adică £,(w) ( ) și are expresia (vezi ) V( ; - f** l sin t arctg ) V ■'fP* ir f» — г -ггПгД J и ( ) Punînd, cu w *= $о(Іп) “ ЙоМ e*** — ^(> ) + j^dOl) ( ) obținem din ( ) legăturile între părțile recile și imaginare ale echivalenților de ordinul și respectiv zero (în circuit închis) ( ) ( ) ч r = ^о(’П) — Д( ) ( ) ( ) Vom căuta acum să evidențiem o procedură care să permită evaluarea parametrilor caracteristici și x ai trapezului reprezentativ al caracteristicii reale &o(tî) a echivalentului neted de ordinul în circuit închis (v fig ) direct din reprezentările ( ) ale echivalentului neted de ordinul zero în circuit deschis Această evaluare va consta din două etape: în prima etapă se determină parametrii caracteristici și x ai trapezului reprezentativ al caracteristicii reale / (t]) a echivalentului de ordinul zero în circuit închis (v figura ), iar în cea de a doua etapă, utilizînd într-o anumită manieră relațiile ( )—-( ), se determină parametrii ?) și x ilustrați în figura Prima etapă este procedural identică cu metoda utilizată în cazul aprecierii duratei regimului tranzitoriu la sistemele netede de reglare automată (v fig ) Pentru cea de-a doua etapă vom proceda astfeh Conside- J ’ ^ A >)o [Hs(s)] c dacă și numai dacă funcția de transfer a s₽ • erorii Яе(з) are în origine un zero cu ordinul de multiplicitate > p, adică avem I He(s) = Sv -Ge(s), V > p, Ge( ) finit ^ Avem atunci din ( ) cu ( ) s»Ge(s) relație care arată că un sistem neted de reglare automată este precis la refe-, Mnța yr(i) == j)-i „ cu p > , dacă și numai dacă funcția de transfer în Circuit descins Hb(s) are în origine un pol cu ordinul de multiplicitate > p adica are tipul p, avînd exprimarea Г- ! p, adică avem ' ф) = Așadar, un sistem zero cu Utilizînd din nou ( ) obținem din ( ) adică, un sistem discret de reglare automată este precis la referința yT(f) — = ^-i -Г P> , dacă si numai dacă funcția de transfer în circuit (# — )₽ deschis Дф) are în z = un pol cu ordinul de multiplicitate > p, adică avem Se vede că G*(l) este chiar factorul (discret) de amplificare al lui H^y} Qintetic, rezultatele de mai sus sînt ilustrate în figurile , a și b Cititorul poate obține direct rezultatele ele mai sus ca o consecință a teoremei fundamentale a reglării (v teorema ) S-a preferat însă o deducere intrinsecă a acestor rezultate pentru a familiariza cititorul cu raționamente care operează direct cu funcțiile reprezentative ale sistemelor de reglare automată O întrebare firească care se pune este următoarea: în ce condiții un sistem de reglare precis la referința poseda proprietatea (A), adică lim e(/) =- și pentru perturbatia ae același /~>o / z X tip v(t) ~~|ф) = jg/’ — j, cu același p > ca și în expresia jui S " l) V Fig ia AwtaixmW» ^l«r« antomatîi >re »^a пееешіі Ц HU/iciMlă M r z f / ' ' ) ' • l J • - ' j •’ i У ’ ■J f ^Y' rfj'MltA c făspunsoi h M un нЫет de re* s î jtaSifrft’ proprietatea (Ц) p гл în expresia hi AWt/ * M»(/) „ В / I / o (; ) «Я sero cu ontciiul de £шм&}^ ' «z pmU Va « «тисам ste «Wla f„ Z * P’ r„(X) =a ya(X)y>i(X) »"«c r!( și A(X) divid ,(X) și respt;t)iv /)(A)> Jl ^i(x) ra(X) unde [J • ’■ ■x(A) p unde punctul de Șea din figura r (x) A(x)A(a) »J ' ' p Punînd atunci ( ) va rezulta Яг(х) A rÂ, H (Ă) A A(X) PM д A(x) că un sistem de reglare care satisface ( ), precis — în cazul neted, sau — cP la referința în cazul discret, posedă pro- z prici alea (R) pentru v(t) = Л" — în cazul neted, sau == jg- -în cazul oP \O discret, dacă și numai dacă funcția de transfer Hff) posedă în origine гт pol cu ordinul de multiplicitate p în cazul neted, sau un pol în X = cu ordinul de multiplicitate >p, în cazul discret Așadar, structura din figura devine cea din figurile , a și b, unde P ( ) în cazul discret, adică dacă și numai dacă funcția de transfer a erorii sistemului neted de reglare posedă în origine un zero cu ordinul de multiplicitate v > p și respectiv, funcția de transfer a erorii sistemului discret de reglare posedă în z ~ I un zero cu ordinul de multiplicitate v > p, adică dacă și numai dacă sistemul de, reglare este precis la referințele ( ) in cazul neted și respectiv ( ) în cazul discret Numerele Я^( )/г! în cazul sistemelor netede de reglare și respectiv în cazul sistemelor discrete de reglare, se numesc coeficienții generalizați ai erorii în cazul cînd v = p, adică este verificată egalitatea ( ), mărimea nemția (/ ) ( ) P! în cazul sistemelor netede de reglare și e, A (Q) ( ) p! în cazul sistemelor discrete de reglare, poartă numele de eroare staționari:, (xie-nulă) ia referințele ( ) și respectiv ( ) Așadar, eroarea staționară nenulă va fi evaluată de prima derivată nenulă £Tp)(O)/p! în cazul neted și Я^р)( )/р! în cazul discret în cele ce urmează vom face permanent ipoteza că (cel puțin) Яе( ) = ( ) în cazul sistemelor netede de reglare și, Я ) (ЛП și deci e, lim e(/), ( ( ) І-фвО evaluare care se face direct cu ajutorul teoremei valorii finale și conduce la expresia ( ), expresie care arată că eroarea staționară (nenulă) este invers proporționala cu factorul de amplificare al funcției de transfer în circuit deschis Așadar, cu referire la specificațiile în circuit deschis ale sistemului de reglare automată, factorul de amplificare în circuit deschis este unicul element care determină eroarea staționară» Să facem acum o evaluare a erorii staționare în funcție de specificațiile în circuit închis ale sistemului de reglare automată, adică în funcție de elementele care caracterizează funcția de transfer HJis) (în circuit închis)» Vom face această evaluare reamintind că condițiile ( ) și ( ) sînt permanent adevărate Ne vom mărgini cu efectuarea calculelor la evaluarea erorii staționare Я,\ ) și respectiv Яе( ) Cum Я£(Х) = — Я (Х), rezultă din ( ) și ( ) că Ho(O) = ( ) \ W ' • • 'К A * ч și respectiv H (l) = ( ; cu m ) și ( ) se vede din ( ) și ( , ) ș& și A speci/icA complet Jtmpia de transfer în circuit închis Din ( , ) rezultă mărime cc ovnluruză eroarea staționară pentru o mărime de referință yr(t) t I -C ’ ; >- !»» ( ^» а« ісЛ pentru o mărime do referință de tip so ii sa rampă unitară în catul discret, expresia ( ) devine (exercițiu pentru cititor) ti tti zr;(i) ( ) ( ) \desca eroarea 'tuționară ( ), ( ) la mărime de referință de tip r; apă se notează cu c,„ adică л ( ITt( ) în cazul neted în cazul discret Cum polii și zerourile /?, și zt apar în perechi complex conjugate și cum Re /> V sînt reale și pozitive în — consecință, se poate trage concluzia din inspectarea expresiilor ( ) și ( ) ca: polii funcției de transfer în circuit închis conduc la o creștere a erorii staționare ev, pe cînd хегогегііе conduc la o scădere a acestei erori Exemplul Să se analizeze, sub raportul performanțelor discutate mai sus, sELe^-l de reglare automată cu funcția de transfer în circuit deschis ( + , s) Ht,(s) = s(l + lOs) ( + , s) ( + , s) f / Caracteristicile amplitudine-pulsație și fază-pulsație sînt ilustrate în figurile , a și b iar caracteristica amplitudine-fază, cu specificarea curbelor de L = Uomax Și t/ = — Goai; este ilustrată în figura Se obțin valorile Z/o( ) = [ ( , ) » , , = , L O( , ) e» , ) ‘ E®« A A ’ ' > U Fig SINTEZA CONVENȚIONALĂ A SISTEMELOR DE REGLARE AUTOMATĂ în paragraful s-a expus modul general de soluționare al problemei reglării în sensul sintetizării unui sistem compensator care să satisfacă dezideratele (S) și (A) Dacă clasa mărimilor exogene devine restrictivă (v ( ), ( )), adică mărimile externe sînt de tip polinomial (în timp), atunci condițiile ( S) de stabilitate și (A) de reglare, în particular în acest context, de precizie, se detaliază, în performanțele regimului dinamic adică în Performanțele răspunsului indiciul al lui / (s) (Я (з)) și respectiv în performanțele regimului staționar, măsurate de un singur indice sintetic și anume eroare staționară - Vom înțelege prin sinteza unui sistem de reglare automată rezolvarea problemei reglării în cazul restrictiv al mărimilor exogene polinomiale cu asigurarea suplimentară a unor performanțe ale regimului dinamic și staționar apriori impuse Mai exact, problematica sintezei convenționale constă în determinarea unui sistem compensator astfel îneît pe lîngă satisfacerea condițiilor obligatorii (S) și ( ?) să fie îndeplinită și următoarea listă (minimală) de performanțe ale regimului dinamic și a celui staționar u% (x) Atunci, așa cum dealtfel s-a explicitat și în ( ) se obține imediat că Я»(х) = Уо(Х) - Я (ЗЛ ) ь ] - ЯП(Х) jbq(Ă) unde rt(X)Ar (X) AW^^)-^) ( ) rezultînd vizibil că fracția гь(Х)/Л(Х) este dc asemenea ireductibili W unde, conform celor de’ mai sus, perechile de polinoame (r»(X), Л(>)) ? (r(X) jfr(X)) sînt coprime Punem m consecința />(X) = р(Х)Дх(Х), A>(X) = Л(Х)Дх(Х) ( ) r(X) = ?(Х)Дг(Х), гь(Х) = гй(Х)Д (Х) ( ) Unde ; Д,(Х) Ac m m d c {/>(Х), А(Х)}, Д^Х) = c m m d c {/(X), г (л)} ( ) Atunci din ( ), obținem cu ( ), ( ), zc(X) ? (X) + r(X)p(X)f|,( ) = р(Х)Д (Х) г(Х)Л(Х) + r(X) Д (Х) P(X) Fo(X) = P(X)r(X) [Ді(Х) p (X) + Д (Х)г (Х)] — PW ^(х) CiMx) + ro(x)] = ; (X) F(X)j (X) unde s-au utilizat ( ), ( ), și ( ) Cum evident (X) F(X)J c Așadar> pentru satisfacerea, condiției ( ) este necesar și suficient ca zerourile și polii funcției de transfer //(X) ncsimplificabile cu polii și respectiv zerourile funcției de transfer impuse în circuit deschis H,,(X) să fie în O", în cazul neted și în Ut( ) în cazul discret în consecință din punct de vedere procedural H (X) = / (>Ш(Х) trebuie aleasă astfel încît fracția >p(X)/(^ (X) — r (X)) să simplifice prin înmulțire cu fracția ^(X)/r(X) toți polii și zerourile acesteia din urmă care sînt plasați în Re X > , respectiv în | X | > în fine, pentru ca condiția ( ) să fie satisfăcută se vede că trebuie să avem, după cum rezultă din ( ), o - și deci rezultă că ; , p e > e - ( ) Deci, pentru îndeplinirea condiției ( ) este necesar și suficient ca excesul polilor față de zerouri ale funcției de transfer impuse H (s) să fie superior excesului funcției de transfer H(s} a sistemului Restricțiile ( ) și ( ) impuse ecuației algebrice ( ) se numesc restricțiile de structura ale sintezei exacte Aceste restricții se reflectă, așa cum am evidențiat mai sus, în alegerea lui în sensul că Д>(Х), pe linge d zi-deratele ( ), ( } și ( ) trebuie să satisfacă și restricțiile de struc tură Exemplul Fie sistemul neted cu funcția de transfer Ф s - p(s) H( ) « ireductibilă % ?i Л presupunem că cerințele ( ), ( ) șl ( ) sint satisfăcute de funcția de transfer circuit închis T-J (t) « P* t——za ~ ) ireductibilă ot' s» - I- « ('n»le așa cum s-a mai precizat / ( ) ~ ) Se vede că e = și ‘Mia ( ) este satisfăcută Se vede imediat că ■ ’• ' ? % «* consecință con' / ^( ) a Atunci rh(s) , *о(д) pb(s) /■«(*) ' ro(s) & * Pentru a simplifica zeroul $ «e d!n c(s} s(s)pb{s) s rezultînd structura din figura Trecem acum la alegerea propriu-zisă a funcției de transfer ЯС(Х) din punctul de vedere al satisfacerii performanțelor ( ), ( ) și ( ) Inspectînd expresiile ( ) și ( ) se vede că alegerea funcției de trar : -HQ(s) revine la precizarea locației polilor pț și a zerourilor a acestri funcții în planul complex O Din acest motiv, procedura de construcție a lui HQ(z) se mai jmxusște și metoda repartiției polilor și zerourilor funcției de transfer în circuit închis Această metodă se aplică relativ diferențial în cazul sistemelor netede și respectiv în cazul sistemelor discrete începem cu • Cazul sistemelor netede Vom presupune pentru început că b H(s) nu are poli și zerouri în O* Atunci, singura restricție de structură este e > Vom alege considera ca „element fundamental", elementul de ordinul , /ТТ /Л\ \ /Л У A I I Din ) și ( ) rezultă că pentru răspunsul indicial al acestui element indicii sintetici principali au evaluările: — £rr , * — suprareglaj ul: a =« e — durata regimului tranzitoriu: lt Atunci din ecuațiile F\g ’ id rezultă ^succesiv c, și apoi ,» In consecință, z va fi selectat astfel incit » z Introducerea zeroului -^-z poate avea însă două efecte neplăcute tn primul rînd, după cum se vede din ( ), excesul este acum ==• și în consecință s-ar putea ca restricția de structură ( ) să fie Dacă așa se întîmplă, va fi necesar ca în ( ) să se introducă un mentar —p, adică ( ) să devină ^^^В^ні^ВВ^^В î BB Bl f П /л I încălcată, pol supli- W (S) e - — — —— Ho == I O )(« -I- ^>„ + o> ) ț unde /> > o este astfel ales incit să nu altereze inegalitatea ( ), adică să avem în continuare (vezi ( )) i: e == ■—- trebuie astfel ales incit °mponenta tranzitorie pe care o introduce și anume de tipu e , e K stingă „repede", astfel încît performanțele dinamice ( ) și ( ) să nu fie alterate Aceasta revine, împreună cu ( ), la o condiție practică de tipvil /> = ( - ) „ Dacă condiția a) este infirmată, adică dacă de exemplu e = , atunci se va alege e$ — și în consecință evaluarea lui H (s) nu va începe cu o funcție de transfer de tipul ( ), ci cu una de tipul #o(s) — № (s+/>')(s + C«„s + « ) ( ) unde, după determinarea lui £ și wn din ( ), ( ), p' se determină ev o condiție identică cu ( ) Detalierile în continuare ale procedurii, în sensul verificării și satisfacerii inegalității ( ) în acest caz, sînt lăsate ca er er-cițiu cititorului O infirmare a condiției b) se tratează ca în exemplul Spre exemplificarea metodei discutate considerăm Exemplul Se dă sistemul cu funcția de tiansfer > și specificațiile ( ) Din ( ) obținem iar din ( ) este convenabilă valoarea wn = Obținem Va fi deci necesar să introducem un zero și corespunzător uu pol astfel încît se avem ou cu ’ ‘ f ѣ W «ț • F * v ii *• % > v • x л e f He(s) = H- (S) HM = ^ojA )(l+ Tu)(l^ s(s + s + ) • ‘ ‘ ) •£' • ♦ I Л c unde К și ' se înlocuiesc în funcție de parametrii procesului Trecem acum la • Cazul sistemelor discrete z> Se poate imagina și în acest caz, o metodă similară celei expuse în cazul sistemelor netede, bazată pe elementul discret de ordinul doi Nu vom expune un asemenea procedeu, lăsîndu- ca exercițiu cititorului și vom da o metodă eminamente specifică sistemelor discrete și anume metoda răspunsului indicial cu durată finită a regimului tranzitoriu Aceasta constă în sintetizarea unei funcții de transfer H (z), care să satisfacă restricțiile de structură ale sintezei exacte și care să producă o funcție de transfer a erorii He(z) = — — Яо(^), astfel îneît eroarea corespunzătoare unei mărimi de referință în treaptă (unitară) și anume е(г) = HĂZ) ( ) r — să șe stingă într-un număr finit N de pași, adică ^ = ( ) Așadar e(z) = + e(l)? + г( )г~ + + e(N - ) z^ ( ) (v fig ), unde evident e( ) = I Din ( ) și ( ) obținem H ( } — [!— У(АГ— )]+[з/(-гѵ'— I)—y(N— ) | ?+•• • +țy( ) ) iLLtMLLl'lt ( ) ( ) • • • u nde У(і) I t f e să observăm că prin discretizarea sistemului neted (Л, b, cT) obținem sistemul capitolul anterior, observăm că pr discret LL> — , unde după cum se cunoaște dm елг& Deci ( ) și în tensednță, generic vom avea conform cu ( ) Cum evident H(») = cjba va rezulta din ( ) și ( ) că dud sj ( ) adică generic, prin discretizarea sistemelor netede, excesul devine unitar Inspectînd ( ) se vede că eQ = și deci se poate considera că condiția de exces ( ) este automat satisfăcută în ceea ce privește condiția ( ), ce explicitează restricția de stabilitate internă, aceasta revine la impunerea unor condiții de tipul ( ) unde ;ȘÎ zt sînt P°iii Și zerourile lui H(z) plasate pe frontiera și în exteriorul cercului unitate Condițiile ( ) și ( ) impun anumite restricții între mărimile y(i), restricții cc se explicitează în fiecare caz concret b, c Interpretarea acestor condiții este elementară și nu o mai comentăm, menționînd numai faptul că condițiile b) și c) se iau în considerare numai după satisfacerea restricțiilor de structură Să constatăm în încheiere că specific metodei expuse este alocarea în origine a tuturor polilor funcției de transfer în circuit închis HxtinyJuI , Dîndn-bc «sternul cu funcția de transfer se cere să se sintetizeze un regulator J/C(i) astfel incit eroarea t ec stingă fotT-un număr minim de puți Se vede că />+ = - și deci vom avea o singură condiție de tipul ( ) Namănd minim A cu care acesta poate fi satisfăcuta este AT = Arcuit deschis începem cu I METODA RAZATĂ PE CARACTERISTICILE SEMILOGARITMICE Această metodă se compune din mai multe etape: Etapa , constă din transpunerea „în forma" părții reale f («) a condiției de precizie la y>(t) = l(t) și a performanțelor ( ) și ( ) ale regimului dinamic Avem astfel: a Condiția de precizie se transpune banal prin Ct( ) = ( ) b CW’jffa de suprareglaj о?/б Ug max egalitate ce are sens numai pentru ( ) min min și obținem din ( ) și ( ) max max л a - a , coeficientul de suprareglaj , , , > , ( ) wo to ^ Atunci, cu valorile cod, w , w (, fixate relativ la (>)(( conform cu ( , ), durata relativa in raport cu co a regimului tranzitoriu a răspunsului indici I obținut dm sumare a algebrica a răspunsurilor indiciale corespunzătoare celor două trapeze va fi parametrizată în raport cu (Jomiw și va avea o de pendența (obținută prin calcul) de tipul celei ilustrate în figura , în consecința, odata Uq„Mx determinat din considerente de suprareglaj acesta determină, prin intermediul dependenței, din figura , , în mod univoc durata relativă т( a regimului tranzitoriu și deci obținem «o ( ) Etapa a -a constă clin transpunerea în caracteristicile semilogarit-mice ale funcției de transfer în circuit deschis Hb(s) a performanței regimului staționar explicitată prin condiția ( ) și a „formei" caracteristice! reale de frecvență Î (co) caracterizată prin parametrii sintetici £ o(O) = l, și Qo obținuți în prima etapă Această operație de transpunere se face practic prin stabilirea zonelor de joasă, medie și înaltă frecvență a caracteristicii asimptotice amplitudine-frecvență [Нй(со)]аБ începem cu a Zona de joasa frecvență, a cărei determinare revine în fond la precizarea poziției asimptotei de joasă frecvență Lb ( S) unde Kb este amplificarea în circuit deschis, iar q este tipul funcției de transfer Hb(s) în consecință, pentru ca sistemul de reglare să fie precis la v, (ri ~ == Ф) se va alege lar condiția ( ) va furniza conform cu ( ) șî ( ) b , b Zona de medie frecvență este delimitată de trei elemente: pulsația de tăiere a caracteristicii asimptotice aniplitudînc-pulsație în circuit de>Uus panta acesteia și așa-numita rezervă de am flit adine a cărei semnHuațte \a ti imediat explicata, ’ Pentru determinarea pulsației de tăiere să observăm ca acrasra со rcspunde intersecției locului de transfer Z^(jw), cu cercul Unitate, adica cu cercul de / , Așadar, ) OO, b'o max (U max ) (^ m«x — O ( ) ! OO = Uq max Atunci se vede imediat că C ПМ Ѵ I OO I I ( ) max rezultînd încadrarea unghiulară a celor două cercuri din figura Mai mult i + I oo, | та v TT — max ( ) J ma s Introducem următorii indicatori L A log | | = - log | OPZ | = log inax S H ОС —-я Uq max max ( ) denumiți rezervă de amplitudine și respectiv rezervă de fază Este acum evident că satisfacerea condiției ( ) (și în mod obligator a celei de stabilitate internă), conduce la o poziționare a locului de transfer Hb (] ), ilustrată în figura (vezi locurile trasate cu linie continuă și întreruptă) și care practic revine la interdicția pătrunderii acestuia în trapezul curbiliniu hașurat Vom caracteriza în consecință zona de medie frecvență prin satisfacerea simultană a inegalităților — E фь(«) > — тс + a ( ) în fine, c Zona, de înalta frecvență este caracterizată printr-o relație între pantele asimptotelor de înaltă frecvență ale funcției de transfer a sistemului de reglare în circuit deschis Hfs) și respectiv ale funcției de transfer H(s) a sistemului (procesului), de tipul eb -= — c b dec dec e — dec ceea ce implică evident ( ) ( ) — e> fiind astfel îndeplinită restricția de structură privind excesul ceea ce темпе la obținerea unei funcții de transfer proprii (strict propui) lîc( ) a compensatorului Cele relatate mai sus referitor la cele trei zone reprezentative de frecvență sînt ilustrate în figura Ixgălura între zona de joasă frecvență și cea de medie frecvență stface în mod uzual printr-un segment de dreapta cu panta de l ^ / ec sau - dB/dec sau, eventual printr-un contur poligonal mai complicai Legătura între zona, de medie frecvență și cea de înaltă frecvență ж realizează de asemenea printr-o linie poligonală cu pantele multiplu de- dB/dec Odată caracteristica [Я»]ю obținută, caracteristica [Я,(«) ш rezulta din diferența (grafică) [Hc(w)]dB = [#»( )] porțuuu a limită, a condiției ( ) se poate accepta ca între panta asimptotei la multă cy, caracteristicii [Яь(ш)]ив și cea corespunzătoare noi caracteristici [Н(ю)] С + s) ( -s) — К s “ s( + , s) ( + , o-) care este funcția de transfer căutată Răspunsul indicial al sistemului de reglare proiectat este ilustrat în figura Se coa* stată o uș Metoda are la bază procedura de rezolvare a următoarei probleme: dîndu-se funcția de transfer în circuit deschis explicitată în forma ( ) unde Hb(s) este cunoscută, să se determine Kb astfel încît valoarea de rezonanță HQmax a funcției de transfer în circuit închis să fie egală cu o valoare prescrisă/ Presupunînd problema rezolvată, aceasta corespunde unei poziții a hodografului Hb(]a) ilustrată în figura , unde cercul trasat cu linie continuă și tangent la locul de transfer este un cerc din familia de cercuri До = ct Ducînd din origine dreapta Д tangentă la cercul Ьо ~ НОтлх, acesta face cu semidreapta negativă a axei reale un unghi a, explicitat prin t' ! unde s-au utilizat relațiile ( ) Considerînd acum o omotetie cu centrul hi O și factorul a tandemului format din locul de transfer Я (іы) și cercul Яо Н н,, , se obține tandemul omotetie trasat cu linie întreruptă format din locul de transfer cunoscut Йь(у*) și cercul C cu centrul O{, plasat pe axa reală și tangent simultan la locul de transfer Яь(] ' Il(s) = H(s) e-^, H(s) e!E« (s) ( ) Trecem acum la III APLICAREA „METODELOR ÎN FRECVENȚĂ" LA SINTEZA SISTEMELOR DISCRETE DE REGLARE AUTOMATĂ Rezultatele obținute în cazul metodei I bazată pe caracteristicile semilo-gariunice ale sistemelor netede pot fi aplicate, cu unele modificări, și în cazul sintezei sistemelor discrete de reglare Astfel, inspectînd diagramele din figurile și , rezultă că pentru un trapez unitar U (p), parametrii caracteristici și x care furnizează răspunsuri indiciale acceptabile sînt: T) e[ , ; ], xg[ ; , ] ( ) Practic, zona de lucru pentru parametrul ( ) este mai restrînsă t)og[O, ; ] xg[ , ; , ] ( ) in#' ipld i adică t) și x vor fi plasați în jurul unității și respectiv în jurul valorii , Din diagramele ilustrate în figurile , a, b și c rezultă că zona de lucru menționată este aproximativ aceeași și pentru trapezul în consecință, transpunerea performanțelor în caracteristicile semilogaritmice ale funcției de transfer în circuit deschis (vezi ( ) și ( )) va fi practic similară cazului neted Astfel: — rezerva de amplitudine se alege L = — dB, cu o micșorare a — L pînă la eventual — dB, în scopul reducerii valorii UQmin și — zona de medie frecvență este impusă de pulsația de taiere situata plicitată de o relație similară cu ( э), practic utili- valorii L^ (oo), ultima valoare nefiind nulă, după cum este cunoscutj în jurul unității și ex zîndu-se o formulă de tipul ( ) și anume — , ( ) / (după cum rezultă din ( ) pentru иОтп,ѵ~ , si x~(,l), —'întrucît Heste o funcție de transfer (ptopiie) de faza neminima (Z/Jl) = ), aceasta obligă la introducerea în v) = a unei variații de panta de + dB/dec, si a unei variații de fază de —tc/ , Din acest motiv controlul simultan al rezervei de amplitudine și al rezervei de țața este obligatoriu, teorema lui Bode nefuncționînd w л , , în încheiere, se impune o precizare Se constata imediat ca ш sistemelor discrete de reglare, intervalul de medie frecvența oști apriori impus (v ), ceea ce nu se întîmplă în cazul sistemelor de reglare netede O puma aplicație, de ansamblu, este aceea că în cazul transformării Q nu ex sta o teoremă a asemănării similară celei din cazul transformam LaPla“ Л consecință o dilatare (comprimare) a frecvențelor nu implica o compi in (dilatare), cu același factor, a timpului O a doua expica,ie, e nua ț transformării Laplace și m mult procedurală, este următoarea în discretizaica unui sistem neted pasul h se poate găsi în două zone O primă zonă, cea naturală, este aceea a valorilor mult inferioare constantelor de timp ale procesului continuu; o a doua zonă, impusă de restricțiile privind rapiditatea de calcul a compensatoarelor (calculatoarelor) numerice, este aceea a valorilor lui h comparabile cu constantele de timp ale procesului Fie, din punctul de vedere al acestui ultim caz, o pulsație de discretizare co = comparabilă cu pulsația de tăiere a procesului continuu, să zicem =• Atunci, din egalitățile evidente z = = = - t)/i — у) = ei arct(^ rezultă = tg ~^și din ( ) se obține cu иЦйЙIкВ о = —valoarea r = l, tg-^- = , adică, o valoare comparabilă Л cu unitatea și cuprinsă în intervalul de valori uzuale furnizat de ( ) Așadar, din punctul de vedere strict aplicativ, o proiectare a sistemelor discrete de reglare efectuată cu metode eminamente specifice sistemelor discrete apare ca indispensabilă în cazurile în care discretizarea sistemului (neted) reglat se face cu un pas de discretizare comparabil cu constantele de timp ale acestuia, în acest caz, orice similitudine cu funcționarea sistemelor netede impusă de „micimea" pasului de discretizare încetează METODA DE SINTEZĂ BAZATĂ PE LOCUL GEOMETRIC AL RĂDĂCINILOR funcția de transfer în circuit deschis a unui sistem de reglare automată pe care o scriem în forma m ( ) Atunci cu ipotezele de inductibilitate cunoscute n in Xo(M — l’î (ь - A) + П ~ ~?) ( ) este chiar polinomul caracteristic al sistemului rezultant în circuit închis, iar lZo(^) l =* -I- К ( ) Vom numi loc geometric al rădăcinilor locul geometric descris de, zerourile poUnomului caracteristic %q(X) explicita! prin ( ), atunci cînd factorul К parcurge intervalul [ , oo) Este vizibil că acest loc geometric conține n ramuri simetrice față de axa reală și caracterizate simultan după cum rezultă din ( ), prin: m n П IĂ — A t = * Ш n == ( k + ) t, keZ i=l i— unde ( ) ( ) ( ) Vom începe cu Regulile de trasare ale locului rădăcinilor Punctele de plecare și de sosire, Făcînd К = în ecuația / (X) = , rezultă din ( ) că aceasta se reduce la П (X -A) = ( ) t=i și deci locul rădăcinilor pleacă din polii pif i = , , , n ai funcției de transfer în circuit deschis, Făcînd acum în aceeași ecuație К oo se obține, pentru rădăcinile cu locație finită X П (X - A) = o ( - > î=i adică m ramuri ale locului ajung pentru К со, în zerourile funcției de transfer în circuit deschis Restul de n—m ramuri sînt în consecință ramuri asimptotice Comportarea ramurilor asimptotice Sciiem acum ecuația /o(^) d în forma Atunci cu substituția ,u = X - ca ( ) ecuația ( ) devine p ’»-” = - К ( ) și deci, acele rădăcini Xk ale ecuației x (X) = O care au comportarea | Xk | —> со cu К —> со vor fi furnizate, conform cu ( ) și ( ), de = , , n — m — ( ) Așadar, direcțiile asimptotice ale acelor ramuri ale locului care tind ia infinit (cu К —■» oo) fac cu axa reală unghiurile (v ( )) O, , n — m — ( ) și trec toate prin punctul ( = - ± J Atunci Punctele de plecare sînt polii i = , , , Neexistînd zerouri, toate ramurue lo cuini sînt asimptotice Direcțiile asimptotice sînt date de adică iar Fig «sie punctul de intersecție al direcțiilor asimptotice, Exista o singura porțiuno a locului plasata pe axa reala și anume intervalul [ — , ) Intersecțiile cu axa reala sînt date de care are o singură rădăcinii reală a = — e[— , ] și deci este o intersecție cu axa reală » Pentru intersecțiile cu axa imaginară explîcltăm XoW» carc estc Xo(X) « X , % = > , % = - К = rezultînd K= / Introducem această valoare a lui К în ecuația caracteristică și facem X= jo rezultînd (o - co -}- / ) - jw(- o - ) = Deci rezultă со — со' ~ / = care nu are rădăcini reale și co — = (со ^ ) care are ca rădăcini wb = ± I/ — și deci sînt punctele de intersecție cu axa imaginară Unghiurile de plecare sînt, cititorul verificînd aceasta, Se obține astfel locul gecmetric ilustrat în figura Aplicarea metodei locului rădăcinilor în sinteza sistemelor netede de reglare Operația de sinteză constă dintr-o succesiune de operații elementare de compensare utilizînd două tipuri standard de elemente compensatoare: a compensatorul cu avans de faza avînd funcția de transfer , Compensatorul cu întârziere defaza, avînd funcția de transfer ВД — Justificarea denumirilor celor două tipuri unde, în ambele cazuri, zc, pc unde pA (”*£ + î Vi — c = = iKl + tg ( ) Din motive „constructive" se impun unele restricții asupra raportului IAI / ІгсІ> acesta luînd valoarea maximă (intervalul [|zj, |^c|] în figura , « fiind de cel mult o decadă); în consecință pot fi necesare, pentru îndeplinirea condiției ( ) compatibilă cu restricția impusă raportului \pc\l\z'c\, utilizarea mai multor compensatoare cu avans de fază conectate în cascada, adica în fond iterarea atît cît este necesar a acestei etape kSe evaluează acum gradarea în factorul К a locului corespunzătoare polului pd, «adică • , m П (Pa - zi) + К -= ( ) ГІ (Pd Pil = unde z și pc au fost absorbiți în zerourile și respectiv polii sistemului în circuit deschis Odată factorul К evaluat din ( ), se introduce în ( ) și se calculează factorul de amplificare Kb în circuit deschis, adică Kb = [s^ (s)%= ( ) (sistemul fiind de tipul q = ) Dacă acum /Kb > zvd, se trece la etapa a doua, de corectare a performanței asociate regimului staționar în caz contrar, compensarea nu este necesară Etapa a -a constă din introducerea unui compensator сгі întârziere de fază ( ), astfel încît să avem îndeplinită condiția Kp pc %c și în același timp performanțele dinamice să na fie alterate, adică ramurile complexe ale locului obținut în etapa , să rămînă practic nemodificate Aceasta se realizează plasînd perechea (zc, pc) „aproape de origine'" ceea ce se realizează prin situarea lui zc, astfel încît distanța acestuia pînă la origine să fie văzută din polul pd sub un unghi de ° (v fig ) Odată zc fixat, pc rezultă imediat din ( ) C APITOLU L TEORIA GENERALĂ A SISTEMELOR LINIARE Spre deosebire de capitolele precedente, în care a predominat punctul de vedere frecvențional în investigarea sistemelor liniare cu o intrare și o ieșire, în prezentul capitol 'punctul de vedele structural va fi preponderent Această schimbare de optică este impusă de însăși tipologia obiectului investigat care, dacă, în primele două capitole anterioare era relativ simplu, fiind vorba de sisteme cu o intrare și o ieșire, problematica dezvoltîndu-se în jurul comportamentului intrare/ieșire, în acest capitol obiectul investigației îl formează sistemele cu mai multe intrări și mai multe ieșiri Este firesc ca la astfel de sisteme metodologia de abordare să se modifice datorită în primul rînd „complexității interne" evidențiată de „interacțiunile" mărimilor de comandă, în consecință, descrierea acestor interacțiuni devine o necesitate obligatorie, -orientară față de studierea comportamentului dinamic propriu-zis, reclamînd cu precădere un instrument de investigație algebric și geometric Această modalitate de abordare a făcut în fond posibilă eficientizarea punctului de vedere sistemic, ca punct de vedere eminamente metodologic, șioperînd cu concepte specifice noțiunii de model, noțiune concepută să reflecte tocmai acele proprietăți sistemice ale obiectitlui (procesului) în acest sens, prezentul capitol introduce principalele concepte și proprietăți structurale ale sistemelor liniare, cu intenția operaționalizării acestora, dezvoltările făcute avînd ca trăsătură dominantă evidențierea unor elemente algebric invariante Să încheiem menționînd că asimilarea conceptelor și ideilor din prezentul capitol este esențială pentru înțelegerea problematicii sistemice și a modalităților de rezolvare a acesteia SISTEME LINIARE CONSIDERAȚII GENERALE Elementele primare cu care operează definiția care urmează sînt: mulțimea momentelor de timp, asimilată cu mulțimea numerelor reale M-sau cu mulțimea numerelor întregi Ж; clasa comenzilor U descrisă de aplicațiile continue (continue pe porțiuni) de la ML în № sau de aplicațiile de la Ж în Ж Se utilizează și notația ; — spațiul stărilor Qj asimilat cu ’ în — clasa ieșirilor Y i ? sau prin aplicații de la Ж în Introducem descrisă prin aplicații cel puțin continue de la Este utilizată și notația ф A j ▲ Definiția Se numește sistem liniar un triplet (А, В С) cu A pTRE»** fi eJESJi Ktn și C ES pxn ce se explicitează prin unde xr = Ax у Bu у = Cx x( ) = xQ ( ( ) unde ?/(•) g U este comanda, x(t) g Ж este starea, iar v(*) еУ este ieșirea Dacă t g JSL sistemul se numește neted sau continuat sau cu timp continuu, iar dacă t Ж, sistemul se numește discret sau cu timp discret Așadar, din ( ) și ( ) sistemele netede se descriu prin x(t) = Ax(t) + Bu{t), x( ) = x y(t^ = C*x(/), t g JLJt iar cele discrete prin • J • ч £ ж x(t -f- ) = Ax(t) + %( ) = x y(t) = Cx(t) /eN ( } Introducem: matricea fundamentală Ф prin ( > matricea pondere prin | СФ( В, Z eH A j f , СФД - ) В, ( - > Precum și matricea de răspuns cauzal Ia impuls ГД/) A !(/) Г(і), t e » sau t ( > Numim matrice de transfer, matricea din Wxw' (X) definită prin T(X) = Я Tc{t), ^{CW} t e ES i еЖ s pt sisteme netede z pt sisteme discrete Ă = n ) și ( ), ( ) se explicitează unic în T(X) = С(ХІ - А)- В ( ) /I x( $) - -Bk(s) ( , sx(s) — -t'o •-= j(s) = Cx(s) și respectiv »#($) — zx = Ax(z) - Bu(z) ( ) y(z) = Cx(z) Rezultă atunci cu ( ) și prin eliminarea lui x(s) și x(z) dm ( ; și respectiv ( ) M уЩ = СФМх +ТѴ)и(\) te reflexivă, simetrică si tranzitivă și deci definește efectiv o relație de echivalență notată ( , В, C) Д C) Din ( ) obținem imediat că Ф(/) == 'ф(/) - 'Г(Х) = T(X) ( ) ( ) unde Ф(/) și T(X) corespund sistemului (А, В, C) Semnificația relațiilor ( ) se obține direct din ( ) prin aplicarea schimbării de coordonate A Definiția Două sisteme (А, В, C) și (А, В, C) se valențe intrarejieșire dacă au aceeași matrice de transfer, adică T(X) = C(X - Ay'B = C(\I - A)~'B = Î(k) ( ) numesc echi- ( ) Se observă că și ( ) este reflexivă, simetrică și tranzitivă și deci este efectiv o relație de echivalență adesea notată și (А, В, С) ~ (А, В, C) Comparînd ( ) cu ( ), se vede că echivalența pe stare este automat o echivalență intrare/ieșire, reciproca nefiind valabilă STABILITATEA SISTEMELOR LINIARE A Definiția Un sistem (А, В, C) se numește (intern asimptotic) stabil-dacă lim Ф(/) = /-> oo ( '> Conform celor din capitolul , condiția ( ') este echivalenta cu a( ) cz C " pentru sistemele t i( ) pentru sistemele netede discrete Semnificația stabilității (asimptotice interne) este dată de Teorema Un sistem (Л, В, C) este stabil dacă și numai dacă lim Ф(/) x = Vv e (%> Demonstrație Condiția este necesară, întrucât ( ') implică direct ( ) Reciproc Presupunem ( ) adevărată și să presupunem că ( ') este infirmată în consecință, conform сц ( ) există oț li, asttel încît KeX> (|X|> ) Fie x un vector propriu asociat lui X Atunci din л X v rezultă imediat că te U t Atunci din ( ) punînd x = x / O (deoarece este vector propriu) rezultă că ( ) nu poate fi îndeplinită, contradicție care demonstrează teorema Să punem acum în evidență cîteva aspecte structurale legate de stabilitate Fie u polinomul minimal al matricei A și fie O“ și O , astfel îneît C® n П C‘ = și O = C“ U O" Fie ji(X) = (jl“(X) p ft(X) ( ) •o factorizare a polinomului minimal astfel îneît [[x“(X)] a O“, ElA*-)] c ( - ) Я? A Ker |і“(Л), O>b Ker f?(A) ( ) Atunci Intr-adevăr, pentru (i) avem din ( ) xe( a => tia(A')x = => Ap a(A)x = = => u®(- ) Ax = => Axe a (identic pentru Об”) Pentru a demonstra condiția (ii) să arătăm că “ П = { }- Q caz contrar, există x astfel îneît x e și x b, adică (X®(A)x = și [л (А)ж = ( ) Fie а (X) anulatorul minimal al lui x Atunci, din ( ) rezultă că ^(X) | p “(X) și ^(X) | j?(X) ( ) Cum |i“(X) și (?(X) sînt coprime (v ( )), din ( ), rezultă că |лж(Х) = ct^O și deci ct x = , ceea ce este absurd Așadar, se poate scrie * ( ) Fie un %g J • ( ) unde / colectează tonte valorile proprii ale lui A situate în C®, iar Ab pe cele din O‘ Intr-adevăr, fie și xn^ , astfel încît Aaxa = \axa Atunci avem evident ( ) A dim și de asemenea Ax = \ax rezultînd de aici că x® este un anulator al lui x fiind totodată și minimal deoarece XMX)] = Deci X — Xa|ti“(X) și în consecință X® eC“ Așadar, x t-rao ( ) (cu /eH, sau Ы în funcție de tipul sistemului), adică o traiectorie tinde către origine dacă și numai dacă este inițializată în Q ~', de asemenea, orice traiectorie inițializată în rămîne pentru orice/> în în fine, pentru orice evoluție liberă x(/) = Ф(/) г(), /> , se obține descompunerea unică x{i) = x~(t) + x+(f) cu x*(t) eț&F, /> și unde x (ț) -> cu t-> oo, iar x*(t) -/■+ cu t -* oo ' Din relația ( ) detaliată prin discuția dc mai sus, rezultă că un '- кіет este stabil dacă și numai dacă Qj -• Ж sau echivalent, % ‘‘ { ( Adesea se numește subspațiul instabil iar Ж ’ subspațiul stabil, denumiri provenite din considerațiile făcute mai sus PROPRIETĂȚI STRUCTURALE Introducem acum principalele proprietăți structurale ale sistemelor, proprietăți pe baza cărora se construiesc toate dezvoltă; ile ulterioare și care condiționează esențial întreaga problematică sistemică sub raportul solvă-bilității acesteia CONTROL ABILITAI EA x a Această proprietate operează cu perechea (Л, B) și vom începe abordarea ei pentru cazul sistemelor discrete, dinamica asociată perechii în discuție fiind descrisă de prima ecuație din ( ), Tranziția unei stări xQ după k pași, în urma aplicării comenzii (secvenței de comandă) u( ), ^( ), , u(k — ) se obține din ( ) cu ( ) (pentru cazul sistemelor discrete) și este %(£) = Л% + Л'^ВЦО) + Л*~ ВЦ ) + + Bu(k - ) ( ) Introducem matricea Rr^[B А В Л*-іВ] ( ) numită de k — control abilii ale, precum și vectorul de comanda în pas ~u(k - )" u(k — ) ( ) ceea ce permite rescrierea lui ( ) în forma compactă x{k) = AkxQ + Introducem ( ) A Definiția O stare există un vector de comandă = în starea x(k) = x x e (X) se numește controlabilă în k > pași daca în k pași, uk, care transferă originea, adică starea Avem imediat e Teorema O stare x g Qj este controlabilă în к pași dacă si numai daca x unde A Im Rk ( ) ( ) Demonstrație Făcînd în ( ) x = se obține ecuația = x ( ) care are o soluție, uk, dacă și numai dacă apartenența ( ) este îndeplinită Explicitînd acum ( ) rezultă, cu ( ) că / Evident pentru k — incluziunea este adeva- > rată Dacă acum * Cele de mai sus se pot prezenta compact și în forma unde iar ( ) ( ) n/W A ‘ ; °№ech* (Л' se numește controlabilă daca întreg spațiul « este controlabil, adică dacă Din ( ), ( ), ( ) și direct din teorema rezultă O pereche (A, B) este controlabilă dacă și numai dacă ran£ R = n ( ) e Teorema sau, echivalent, matricea R este epică Vom spune că sistemul (А, В, C) este este controlabilă Fie acum două sisteme echivalente ( ) obținem controlabil dacă, perechea (A, B) , C) Utilizînd relație care explicitează conservarea prin echivalență a controlabilitații Fie (A, B, C) un sistem și subspațiul său controlabil Punem Ж = ® ( - S) • ' / • » » •: Utilizînd teorema de descompunere controlabilă, avem rang [X/ — A B] = rang [)J — A B]= = rang = rang [X/ — A j ВЛ d~ rang [XZ — A$>] ( ) Dar a(A) == a(A ) ♦ а(Л ) și, pe de altă parte, conform condiției de necesitate, rang [X/ — Aj Bf ~ nY V X o-(A ) Rezultă din ( ) atunci ca rang [XZ — Л ] ~ n — n, — ѣ,, VX a(/l ) ceea ce este abursd dacă « > De observat că oricare ar fi perechea ( , B), condiția de rang din ( ) este automat îndeplinită VXeO\n(/l) Utilizînd criteriul lui Hautus se ohtinp nn A perechile controlabile în care matricea В este momeă дД^ТЛ^г privind • Teorema Fie perechea (/ , B) (i) există o transformare este monică Avem astfel сгь В monică Atunci: ccii'c aduce perechea la forma m n — тп з } n — m m J cstfconFolabUă eSt unlf •••> ^nnî) ( ) rezultînd prin „citirea" pe linii a matricei A și apoi a matricei B Cu această ultimă precizare, fiecare pereche (z , B) se identif ică în mod unic cu un punct G t?’+ew și reciproc Evidențiem coloanele matricei В scriind * • f\ V » » momente de timp este, după cum rezultă din ( ) Л'( Bît(-v) dr ( ) Introducem, * Definiția O stare xefâ se numește controlabilă la momentul t> dacă exista o comanda u(*) £ U care transfera starea Xq = , adica originea, în starea x(i) = x > l • * > -J J И î И? • • i \ Л i • -K • / V ' • • * * , L V • f» Făcînd în ( ) x = și %(/) = x, obținem ecuația \ ■ î • • * • Ж j I ’ * c C eA(^ Bu(x) dv= x ( ) Jo Pentru studierea ecuației ( ) introducem (pentru £> ) matricea -i t* ♦ ’ *V- Л' '' ) \ $ A V* ’ '** V ‘ ' ' '’ ' - t * CJ J h Ж( дГ eA(i-T) BBT &AT(i-x} dT ( ) *’o Din ( ), se vede că F (Z) este simetrică Mai mult, pentru Vx Ж avem W{t) x = C i'O BToAT^ x dv> ( ) adică W(r) este pozitiv semidefinită W(t)^ Caracterizarea algebrică a matricei W(t) este dată de • Teorema Pentnb orice t > (de fapt orice t ) avem ІтІУ(/)=бе / ( MO) (unde (Ц are semnificația dată de ( ), ( )) Demonstrație Deoarece W(t) este simetrică Im Г/( J Ker IP( și Ж = Im PF(Z) © Ker IF(/) în consecință, ( ) este echivalentă cu Ker ИДО => Ker RT ( Ml) și vom demonstra ( ), Fie я Ker И /) Avem lanțul de implicații duble (v și ( )) Иф) x = rz>F(/) x == k || Вт cAT^) x Ц dT Jo В елТ^”т) лѵ= , cu Вт елГ° = , си О Вгсл'‘йх — ѵО е ) Din ultima egalitate obținută în ( ) rezultă că — [Вгс- Г г']ѳ= = , k = , , , ( ) ceea ce conduce la xTAkB = ? = , , , ( ) și în particular la xTR = , adică x e Ker RT Deci : Ker ( ) c Ker ( ) Fie acum леКегѴ, adică xTR = și deci xTAkB = e = , , и — ( ) • * • • * • • • • v ?’* v • x • o stare x este controlabilă la momentul t dacă și mimai dacă ' ■ ■ • • Demonstrație ' ' ; ‘ ’ Condiția este necesară într-adevăr, fie un t> și x o stare controlabilă la acest moment Vom arăta că x e Im W(t), ceea ce este echivalent, conform teoremei , cu ( ) Conform ipotezei există o comandă U pentru care ( ) este adevărată Întrucît Im W(t) și Ker ( ) sînt o descompunere ortogonală a lui X, se poate scrie unic pentru orice x e Ж x=Xi + X i %leImP ( ), x e Ker ( ) ( ) Cum xr I x , se obține din ( ) |x i| =;vfx ( ) Atunci din ( ) rezultă cu ( ) că ел(/-т) J?»(t) d-r ( ) Dar x, e Ker ( și deci ( * = O, de unde cu ( ) т) S^(t) = O cu O și deci înlocuind în ( ) obținem || t j| = Atunci, din ( ), rezultă că xe Im ( - ' v Condiția este suficientă, intr-adevăr, presupunem ca o stare x satisface ( ) și fie un t>Q Atunci, conform teoremei ( ) xg Im W(t) ( ) Din ( ) rezultă că există ze^, astfel îneît x=W(t)z ( ) Introducînd ( ) în ( ) obținem x = ( e-W-t} B[BT елТѴ—д z] dr ( ) • Comparînd ( ) cu ( ), rezultă că comanda • ад(-г) Автелг r iz Ъ * » L V »"*’•! « - ♦ , • »’’> *' • * T — W Din teorema , rezultă că dacă o stare x este controlabilă la momentul t> eă este controlabilă la orâi? moment £> Introducem acum A Definiția O stare se numește controlabilă dacă existat astfel incit x să fie controlabilă la acest moment Din teorema și din remarca făcută la această teoremă obținem rezultatul final dat de • Teorema O stare x este controlabilă dacă și mimai dacă ( este adevărată Dacă ( ) este adevărată, starea x este controlabilă la orice moment t > , o comandă fiind furnizată de ( ) și ( ) Cu acest ultim rezultat, se constată o perfectă identitate algebrică și geometrică între noțiunea de controlabilitate introdusă în cazul sistemelor discrete și cea introdusă în cazul sistemelor netede, dezvoltările făcute începînd cu definiția încolo fiind perfect identice și în cazul sistemelor netede O deosebire totuși există din punctul de vedere dinamic și anume dacă în cazul discret transferul originii într-o stare r^ controlabilă se face cel puțin într-un pas (și cel mult în pași), în cazul neted, acest transfer se face oricît de repede Față de această ultimă afirmație, să considerăm ecuația în algebra cu convoluție Ф', (*) Ax - ( ) unde este coloana i din B, iar este distribuția Dirac Ecuația se mai scrie și (*> ( - Л) *x = bt avînd soluția (*) (*) • x = ( - Л)- * ; = *( - Л)’ b = [( — ЗЛ)-і b^W = [!(•) e^ b ^’ = (*-l) (fc- ) Din ( ) rezultă că partea netedă a traiectoriei sistemului este xdf) = eAtAkbt cu t> , de unde rezultă că ( ) ( ) ( ) ( ) ж ( +) = Л* &,• ( ) | ’ J ’ • * • - * • С' Г * * • ,/■ w J A j J • Așadar, aplicînd pe comanda i~a a sistemului un „impuls" derivat de k ori, origina este „instantaneu" transferată în starea Лл \ Fie atunci (*-i) m cu м,- e k — Din ( ) rezultă atunci că Uo ku comandă -uk-l - unde u’; eJER mk- Din ( ) rezultă că pentru orice stare xe(Qk există o distribuiională de tipul ( ) care transferă „instantaneu" origina în starea x Considerînd acum conceptul de nul-controlabilitate pentru cazul sistemelor netede, acesta va fi asociat cu particularizarea lui ( ) pentru x = x și = , adică cu ecuația ( ) Cum ел' este nesingulară, rezultă că în cazul sistemelor netede, conceptul de nul-controlabilitate se identifică cu cel de cont? olabilit ale OBSERVAB UITAT EA j Controlabilitatea, discutată în paragraful anterior, evidenția posibilitatea determinării unei comenzi în vederea efectuării unui transfer ш spațiul stărilor Ohservabilitatea, pe care o discutăm mai jos, este o proprietate care pune în evidență posibilitatea determinării unei stări din „prelucrarea" mărimii măsurate y ® începem cu cazul sistemelor discrete Fie sistemul discret în regim liber %(/ + ) — Лх(/) jy(/) t= Cx(t) ■ ( ) Rezultă imediat că pentru t =* O, , , y( ) = Cx y( ) = CA х k — , avem y(k— ) = CA*' x Introducem ieșirea de lungime k prin 'З'(О) -(l) I ■■■*■ - o precum și matricea, numită de k-observabilitate, prin Cu ( ) și ( ), ( ) se rescrie compact în forma У t = Qkxo ( ) д Definiția O stare x este neobservabilă 'n k pași dacă ieșirea liberă pe care o provoacă cu x = x (și «(•) = ) este nulă pentru t = , , , k — ; adică Уц — f> Punem * • ’ f ( ) numit subspațiul neobservabil în k pași Avem atunci în mod eviden • Teorema O stare x este neobservabilă în k pași dacă și numai dacă ( ) • * »Vf A Din ( ) și ( ) rezultă că k k și de asemenea ^Jc+ Din ( ), rezultă că cu Wo - ♦ Ker CA = A Ker C adică subspațiile ^observabile Qlk formează un lanț descrescător Rationînd rența (t iv), lezulta ca exista un întreg astfel încît e л = ( ) sau pe scurt lim ^=^ = , ( ) r • Se poate introduce acum д Definiția care o provoacă = , Vk O stare x g Ж este neobservabilă dacă ieșirea liberă pe cu x = x are primii k pași nuli, pentru orice , adică Pe scurt, dacă xeQtk Vă> Din ( ) din teorema rezultă o Teorema O stare x (X) este neobservabilă dacă si numai dacă ' ; ?? ч - - ‘ r : / Punînd ( ) rezultă că Qt este subspațiul neobservăbil al perechii (С, Л) Notînd п ( ) care poartă și numele de matrice de observabilitatc se vede că rft = Ker(T' » • ГТ • ♦ '} : > * * , \ • j Л л în limbaj matriceal, semnificația șirului ( ) este descrisă de : rang = rang = ( ) •ț Din ( ), ( ) și ( ), obținem că Din ( ) rezultă că adică Qt este un subspațiu A invariant conținut în Ker C Mai mult, cititorul este invitat să probeze că CJL este cel vicii mare subspațrii - invariant conținut De asemenea, considerînd ecuația este soluția suprcmală ( ) ) Definiția Perechea (С, А), se numește observabilă dacă singura stare neobservabilă este a: O, adică PfL erecJi&t {BT, AT) este observabila Demonstrație Demonstrația rezultă c precum și din faptul că r ■, ' ii f • к ЧП este matricea de observabilitate a Utiliziod pnnripiul duabtățu noțra“a duaiă cclei de "" -controla-bilitate a Să trecem acum Fie în acest sens sistemul neted, m regim liber x Ax, x{ ^ == Xq i Â v J r ' • ' la cazul sistemelor netede ( ) rezultînd х , Introducem д Definiția pentru x( = x avem у(т) = O, O Pentru a studia egalitatea Starea este neobservabila la топіепігй t> , ( ) introducem matricea ( ) care, similar matricei W(t) introdusă prin ( ) este pozitiv semidefinită, MO ■ Utilizînd o manieră de demonstrare identică cu cea expusă în teorema avem • Teorema Pentru orice t > (de fapt W> ) avem Ker N(t) = Ker Q ( ) (unde Q este matricea de observabilitate introdusă prin ( )) în continuare echivalentul teoremei este (demonstrația este lăsată ca exercițiu cititorului) ‘д V X • *' T ' , o Teorema Pentru orice t> o stare x este neobseyv abila la momentul t daca și питаг daca xeOt ( ) • t I *• • • Z * v ’J * * * (cu Qt introdus prin ) * r ’* • • , ea continuă să rămînă neobservabilă pentru orice alt moment t> Cu această precizare introducem Definiția O stare xe(p este neobservabilă, daca exista гт t> pent/ru care x este neobservabila la acest moment în consecință, teorema considerată în contextul definiției exprimă un rezultat final formulat prin ® Teorema O stare xe(Xj este neobservabila dacă și numai dacă (i lol) este adevărată Ca și în cazul controlabilității, cititorul poate constata o perfectă identitate între conceptul de observabilitate considerat în cazul sistemelor discrete și cel introdus pentru sistemele netede, dezvoltările făcute începînd cu definiția încolo fiind aceleași și în cazul siștemelor netede DESCOMPUNEREA STRUCTURALĂ Proprietățile de controlabilitate și de observabilitate se reflectă într-o anumită structură a tripletului (А, В, C) pe care o punem în evidență în cele ce urmează Fie și T/f subspațiul controlabil și respectiv neobservabil al sistemului (А, В, C) (asociate cu perechile (Л, B) și respectiv (С, / ) Punem ( ) și introducem subspațiile в ȘÎ ^з> P ™ ( ) respectiv Qt — Ф Жз ' observabil (partea, controlabilă și observabilă a sistemului inițial) IT? «L ih > i Cbj — Л) = Avvin de asemenea (exercițiu pentru cititor) T(X) C'(Xr -Zl) 'B » C»YXJ — ( ) Cu ( ), ( ) Și ( ) obținem • Teorema (de structură) Pentru, orice sistem arbitrar (Л, В, C) există un echivalentul său (â, fi C) cu structura ^AlO) încarc B ; Q «sic observabil și, echivalent wtrare/ieșire cu sistemul ini ial REALIZABILITATE Fie T(X) о pxm matrice rațională strict proprie Introducem a Definiția Un sistem (Л, В, C) este o realizare a lui T(X) T(X) = » C(X Considerăm A {(zl> в' с^> IB ) Pentru aceasta scriem '(X) în forma T/A\ ^fl i' J (X) = -—“ Yo - YA + ( ) /і\ A v v А -I- -I- y A; + X* este c m m m c monic al numi- matricea raționala (X) și unde л(, г — и, ь е Atunci b,„ l ni О» Yo^ пГ ’Y i/ iit este о realizare а lui •іп I * III т і ІП ■-adevăr, punînd ПИІЖ Ѳ(Х) /п х W—* ( ) obținem, după un calcul simplu, (X/ — А) ©(X) ■= ix(X) ( ) sau încă (X) }л(Х) = (XI - А)- В ( , ) Cum С (X) = Ko - -KjX - + înmulțind la stingă ambii mem- brii ai lui ( ) cu C și ținînd cont de ( ), rezultă că T(X) = С(ХІ - Л)" B ( ) Construind matricea de controlabilitate R a perechii (Л, B) din ( ) rezultă imediat că rang R — Im, adică perechea (A, B) este controlabilă Realizarea ( ) poartă numele de realizarea (standard) controlabilă a lui T(X) Similar se arată că ( ) este o realizare observabilă a lui T(X) numită și realizarea (standard) observabilă Dezvoltînd ( ) în jurul lui X = °o obținem unde R este \, i =■ , , "(X) = jX i ~ CgX - , I X R ( ) ales de exemplu astfel încît з(;л(Х)) c UR(Q)* Coeficienții pot fi imediat obținuți prin identificare din ( ) și anume: ( ) = Гі+ - Yj-i^i-i+i - ••• * Yo^ Rx ( ) unde Ra este raza spectrală a matricei A Coeficienții САІ~ В, i — , , k se numesc coeficienții Markov ai sistemului (realizării) (А, В, C) Identificînd ( ) cu ( ) (pentru | X | suficient de mare) obținem Г\ = CAl~lB, i = , , , ( ) adică toate realizările din au aceeași coeficienți Markov, identici cu cei ai lui ’(X) *Ув( ) «{ХЦХІ B t C ) și avînd dimensiunea n • Ă>rang B», rang > rang = n ( ) ceea ce contrazice ( ) Rezultatul remarcabil este dat de • Teorema Oricare doua realizări valențe minimale din ec^~ Demonstrație Fie (А, В, C) si (Л, В, C) două realizări minimale din CfT(-y Evident, ambele au aceeași dimensiune notată n Conform cu teorema ambele realizări sînt controlabile și observabile Din ( , ) avem sau încă ( ) ( ) de unde rezultă ( , ) ( ) Cum este monică, Q„Qn este nesingulară și deci din ( ) obținem ( ) cu Similar, sau mea I ( ) deoarece Rn este epică, este nesingulară, rezultînd din ( ) unde Tot din ( ) se mai poate scrie și - v+; ( ) atunci v și [i sînt indicii de controlabilitate și respectiv observabilitate ai realizării minimale Demonstrația este o consecință banală a lui ( ) LEGI DE COMANDĂ ALOCABILITATE STABILIZABILITATE DETECTAȘI LITATE л •, r i•'ггі ? • - r • -*? ?f» r, T : »•:- • Numim lege de comandă o dependență de tipul ; \ ,r • , t \ ■ I * exprimă modul de schimbare al matricii de reacție la o transformare de coordonate Dacă G în ( ) este nesingulară, legea de comandă ( ) nu modifică proprietățile de controlabilitate ale perechii (A, B) Acest fapt este pus în evidență de o Teorema (i) Dacă perechea (A, B) este controlabilă, atunci perechea (A -f- BF, BG) este controlabilă, oricare ar fi F e și G g nesingulară (ii) Dacă xn definiți prin = + Buk> k — , , n — , xt = b^O ( ) sînt liniar independenți într-adevăr, să presupunem că am construit vectorii liniari independenți xlf x , , xk și că oricare ar fi uk e JES W vectorul xk k obținut prin ( ) acesta este liniar dependent de ^, x , xk, adică ¥ » -• • *%*•!' k > * ~b Buk g =^= т[# x xk cjRw ( ) Atunci făcînd în ( ) = obținem: ЛѵЛ ( ) și deci din ( ) și ( ) rezultă că: ’ Buk g Sk Уик e ( ) Cum xit x f xn sînt liniar independenți, rezultă din ( ) că rang [ , XF , , n Și deci perechea (Лр, b) cu F selectat prin ( ) este contro- labilă Teorema ^ , deși nu furnizează un algoritm de calcul al matricei F, are o semnificație procedurală deosebit de simplă Intr-adevăr, fie un belmB, b^O și fie FI(F) A det [b AFb AnF~l este un polinom de mn feR ”, introducem ■ • - ■ FeJ Lrax” ( ) variabile Asimilîndmatricea F cu un punct din K •, |П'Г) = О} ( ) = ; jgmn pentru care I (F) Ф O și Ultima constatare arată că M este Conform teoremei , există un „punct F £ mn deci II(F) în p p j^mn pentru care densă în M " " vt w'w* "••v” / v îp si în consecință mulțimea punctelor o suprafață algebrică «mținu ta^ K este contro)abilă> este o mulțime Obțineau astfel ca o consecință operantă a teoremei e Teorema Considerînd perechea (A, B) controlabilă, controlabiliMea perechii (AF, b) este generică în raport cu F е^ техл și b= Bg, ^ e Altfel spus, „alegînd aleator" perechea (F, g), perechea ( F, b), rezulta con-tr° Cu acest ultim rezultat, se poate formau și demonstra rezultatul fundamental al acestui paragraf pus în evidența de - • Teorema O pereche (A, B) este alocabilă dacă și numai aaca este controlabila Demonstrație Condiția este suficientă începem cu Be» a“a o pereche (A, ) controlabilă și Л un set simetric de » complexe X„ i = » Arătăm că extstă un/e JEf astfel incit numere n n Z(X) ™ JJ (X — XJ °Ч) i аіЛ ‘ •” / XO* polinomul monic ce are ca rădăcini mulțimea A Atunci condiția ( ) este echivalentă cu ( ) unde , г este polinomul caracteristic al matricei A -ț-bf vtitnh» Întrucît perechea (A, b) este controlabilă, matricea de contnJab^^ asociată R [b Ab este nesingulara în consecința, exista $ astfel îneît ll ( ) sau explicit qr[b Ab Л^Ь] = [ (se vede că este ultima linie din A’ ’) Din ( ) obținem — , -^i^n— , -l’ lb ~ relații cu care se pot scrie în continuare, pentru orice feVS ’ următoarele egalități qT = f qT(A - bf) = V/gKU ( ) f t f t \ t • • > I * | • Utilizînd teorema Cayley-Hamilton și considerînd ( ) adevărată, va rezulta că /(Л + bf) = y^^fA + bf) = , ceea ce face ca din ( ) să obținem f = ( ) Aceasta este chiar f căutat într-adevăr, cu /ales prin ( ) rezultă conform cu ( ) că qT -(A + bf) = ( ) Din ( , ) și din ( ) obținem ГХ(Л - bf) = ( ) ) И qTA*~ aceasta fiind o matrice nesingulară întrucît TR este nesingulară In consecință, cu T nesingulară, obținem din ( ) că z( + bf) = o ' I • ‘ ( ) Cum perechea (Л - bf T, b) este controlabilă, iczultă că + bf este o matrice ciclică cu generatorul b și deci condiția ( ) este echivalentă cu faptul că y(k) este polinomul caracteristic al matricei A - bf Mai mult /, furnizată de ( ), este mică (exercițiu pentru cititor) Trecem acum la Cazul w > I , perechea (AEt,b) va fi con- Alegînd aleator un FQ teoremei , controlabilă, unde Apf~ A + B ?o, b Д Bg fOFL Atunci, aplicînd rezultatul obținut în cazul I, rezultă că alegînd/ = = — (ffA pj vom avea ( ) Dar XF, +bfT = A + BF + BgfT = A + B(F + gf) ( ) și deci F = Fo - gfT ( - ) \ • * t ф • este matricea de reacție căutată , A , > Condiția este necesară într-adevăr, considerînd perechea (A, B}, alocabilă, să presupunem că aceasta nu este controlabilă Atunci, aplicind teorema de descompunere controlabilă, obținem cu ( ) ( ) ж • • ■ * ф ж ® w *■* *• Atunci ținînd cont de ( ) și ( ), obținem din ( ) că C(A + BF) = g(A - BF) = g(z j - B F^ U Ш - S • * > f * * * Г • * e l W * “w / ь • Л л * • Incluziunea ( ) pune în evidență așa-numiții poli ficși față de legile de comandă, constatîndu-se imediat că aceștia se identifică, conform definiției , cu valorile proprii necontrolabile ale perechii (A, B) * Combinînd teorema (de descompunere controlabilă) cu teorema , rezultă că perechea (Alf BJ este alocabilă prin alegerea convenabilă a matricei Flt fapt ce se formulează și prin propoziția: restricția A - BF j & la subspațiul controlabil este liber alocabilă prin alegerea restricției F î n ( ) dacă și mimai dacă; (i) vectorii xit i^ , , , n suit liniar independenți; ' (ii) \k =' dacă și numai dacă xk Xj} ^iii) xk € (A — ImB, i = , , , л Demonstrație Necesitatea primelor două condiții rezultă imediat din faptul că A f BF este o matrice reală cu structură simplă, valorile proprii fiind distincte pin ( ), punînd n ( ) ■л ) eF uns :gile iefi- tfll unde obținem (- — V) x( == — Buit i=l, , ,n ( ) adică, (Л X( x)^j £ Im B, i = , , , n, ceea ce este echivalent cu condiția (ni) Зи/ісіеф condițiilor este procedurală într-adevăr, din condiția (iii) rezultă că există ut eC'“, i= , , ,n astfel îneît ( ) să fie adevărată Din condiția (ii) rezultă că în ( ) vectorii ut se pot alege astfel îneît — Xj dacă, și numai daca ( ) Considerăm acum sistemul -^ ^ •••• ( ) care are soluția ( ) Soluția furnizată de ( ) este reală într-adevăr, din ( ) avem = U M„] П ([* X ■■■ Xn\ П) = — / /re) [a-'j x x^\ — F ( o ) unde s-au utilizat (ii) și ( ), П fiind o matrice nesingulară adecvat aleasă Cu F furnizată de ( , obținem, utilizînd ( ) și ( ), direct ( ) Dacă В este monică, atunci, vectorii щ dacă există ,smt unic determinați și în consecință ( ) furnizează o unică matrice de reacție Din ( ) și din condiția (i) rezultă că soluția generală a ecuației x = (A + BF) x ( ) este de forma x(t) = Cj еВ'ху + + c„ex»'x» eK ” ( ) unde X( = Xy dacă și numai dacă constantele arbitrare с,- eC satisfac cf = ,- Atunci ieșirea „liberă" a sistemului va fi y(l) = Сі^Уі + + cn ex«*y» еЖ” ( ) Э Уі C Vp — , , , Din ( ) și din ( ) rezultă că avem ( ) ( ) si a vectorilor yit i — > ••••» * liber ( ): dîndu-se setul simetric A de numere n complexe distincte și vectorii \ ' -r • T? Г~И Ш \z ? — - - Z * / - / J л к г • Cu aceste precizări se poate formula și problema alocării valorilor proprii — , n, cu ajutorul cărora să se scrie orice răspuns » v v >W să se determine F eK wxn și vectorii liniar independenți xn> astfel încît ( ) și ( ) sa/te simian satisfăcute Obținem astfel • Teorema , Dîndu-se setul simetric Л de n numere complexe distincte și vectorii Уі, y%, ••••; Уп problema formwată mai sus are o soluție dacă și numai dacă: (i) există vectorii liniar independenți x±, x , xn^Op precum și vectorii ult u f , um e O\ astfel încît ( ) să fie îndeplinită; Xj, = Șî Уі = У г Xj daca și numai dacă xt Pentru ca ( ) să aibă o soluție, oricare ar fi vectorii yt, este necesar și suficient ca rang = ranv O == p rang [A — B], i = , , ,n ( ) Presupunînd perechea (A, B) controlabilă, obținem din ( э), cu ( ), adică trebuie să avem m p, iar X/să fie, pentru i = , , ( ) r | r n, diferite de „zerourile tripletului (A, B, C)“ (v § ) Față de cele de mai sus, se va proceda astfel' — se aleg numerele complexe distincte cu , bm sînt primii m (din cei n) vectori liniar independenți selectați ; — dacă Ak~l bt figurează ca o coloană selectată în blocul Ak B, dar Akbt nu figurează ca o coloană selectată în blocul AkB, atunci AJ ^ nu va figura în continuare în nici un bloc AjB cu j>k Cu aceste constatări, cei n vectori liniar independenți selectați se pot scrie într-o ordine lexicografică de forma •••• Ani~lb, b Ab ••• bm A*»»” bm] ( ) unde întregii n,( > , i = , , , m cu -|- n > + nm = n se numesc indicii structurali ai perechii (A, B) Se constată direct că v А щах щ este chiar indicele de controlabilitate al perechii (A, B) Introducem de asemenea și numerele ст ; = n, -p n ф nit i = , , , m ( ) unde = nL, iar == n A I Inversăm matricea Se și selectăm din S liniile cu indicii C rezultă imediat din ( ) si ( ) că * • j Aq = A “ este o pereche controlabila, iar /л, (X) = Xn Se pune în evidență astfel o procedură efectivă de calcul a reacție introduse prin teorema Perechea ( ) obținută mai sus este o formă canonică de tipul ( ) Avantajul acestei ultime forme canonice, pentru m = , este vizualizarea efectivă a matricei de reacție, împreună cu unicitatea ei, care alocă un spectru arbitrar într-adevăr, dacă yw(X) este polinomul monic ce are ca rădăcini setul simetric arbitrar A, fiind exprimat prin ( ), atunci alegînd A ” [ « al •” a- -l ! ( ) obținem direct din ( ) că ( ) și deci rs(A + bfT) ~ Л m matricea Intioducorn, utilizînd indicii structurali îl polmomulă eM”X”[x] definita prin ( ) Obținem atunci egalitatea remarcabilă (XI — Лс) SC(X) = BcP^PnM ( ) unde P M - ЛИВ(Х) ( ) demonstrația făcîndu-se prin calcul direct Fie acum un sistem controlabil arbitrar ( , В, C) și fie Tc transformarea ( ) care aduce perechea (Л, B) la forma canonică (Ar, B, ) din ( ) ( ) Punem ' СсАС / ( ) Atunci rezultă din ( ) și ( ) că avem T(X) A C(Xl - ri)^B « Co (XI - Л,)-* Bc = = CcSc(Ă) Z^x(X)Bw = PC(X) P;x(X) ( ) unde ț Л (Х)'А СД(Х), P (X) A B^PM(X) ( ) ■ Dacă notăm cu gradul coloanei i dintr-o matrice polinomială, se constată imediat din ( ), ( ) și ( ) că avem /X) = / АЖА Atunci Uzări Factorizarea obținută ( ) permite scrierea imediată a unei canonice de tip bloc companion corespunzătoare unei matrice raționale stiict proprii T(X), utilizînd următoarea Procedură de realizare — Se factorizează T(X) în forma ( ), ( ), reținîndu-se gradele = I j ~ , ,, w — Se construiesc matricele ■radele m ( ) j unde cu ( ) ( ) precum și matricea C astfel încît Rc(k) = CSJX) Realizarea obținută (А, В, C) este sigur controlabilă Obținerea rea-та П та- е se face aplicind eventual teorema de descompunere observabilă (de fapt de descompunere controlabilă aplicată dualului) Revenind din nou la forma canonică ( ), ( ), să considerăm le^ea de comandă ( ) un ce „A w = Fcx + G, ' ' fi ’ ** ’ A AL E>- Д C T>- c k:rc~^ iu Se constată imediat prin calcul direct că -Г В,Вc d^ag [- CJj, B^GC = diag m CU ( ) definind forma canonică bloc companion d& sinteză Să explicăm semnificația structurii ( ), ( ) Pentru aceasta considerăm mulțimea tripletelor î \ P * • ’• i \ G = {(Г, F, G) | 'gI*»x", Gel-«”*xw nesingulare, FeBL^”} ( ) Voni spune că două perechi ( , / ), ( , B) sînt if echivalente, nofîndu-se (Л, В) ~ (A, B), dacă există (’/', F, G) e ( { “ | j t x «Я i • l} L i ‘ Dacă V este o matrice bază pentru (V, atunci ( ) este echivalentă cu existența matricelor J si H astfel îneît să avem AV-VJ^BH ( ) Notăm cu - • ■ ■ * * ™ - * Z • ■ • С • Fi •••# a • Э(А,В) ^{CU с Ж| Аф с: Ф + Im B} ( ) O proprietate remarcabilă a subspațiilor ( , B) invariante este pusă în evidență de • Teorema , (B e (A, B) daca și numai dacă exisiă F astfel îneît AjXV c astfel încît Im E а ф c Ker D ( ) Demonstrație Condiția ( , ) este suficientă Intr-adevăr, alegînd unF e ( ) va fi adevărată și în consecință vom avea ПФр(/)£ == V/ ( )' ceea ce echivalează cu ( ) Condiția ( ) este necesară într-adevăr, dacă ( ) este adevărată, atunci ( ) este adevărată» Punînd, în consecința СО CfcF, unde este subspațiul controlabil al perechii (JF, E), rezultă direct din ( ), că dubla incluziune ( ) este îndeplinită Fie acum Ol c (£• un Subspațiu arbitrar și fie introdusă familia ЗГ(А, В; tft) А [ф а un lanț crescător (descrescător) de tub-spații din Rtn, adică CUQcz CU^ CU a , ( ) (Фо o ф, => Фг о ) ( ') • ' • * * Atunci dacă există о aplicație izotonă F, definită pe acest lanț, adică = ( ) există un întreg astfel încît să avem C* • t »’ I ‘ И ’ У ) ci) cuv — cu^+\ — сиѵ — •••) л l \ Ti ( ) ( ') Demonstrație Deoarece t este de âdmenswne finită egală cu n, există sigur un întreg pentru care avem - ; “ CUq , este un lanț descrescător și lim фк = ф* = sup g(A, B; Ql) ( ) cu ф*ед(А, B; Ql) , Г** ' i * i * J * » • / Demonstrație Să arătăm prin inducție că ( ) definește un lanț descrescător, într-adevăr Фк = Ql cc (£> = фй ■ Dacă acum Фк g Фк^, obținem din ( ) că ^t+ = ФІ П Л- ^! + Im В) g Ql П А~\фк + Im В) = к А * ț * f * } ? ) • * Г'к ’ ' * * f ) ' ' - * / > și deci Фк[ Aplicînd teorema rezultă că există un întreg = Фо Cum ф e DțA, В; Ql), evident avem Ф cz Ql П А г(ф + Im B) ( ) Dacă acum Ф g фк atunci, utilizînd ( ), obținem Фш = П А-\фк + Im В) g П А~г(ф + Im B) g ф și deci ф g ф/с V k > , adică ф c w • ț Din inspectarea, pentru , a lui ( ) rezultă o Procedură (elementară) de construcție a lui Ф* Se alege o matrice bază V\ pentru M, se face К și se determină matricea C ca soluție epică maximală a ecuației CV = - ( ) adică Ker C = Im V == Ql ( ) Se calculează soluția epică maximală W a ecuației W[V • B] = adică Se calculează soluția monică maximală F' a ecuației л WA ( ) adică = M П А~\ф + Im B) (CV A Im И') ( ) unde s-au utilizat ( ) și ( ) Relația ( ) explici tează în fond recurența ( ) (cu /(>!) Este rang V — rang V'? Dacă „da", atunci Ф* = ImF' și STOP, Dacă „nu", atunci se face V F' și se reia procedura de la pasul °- Noțiunea de (A, B) invarianță poate fi dualizată couducînd la noțiunea de (C, A) invariantă In acest sens, luind complementul ortogonal în ambii membri ai lui ( ) obținem Ker BT ( «) sau încă A Ker Вт) с ф-L ( ) de unde făcînd schimbările duale AT -> A, BT -> C și ф-L -* obținem din ( ) A{ () A Ker С) g ( ) Introducem atunci A Definiția Un subspațiu cz Qj se numește (C, Л) invariant dacă satisface ( ) >, a De asemenea, dualizînd ( ) obținem {А^фі- ( ) sau (AT'+ FTBT) cz ( ) de unde cu substituția duala A T A, BT—>C, (V^—^TO și FT —> К obținem (A + KC) cz У ( ) Așadar, obținem duala teoremei care este e Teorema Un subspațiu cz Ж este (C, A), invariant, dacă și nuynai dacă există astfel încît ( ) să fie adevărată ■ ■ Similar cu ( ) și ( ) introducem g(C, А) А{ту c |И(^АКегС) g } ( ) și respectiv | ■ / ЭЦТО} ={K I (A + К C) g } ( ) Я^Ві ~ ‘ \ **}['* / î? ! / - ’( І ’ I ici £; pentru orice e Э(С, A) Dacă acum & cz Ж este un subspațiu arbitrar, dualul familiei ( ) ' este • ~ ’ S (j; С, A) = { c (£ | eg(C, ) și J c } ( ) Vom spune că * este un (subspațiu) familiei c ( * V e C, A) (jj) dacă => ТО' V e S?(J; C, A) atunci ' g * și se notează * =inf e ( , O = { } ( ) Atunci k, /с> este un lanț crescător și lim k = **= inf g(S; C, A) ( ) cu *еЭ(£-,С, A) în fine, duala teoremei este • Teorema Fie un triplet (A, B, C) și fie pentru un X e X£ x- arbitra? Tf subspațiul controlabil al perechii (A + КС, В) Introducem familia de subspații , )* etp ( , , adică * este cel mai mic subspațiu control abil fața de legile de injecție definite prinKeW^ Așadar, recapitulînd, dacă (А, В, C) este un sistem oarecare, subspațiile СО* = sup S (A, В; Ker C), * = inf £ (Im В; C, A) introduse mai sus semnifică cel mai mare subspațiu neobservabil al perechilor (C, A + BF) si respectiv, cel mai mic subspațiu controlabil al perechilor (A - КС, Bf ’ SUBSPAȚII DE CONTROLABILITATE • l * fi •• w ■ I t \ • k\ f • • • Dîndu-se o pereche (A, B), subspațiul controlabil al acesteia este, așa cum s-a văzut ч , V о unde evident n n -f- i У) Ahc^ = A 'S cu i , , h k = i Atunci, cu ( ), ( ) se rescrie compact — )) (p /I -P BF I Im BG > ( ) , b' c M , VFer/(«) ( ) Demonstrație r ■ Să arătăm mai întîi că egalitatea ( ) este adevărată pentru matricea F, cu care s-a definit prin ( ) subspațiul , unde F este cel introdus prin ( ) Cum (Л BF) (A + BF)k (Cp fi Im В) a cp yk Cum Im BG e cp П Im В rezultă că (A + BF)* TmBG c (A + BF)“ (cp fi Im В) ѴЛ> Din ( ) și ( ) rezultă că cp = q£ Arătăm acum, că egalitatea ( ) este adevărată pentru orice Ге Fie în consecință un F e(/ (cp) arbitrar ales și fie (P ( -|-+ ^F| CP f) Im JB) Cum (A -f-BF) A (A + BF)>‘ (CP П Im В), Л > ( ) Atunci jt+ = (А -I- BF) tSk, k> ( ) Și == ( ) ! Evident Л, == Ф, Л îm В с (Л- Presupunem că J* o și deci g^ ( - ) • r ■ • ' , 'ț , / ’' ’ д'" '■ ■ Din ( ) și ( ) rezultă că = și teorema este complet demon-strată li Notăm cu @,(A, B) totalitatea subspațiilor de controlabilitate ale pere- O descriere sintetică a subspațiilor de controlabilitate este dată dc • Teorema - ' A, B) dacă și numai dacă (i) b ț definit prin е@\А, B) Atunci conform teoremei , (Л se explicitează prin ( ), adică , ■ \ ■ ■ ■ 'Л (A + BFf ( П Im B) ( ) i> pentru un F e ( ) ( ) i ar ( ) Să arătăm că lim * "* ( ) = , ( ) este adevărată cu * definit prin ( ) Pentru Z Presupunem că * * Atunci, deoarece JF g (^( ( ) Atunci, utilizînd legea de distributivi tate permisă de ( ), avem din ( ) g'+ = П Im В -I- (Л + BF) (Rf^CR, A (Im В - (Л + BF) (R'k) = = A (Im В - *) = CR A (Im В - A(Rk) = rg*+ ( ) unde s-a utilizat și ipoteza de inducție Astfel ( ) și nsCesitațea cornii, lor (i) și (ii) este demonstrată Reciproc Presupunem condițiile (i) și (ii) îndeplinite Fie Fe( ( ) Atunci, parcurgînd ( ) de la sfârșit către început (curgț = = *) obținem CRM = A Im В - (Л - BF)fRk, adică ( ) este line definit {independenl de alegereft- Ini F Utilizînd în mod i dentic aceeași manieră de demonstrație ca în teorema , rezultă că ' ‘ * = lini * ' * î ( ) cu • »+ == ( ri‘ A (Im В -Л *), A> , == ( ) I v ( și deci ( ) este independentă de W®*) Să arătăm cădi* definit prin ( ), satisface ( ) Pentru aceasta fie Гй e , F' e Cf-{dH) ( ) Dar fg' a ( ) V A \ \ І • Kt* ’ • ' ' * ’ ' ' •* • ' și în consecință se poate scrie (O* = (îi! ф СО' ■ ( ) Construim acum un Fo astfel îneît Fo| di' = F' și Fo| d)’ = F ( ) unde F si F' au fost introduși prin ( ) și respectiv ( ); Construcția unui asemenea Fo este posibilă datorita lui ( ), m sensul ca di Л ' = (A + BF ( ) : ’ ? l V Д I • ч C- iar din ( ) avem СЦ! = ( ) Comparînd ( ) cu ( ) și ținînd cont de ( ) rezultă că dH c * V dH e @(A, В; dl) ( ) iar pe de altă parte ( ) explicitează faptul că di* e g(A, В; City ( ^) Din ( y și ( ) rezultă că di*, definit prin ( ), satisface ( ), teorema fiind complet demonstrată > ■' ' , Oo = { } ( ) și respectiv ^+ = ^*n(ImB + zl A), Л> , '^ = { } ( ) Dc asemenea (V* = Ker C n ' + Im B) ( Al I) r * Г ' Ш к * • (ca soluție supremală) Arătăm prin inducție Ca ( ) este adevărata în sensul că % = ( ’ ) > v * J A • J f / * f i • z ' *• \ a • ? J într-adevăr, pentru h = , ( ) este adevărată Presupunem că ( ) este adevărată pentru un k și arătăm că este adevărată pentru h - l Avem atunci din ( ), ( ) și ( ) = Ф * П (Im В + А k) = П (îm В -A(^* П Ш) = fîl* (ІтВ + (КегСп/ ->*-|-ІтВ)П О ? = СО* П (Im В - {СО* - Im В) п А{ к А Ker ) = С()* п {Ф* + Іт В) А (Іт В -|- А\ к А Кег С)) = GO* к+ Trccîrid acum la limită în ( ) obținem ( ) Aspectul operant remarcabil legat de subspațlile dc controlabilitate (sie furnizat de • Teorema B) clacă și numai dacă pentru orice set simetric A de p numere complexe, unde- p л dini , există I'e t tm " astfel incit (Л | BB) С Și ст(Л I BFI CU) A ( ) adică e (d, B) dacă și numai dacă este liber cloca’ ’il DewtoHslvafie Vom efectua clemonstrftția în doua etape Fie ГЙ $(Л П) și fie descompunerile Im U « ImBfl бйфіт# ( Л ) împreună cu ) =» Ф S ' astfel încît Im В c *Asup A sup к * » • % * â J • \ • \ J l * J * > = {F + GF'T\ VFn, F , F , F , F și F L = } ( ) pentru uir F e Г ((У *) arbitrar ales- Matricea Az este independentă de alegerea matricei F e Cp( V' ♦ / r : P r • f * ♦Л • ' ; ’ Л г І V ) Л / b inversabil la dreapta dacă Rang T(X) = p ^ / Un sistem simultan inversabil la stingă și la dreapta se numește inversabil ♦ ’ ' -?î ' > Д / ; ’ A ’l Z - ' ‘ H Avem atunci Tf » • * A MV Я a — • '* • Teorema Un sistem (А, В, C) este a inversabil la stingă dacă și numai dacă * В monică și * = ( - ) b inversabil la dreapta dacă șinumai dacă ' C epică și Й* + == Ж ( - ) У f ' t Ar • \ І" \ Demonstrație Din ( ) rezultă Rang T(X) = m rang В = m și (V * Q Im В — { } ■ ‘ i ‘ i' , f z /,*, ; ? V к >; * /* ) ț • ■ В monică si * = • J b rezultă prin dualizarea lui ( ) • % = (ф* n + N, P) este un cstimalor de F (sau al dependenței Fx^), I sau M) daca lim (^(/) t >-? r/i unde x(l) este starea sistemului ( J), iar w(/) este ieșirea (măsurabilă) a sistemului ( ) ’ — Dacă în plus (ii) ICJ- sau, ГЛ(О) ( ) vslimalot'ul se numește stabil, iar dacă, (iii) °(J) — cu Ac C, set simetric arbitrar de numere, ‘ estimatovul se numește (liber) alocabil Schematic, cele relatate in definiția de mai sus sînt ilustrate in figura , Ca terminologie, daca Z* = dn, estimatovul se numește de slave, iar daca l =■ m se numește de lege de comandă', în particular, pentru l — se numește tsiinurtor de funcțională (liniară, ==/ ), z i X’ Re elRtvx*, o matrice necunoscută și introducem eroarea (dc estimație) , ( ) unde Ф (/) este matricea fundamentală a lui J Este imediat vizibil că ( ) este echivalentă cu (J, К, // Л , N, r) lim K >j(t) = ( ) Dăcă estimatorul este stabil sau alocabil, atunci condrpa ( ) este totdeauna satisfăcută Cu detaliile de mai sus și ținînd cont de ( ) și ( ), definiția devine mai specializată prin д Definiția Sistemul (v, J, К, H, V) (M și P au fost suprimate datorită lui ( ), ( )) se numește un eslimator de F pentru sistemul (А, В, C) dacă există V e JRt x" astfel încît ecuațiile ( ) să fie satisfăcute (fâxă ( A )) Dacă în ( ) se impune N = ( ) estimatorul se numește de tipul ; in caz contrar, dacă N figurează ca necunoscută, estimatorul se numește de tipul Dacă în plus ( ) șițsau ( ) sînt îndeplinite, estimatorul (de tipul sau ) se numește stabil, respectiv (liber) alocabil Investigarea ecuațiilor ( ), amendate cu ( ) sau ( ) se poate face direct Se preferă însă o „geometrizare" a acestor condiții, geometrizare care se face cu ajutorul noțiunii de acoperire dinamică Avem astfel • Z » * ‘ * * ' ' * i A Definiția Fie o pereche (A, B) și E eJE^nxr, D n subspațiu cz se numește acoperire (dinamică) de tipul (tipul ) a lui Im E în raport cu perechea (A, B), dacă este (A, B) invariant și acoperă Im E {(V + Im В acoperă ImE) adică A ' * ■ w * •• * /* Jr w ■ Deci, ( ), ( ) definesc o acoperire de tipul , iar ( ), ( ) definesc o acoperire de tipul , a lui Im E Dacă V este o matrice bază a lui CV atunci, conform cu ( ), ( ) și ( ), există matricele J, H, К și G astfel încît sau Dualizînd acum ecuațiile ( ), pentru N = , obținem sau = VTJr + / г утКт , CrNr pentru N liber Punînd în ( ), ( ) Ar / , Cy -* Ii și V -» И У' ~> , jjt H, FT -* E și NT astfel îneît, pe lingă ecuațiile ( ), ( , ) sau ( ) să avem д(Лр ч ч Ч л • * * а X іѣ Л ' І ' • \ ’ ’S i \ ’ f \ V • I * Л * Л \ • J - ' • • * ч > » Comparînd, în contextul dualității ( ) și ( ) cu ( ) și ( ) respectiv, rezultă că acoperirea dinamica stabilizabilă este noțwnea duald estimatorului stabil, iar acoperirea dinamică alocabilă este noțiunea duală estimatorului alocabil (v tabelul , ), * i — -j, ACOPERIRE DINAMICĂ Dualitate ESTIMATOR Tipul AV; = VJ + BH Tipul VJ E l' II Stabili zabi Alocabibl cu corespondențele Stabil Aloc abil Cu precizarea dualității intre noțiunile de acoperire (dinamică) și cea de estimator trecem în continuare la studiul acoperirilor dinamice Mai întîi să evidențiem faptul că determinarea unei acoperiri dinamice de tipul se reduce la determinarea unei acoperiri de tipul într-adevăr, considerînd matricea В monică, presupunem că (Л, B) este f°rnra ( ) ’ în consecință, ecuațiile ( ), ( ) perechea n t zx ^ se scriu Г O ( > Г O ■ш {[ E J [ v sau încă, din primele n — m linii ale lui ( ) și ( ) obținem iar, din ultimele m linii ale lui ( ) și ( ) rezultă • j" ' ■ • ••• • -J • ț- * ■ • ’ л Ѵ; + л^ = vj + н Din ( ) se vede că V\ este o acoperire de tipul a lui Im E, perechea (Alf A ), afirmația făcută fiind demonstrată Un caz remarcabil este cel în care Acoperirea de tipul se scrie atunci din ( ), ( ) în îoima V Pentru ca V să fie monică alegem și deci ( ) se reduce la este stabilizabilă, acoperirea Din ( ) rezultă direct că dacă (/ , B) , * •• • * F Ct c '■*> Fie v o baza a lui ф si A^mci ( ) este adevărată dacă și numai dacă perechea (/ К i și ІтККд = un subspațiu oarecare Atunci, daca exista F e CțF( , (Vo = Fj ( ) Atunci evident din ( ) rezultă că (Vo , de unde cu ( ) rezultă că ( ) este adevărată Atunci, din ( ) și ( ) rezultăcACO = CO' Ș deci F - СО П Im В este un generator al lui ( Cititorul este invitat să arate că, bazîndu-ne pe teorema avem e @( , B) e p(A, B) și f) Im В este un generator al lui ' Introducem acum A Definiția acoperire * Rezultă imediat că dacă (Z) este o acoperire de tipul de dimensiune minimă a lui Im E, atunci Im E este un generator al acesteia, CV fiind în consecință o acoperire minimală Reciproca nu este în general adevărată, adică o acoperire minimală a lui Im E nu este și o acoperire de dimensiune minimă Alai exact, acoperirile de dimensiune minima se găsesc printre acoperirile minimale, avînd toate aceeași dimensiune și anume cea mai mică posibilă în consecință, fie ) [E BJ este monică, reamintindu-se că E are r coloane, iar В are m coloane w Fie T transformarea care aduce perechea (Л, [E B]) la forma canonică bloc companion (v § ) Atunci, evident coloanele matricei S " formează noua bază a lui în care perechea (Л, [E BJ) are forma cano acă bloc companion Cum [E B] este rnonică, cu r -|- m coloane, indicii structurali ai perechii (Л, ГЕ BJ) vor fi w Punem ( ) si scriem forma canonică bloc companion a perechii (Л E В \ conform partiției ( , ) unde ЛиеИ "‘хиЛ i,j= , , £ eJE tn‘xr, Bs e Ып*хи iar perechile (Ju \), (Л , B ) sînt fiecare în forma canonică bloc companion Fie acum matricele Ar+m și [E B}r rTn formate cu liniile ce au indicii * •’ c, = V «*> ' Q^i^n}—\, j = i, , ,r, e, fiind coloanele matricei coloanele matricei unitate In,, iar Лиё, = , n} fiind cei mai mici întregi care posedă aceste proprietăți în consecință, (v ( )) (Vx este un subspațvu de controlabilitate de dimensiune minimă al perechii țA, [E BJ) care conține Im В și deci este adevărată > Фк ’ b * % » * | b r i ь f ' * Jx C • • • * • Teorema esfe o acoperire de tipul și de dimensiune minimă a lui Im E în raport cu perechea (A, B), Dimensiunea minimă este = ^ni (cu ni indicii structurali ai perechii (Л, [E BJ), ipotezele a) și b) i=l formulate'la început fiind adevărate) ® Să explicităm acum, în maniera ( ), ( ), acoperirea dinamică de dimensiune minimă (If, introdusă mai sus Avem evident din ( ) o , r e, fiind coloanele lui E, iar F eCf-fty Or vectorii Ață cu ■ • • " formează о / J ' i x I rezultînd estimatorul («j, J, К, H, KB), care este un estimator de F de tipul și de dimensiune minimă B ALGORITMUL DE OBȚINERE A UNUI ESTIMATOR DE F DE TIPUL Șl DE DIMENSIUNE MINIMĂ Se fac aceleași operații de dualizare ca la pasul al algoritmului A, după care se fac schimbările E В și В ч— E Este identic cu pasul al algoritmului A Se determină cu ajutorul relațiilor ( , ) matricele V, J, și C Z Se dualizează rezultatul do la pasul prin V vr, ,/ ', J llT, К o V* е г О d О о О О о с о о о e -l I -Л l о o e e o o o o o o о o o o o o Indicii structurali perirîi de tipul este лх «о «s s г ч l și deci dimensiunea Взи& canonici este тіпітЛ c e e o o О ferma canonică bloc companion fiind г i I ш w w жж м жж i жжжжжжжЖ ^жжмжжжжжжч - [£Дк- * - - ~ *• — ч! чЧ ѴЧЧ%\К к Ш ж иажх'сж L o і > а Rezulta de asemenea Г-з І - ~ ч ч ч - ! п о — » і *• •*■»•«■»****>*** •>»♦»>•»• % W**'ikW*> **’»’* X* W’ ч* *%ХЧ*Ч - ) ■ф - • « ** v i - - r - - , R,|r | w =s — - „ i ; de vvdc Apficînd formulele ( ) rezultfi Oi - O — J , Dualizînd rezultatele de la pasul obținem I f ■« * și respectiv A = {— -£ ^/ } ceea ce arată, că cei doi estimatori obținuți sînt ambii instabili, priviți atît ca sisteme Betede cît și ca sisteme discrete, ACOPERIRI DINAMICE ALOGATE în subparagraful , prin introducerea noțiunii de acoperire de dimen-siune minimă (de tipul sau de tipul ), am rezolvat în fond următoarea problemă: să se determine cel mai mic întreg v, cu proprietatea ca există o matrice V eH"xv împreună cu matricele J, К, H, (G) astfel îneît ecuațiile ( AS), ( ) (( ), ( )) să fie satisfăcute Soluționarea acestei probleme rezolvă prin dualitate problema estima torilor de F de dimensiune minimă, dimensiunea acestor acoperiri fiind v = pentru tipul și v = pentru tipul , unde ?! este suma primilor r indici structurali ai perechii (A, [E B]), iar n! suma ultimilor r indici structurali ai perechii (Л, [В Ef) în cele ce urmează, alături de minimalitate sîntein interesați însă în obținerea de estimatori liber aloc abili, în acest sens, formulăm următoarea problemă: I, Să se determine cel mar mic întreg v pentru care oricare ar fi mulțimea simetrică A de v numere complexe există o matrice V eRMXV împreună cu matricele J,K,Ht(G), astfel îneît ecuațiile ( ), ( ) (( ), ( )) să fie satisfăcute și în plus să fie îndeplinită condiția de alocare a(J) = A ( } Cel mai mic întreg cu proprietățile de mai sus se numește indicile de acoperire al lui Im E de tipul (tipul ) și se notează cu va Evident, ( ) & = Vk ( ) și unde vom presupune perechea (J, k) în forma canonică controlabilă, adică ( ) unde v este pentru moment nedeterminat Punem în evidență coloanele matei-celor V și H scriind ■ P = v„], H = [A ( h, Лѵ] ( ) înlocuind ( ) în ( ) și ( ) obținem ecuațiile Avt ~ pv + Bkt Л Ѵ = Bk$ ( ) Av* — vv y -p vvv + și respectiv e ( ) Eliminînd и «ii’v din ( ) și ( ) obținem imediat ( ) Dar arbitraritatea alocării spectrului matricei J revine la arbitraritatea coeficienților (, lo în consecință, ecuațiile ( ), ( ) sînt satisfăcute și deci întregul minimal v determinat cu ajutorul incluziunii ( ) este chiar indicele de acoperire vu (de tipul ) В Similar avem și " • Teorema Indicele de acoperire v' de lipul al lui Im в este cel mai mic întreg pentru care Imftf Ae A^e] cz Im [В А В A^B] ( ) Demonstrație Ecuațiile ( ), ( ) devin acum AV = VJ + BH ( ) z КЪ Л \ Ж ша w ш e = Vk + Bg ( ) Utilizînd ( ), ecuațiile ( ), ( ) conduc la ecuațiile ( ) și la ecuația e = vv + Bg ( ) care se substituie ecuației ( ) Eliminînd ѵг vv din ecuațiile ( ) și ( ) obținem ecuația [В А В A'B] ^ ~~ g a • ♦ ♦ [г Ae ЛМ * a Л» — vg — Ѳѵ - s Din inspectarea ecuației ( ) rezultă că dacă v este cel mai mic întreg care satisface incluziunea ( ) atunci acesta este cel mai mic întreg pentru ■care ecuația ( ) are o soluție oricare ar fi numerele , Odată v determinat, fixăm un set simetric arbitrar de v numere complexe X» v cu care determinăm (v ( )) numerele , pe care le înlocuim în membrul drept al ecuației ( ) Soluția acesteia furnizează, cum direct se vede, matricele g și H Utilizînd afum ( ), cu inițializarea vv — =, • / t • * V * • ;• j /T) -adică a acelor estimatori de F = fT Avem astfel Algoritmul de construcție a ținui estimator alocat de tipul I (tipul ) pen-iru j » Se fac operațiile de dualizare Л A Tt В CT, C BT^ Sq determină cel mai mic întreg v pentru care incluziunea ( ) (( )) este îndeplinită Dacă perechea (A, ) este controlabilă, efectuarea pasului l este totdeauna posibilă Din specificarea setului Л se determină prin identificare coeficienții Ѳ,-, (v ( )) și sc rezolvă ecuația ( ) (( )), unde v a fost determinat la pasul In cazul cîhd nu se urmărește efectiv o minimizare dimensională, se poate alege în ( ) v = v (iar în ( ), v v ■ I) Se obține matricea H (matricele g și H) ” Se aplică relațiile ( ) (pentru tipul de acoperire inițializate cu = e — Bg) și se determină matricea V Se construiesc matricele J și k furnizate de ( ) Se operează dualizările TZ , ' ' M = VB Cititorul este invitat să constate că dacă perechea (C, A) este observabilei cu indicele de observabilitate v , se poale construi întotdeauna un estimator alocat de funcțională de tipul (tipul ) avîrid dimensiunea v (v„ — ) Exemplul Fie sistemul Г descris prin o o o o o o o Г ()■ [o o pieccm și a locabili Se cerc să se construiască doi tipul Aplicăm algoritmul de mai sus Ѳ Efectuăm operația de dualizare și obținem ? l • b primul de tipul , iar ос de aî doilea de Determinăm indicii de control abilitate ■o " , f — acoperire și respectiv vR Pentru aceasta construim matricea de o °-| - - o o O prccuna și matricea - • ,[e Ae Лае] =*- ■ - Se constată imediat că Im [e Ae A e] a Im [В ЛВ] șt deci precum și Im [e Ae] a Im [В AB] ș deci Așadar, indicele de acoperire de tipul este , iar cel de tipul este egal cu L în cazul acoperirii de tipul , pentru o specificație p a ecuația ( ) se scrie rezultînd - - Г - - - h\ ^ • l hl • ■ o * ~ Ѳ - +- " H = [* л ] = “ - i V - г Og - - - Л g ° g l g* „ cu soluția -Ѳ Utilizînd ( ) obținem pentru acoperirea de tipul = •— / іЛд L o J deci V « în cazul acoperirii de tipul t deoarece Matricele J fi U sînt — pentru tipul de acoperire — pentru tipul de acoperire » Dualizînd rezultatele obținem — pentru estimatorul de tipul Г “ V = , Og ^ o o o г —ѳ» H — pentru estimatorul de tipul V = [ ], : У - ч « o o o o , o ' « avem numai OJ О, К = , H = [ - O OJV ], M = VB Abordăm acum Cazul r > , LTn prim rezultat este obținut prin aplicarea procedurii din cazul r = К și este furnizat de • Teorema Dacă v (v') este primul întreg pentru care Im [E AE AVE] c Im[B AB A^B] ( } (Im [E AE А'Е] cz Im [В А В А* В]), ( } atunci există o acoperire alocata a lui Im E de tipul (tipul ) avînd dimensiunea Uf) Л, I Atunci punîiul I Гг] • • •• ' • • l t Л ■ ■ gf •ceea ce conduce la considerarea, alături de matricele ( ) și a matricei ( ) ( HO) și ( ) satisfăcînd ecuațiile ( ), ( ) în fine, să observăm că •dhnen iunea acoperirii este inferioară întregului rv (r/) deoarece matricea V 'introdusă prin ( ) nu este în mod obligator monică Să precizăm că dacă v este cel mai mic întreg pentru care incluziunea ( ) , perechea (/, Z > ( ) Ui ; O л\ ii) U'» ^o) es^e controlabilă atunci orice set simetric A de v numere complexe este alocabil în (V, adică există FeCfdfty astfel încît | а(Л = Л ( ) Demonstrație „ j j Fie A un set simetric arbitrar de v numere complexe Conform condiției (ii) există L astfel încît * ; f \ л і% ДоД) ^= Дг ' i ( ) Din ( ) și ( ;' ) obținerii atunci ІіГЛПЯі СЮЭ- j kr X- »ч z AV-V(J + K L) - FA L + BH = V(J + ед - BG L + BH ( ) ■ ' == V(J + K E) + BfH ^ GoL) ■ A > Cum V este monică, ( ) și ( ) sînt echivalente cu ( ) •Л l U -• Г / Л A Cx *i л X V* în fond, demonstrația df măi sus are la bază faptul с&ф — Im Л, В) cu V monică și satisfăcînd: condițiile (i), (ii) clin teoremă, este un subspațiu de controlabilitate, adică (V е б(Л, В) și reciproc, aceasta decurgmd direct din teorema De asemenea, cititorul poate constata imediat că egalitatea ( ) cu J + K L alocată, nu reclamă expres cerința ca Г să fie monică axeasta fiind mai mult o cerință de natură geometrică, care să permită intet-pretarea de mai sus »v Cu constatarea făcută în teorema obținem următorul rezultat remarcabil furnizat de ' • : • Teorema Fie V e o matrice (nu în mod obligatoriu monică) pentru care există J elESyx^ și H с ЯІ,ЯХ'Ѵ astfel încît ( ) este verificată S CW ItlfâtT'ir'/CirC > ***/ ^ * йСчАо atunci Daca există un întreg îneît q^O împreună / - > > -t'r ■,' ' • ' - • J • f ' ; Г' '■ • ri-'i JO J : L T Teorema pune în evidență un mod de extensie a lui V (cu un număr de coloane egal cu jxnr de coloane ale lui Ko), astfel încît condițiile teoremei să fie îndeplinite și în consecință operația de alocare explicitată prin ( ) O reducere sensibilă a extensiei lui V se obține în cazul, în care în perechea controlabilă (J, KQ) matricea J este ciclică, ceea ce în fond este generic adevărat Atunci (demonstrația iese din cadrul prezentului curs), există adică există / cu • I > ’ nP — - ж f f ! ♦ jit '/ astfel încît ( , Ao) este controlabilă înmulțind ambii membrii cu l obținem * Г ( A ? Л \ i J -U ii*I ■ Vk = [ gl ț , q Particularizînd procedura din д і teorema pentru ( T), rezultă că extensiunea matricei V se face numai cu q coloane Să vedem acum în final care este interpretarea geometrică a teoremelor și de mai sus Presupunem că matricea V este monică și punem astfel incit CV П Yț Л Im B, este un generator al lui CV i= (ii) c matrice F e CpifV) astfel încît A + BF | CV este ciclică, atunci se poate construi un subspațiu cu dim Ю /-І •W * f C * f л і * X • л i •* j »’ • • ■ ■ » " ■ ■ ■ ■ * r • » Procedura de construcție este următoarea L Se determină prin procedura prezentată în subparagraful o acoperire de tipul de dimensiune minimă a lui Im E împreună cu matricele Presupunem că J eI vxv este ciclică Se determină cel mai mic întreg q pentru care ( ) este adevărată și se calculează matricele KQ, GQ, G^ Aceasta revine la a determina cel mai mic întreg q pentru care matricea Explicităm ( ) prin rezultînd ic; ' o o o Pentru = avem Ker[-FB]= { } Pentru = rezultă - ■ - - - a o o unde ( , Ko) este controlabilă Selectînd a doua coloană din și punînd - r rezultă perechea (/, /? ) controlabilă o - Г o - „ — ^ — "■'* Punem acum ) Г Si = , Rezultă din ( ) Г o o o o o prequn șl din J Г г f "V • Estimatorul / Tipul / (r л Di m en siu nea ■ ; L -II I n j de stare i Э # l ’ I * • ■ r *•/ • t n — p , • • /' Vo de funcțională — • * de F (uzual) • • * ^(^ — î) • ■ • de F (relativ minimizai) • * ( J > пг + • ( ■■ STABILIZAREA SISTEMELOR compensare dinamică) se formu- Problema stabilizării sistemelor (prin lează astfel : dîndu-se sistemul (neted sau discret) x' = Л x Ви, x e К", x( ) = x ”” Л'W Сат^Ло, adică sau Sd se construiască fie intern asirvi • ^mptohc stabil stabilizator ' * Саге asigură - G ГѵИЙ агеа matricei Л î==* 'Y (v- figura ), astfel încît sistemul țîn r • , /Т ^ndițm ( d ) Pnma ecuație "din^ л ^Se- b^ne i compens(aor ( ^ ftcînd « = F r ( - ^ de unde rezultă - + Să vedem {i îndeplinit™^ conditH sînt necesare iw P "«■ pentru aceasta s co-g № ea із (Э І ) se rescne în forQ^ det(AZ-^fl) £ yĂe ( ) obținem F * r • ' • > * t V * \ \ • t A ’ Jr Л = rang [XZ~A ~ sau VĂ e sau sau încă rang [AZ — Я II J/p‘ VAc de unde rezulta imediat că sau С\ЗД > c > A rang [AZ sau й аХГ“"âXelea™" ’f'rim “"*» ua, xa e evidență faptul că dacă alături adițional de sistem ul ( Л ) У a = ** rezultînd sistemul dinamic extins comandă după ieșirea sistemului extins ( ) atunci introducînd legea de ( ) se obține ( Л ) Egalitatea ( ) arată că problema stabilizării revine la introducerea unei reacții după ieșirea sistemului extins ( ) astfel îneît sistemul rezultant ( ) să fie intern asimptotic stabil Condiția „zonală" ( ) se poate înlocui cu o specificare precisă a locației polilor sistemului rezultant, rezultînd astfel Problema alocării prin compensare sau extensie dinamică, care se formulează astfel: i I Să se determine un întreg n (dacă este posibil cel mai mic) astfel incit pentru orice set simetric Лл de n + nc numere complexe să existe un compensator ( ) astfel îneît a( J яя Aa, ( ) II, Odată An specificat, să se determine efectiv sistemul (Лсл BCiFv,Gf-Lăsăm în sarcina cititorului sa demonstreze că condițiile necesare de soluționare a problemei alocării prin compensare sau extensie dinamică sînt: (j) (A, ) controlabilă: (jj) (C, A) observabilă (indicație: se utilizează aceeași tehnică dc demonstrare ca și aceea utilizată pentru obținerea condițiilor (i) și (ii) în cele ce urmează, ne concentram atenția asupra rezolvării acestei ultime probleme Pentru aceasta, avem nevoie de un rezultat pe care îl exprimăm fără demonstrație, aceasta necesitînd considerente matematice suplimentare • Teorema Fie (Л> JR, C) un sistem controlabil și observabil avînd indicii decontrol abilitate și observabilitate vc și respectiv v Atunci a pentru aproape toate perechile (J В, C) un sistem controlabil și observabil avînd indicii de controlabilitate și observ abilitate vc și respectiv v Atunci ♦ a bentru aproape toate perechile (K,g) (adică aleator alese), К gelHLw sistemul cu o comanda (Ло, b, C) ' unde ЛДЛ -f- BKC ( )- b^Bg' ( ) este controlabil și observabil, avînd indicele de observ abilitate (cel de cortro- labilitate fiind evident n); dualizînd avem și î pentru aproape toate perechile (K, h) (adică cu o ieșire (Ao, B, cT), unde Ло a fost introdus este controlabil și observabil, avînd indicele de control abilitate vc (cel vabilitate fiind evident n) Conținutul teoremei este sugerat grafic în figurile , a și b ~ în figura , a este ilustrat punctul a) al teoremei, în sensul ducîrd legea de comandă (după ieșire) că intro- Fig ■sistemei rezultant (/l,t, b, C) este controlabil Și observabil în figura , b cete ilustrat punctul b) al teoremei, care de fapt se poate obține prin aplicarea unei legi de comandă de tipul ( ) sistemului dual (Ar, C ', BT) Cu rezultatul de mai sus avem • Teorema Problema alocării polilor prin compensare (sau extensie •dinamică) are o soluție clacă și numai dacă sistemul (А, В, C) este controlabil și observabil Dacă aceste două condiții sînt îndeplinite, o dimensiune eficientă nc a compensatorului ( ) care asigură rezolvarea problemei în discuție este nc = min (ve, •>„) — ( ) Demonstrație Necesitatea rezultă din discuțiile făcute mai sus Suficiența condițiilor va fi demonstrată procedural Vom considera eta-pizat următoarele cazuri: • Cazul Sistemul ( ) are o singură comandă fiind deci de forma i ’ I - V * ( ) •Considerăm un set simetric arbitrar A de n numere complexe Există atunci conform teoremei o lege de comandă w = fTx J • • ’ * ‘ І Л £ / l ьV Să evaluăm spectrul matricei ÂR furnizate de ( ) Rezultă atunci det (XI — Ap^ = \I — A - bnTC —(Я + mn ') C -bkT XZ — J — mkT ( } Scăzînd în membrul drept al lui ( ) din a doua linie prima bloc linie înmulțită la stînga cu V obținem cu m= Vb din ( ) det (XZ — Лл) = XZ - A - bnTC -Ы • / V ■r c ♦atunci se aplică procedura descrisă în cazul Dacă • - c'; ' л; :■ ’Tdy compensatorul ( ) avînd dimensiunea ,v — Teorema este astfel complet ■demonstrată Rezultă în consecință imediat SI Algoritmul de sinteză al compensatorului Se specifică două seturi simetrice A și A„ continînd w numere complexe fiecare Se calculează indicii de controlabilitate și observabilitate vt și respectiv v Dacă v Oc, se face v = v — și indicatorul = și se trece la pasul în caz contrar, se face v — v„ — și I — și se trece la pasul Sc restrînge tricei В este > A„ la primele v numere Dacă numărul de coloane al m se alege m mod aleator o pereche (/(, g) și Sc pune Aa = A -(- BKC, l) as Bf caz contrar, Se calculează fT (vezi procedurii) astfel îneît Pentru sistemul (Aq, b, C) se construiește un estimator de tipul și alocat cu Ao de dimensiune v, (J, kT, H, M, n ') Se pune 'C = + ™kT, Bc= H + mnT, Fe = gkT, Gc = К + gnr Dacă = , STOP în caz contrar se face e STOP Se dualizează sistemul inițial punînd se revine la pasul Compensatorul, sintetizat realizează alocarea AR — A (j A„ In consecință, alegînd A(jA în sau С\( ), în funcție de tipul sistemului ( T, В, C), procedura furnizează un compensator stabilizator Exemplul Fie sistemul Г o o • * Matricele de controlabilitate și respectiv observabilitate sînt ' - i -Г - - • Q = - - - rezultînd că sistemul este controlabil și observabil în consecință, problema alocării prin compensare poate fi rezolvată , I Avînd în vedere că v-~v o — vom face următoarea specificare de poh Ал == A U Ao = {•— lt Deoarece m ~ fi , avem К =- , # = Determinăm f ' pentru care a(A J- bf r) Л Din Л rezultă )= + X | W + Rezolvăm ecuația A % • Г o П ' roi J , к = * w • • • -б - • Ъв «J estimatorul obținem , Dualizînd rezultatele A o ■ Я F ’ ■ ' ТП = ' țunde b »огегржпЦе sistemului inițial nednaiiznt) - Obținem anfcl fn final compensatorul Ec *= A * == [ D(~H + Г ~ Л L- O J n G ’' i ІІ ■ : , •• , • * / * /• c -*• • t - - de funcțională min (*n ’o) - • A J * •“ Як лгй z со înlocuind ( ) în ( ) obținem sistemul; ( ) = —Dx + Drxe care descrie complet dinamica sistemului inițial în contextul exogen precizat mai sus " • ' с - ■ Introducem acum o sistematizare de notație, în sensul că se indexează cu mărimile și parametrii propriu-ziși ai procesului, iar cu mărimile și parametrii sistemului exogen Astfel, vom pune x^x, A^A, A ^EEV, B^B, C^-C, D^-D ( ) De asemenea, punem (abuziv ca notație) * -• • * • î* ■- , = EÎ, z^zz ■ ( ) Cu ( ), ( ) și ( ), ecuațiile ( ) se scriu acum în forma definitivă (și de altfel consacrată) — Dixi D x ( ) CU > I s(i!)xSo = O, VxM O S ( OS) Precizăm că proprietățile interne de care am amintit nuanțează problematica reglării, așa cum se va vedea în cele ce urmează Sc impune de la început o observație remarcabilă înlocuind în ( ) și ( ) pe ( ) obținem precum și t ( ) ( ) ( ) Introducînd transformarea — 'ț ( ) care efectuează în fond o permutare, obținem ( ) ( ) ( ) ( ) ( ) Inspectînd ( ) se constată că /îs ч j • • \r \ \ I • * * Cititorul este invitat să probeze echivalența condițiilor de tipul (R) din ( ) sau ( ) cu ( ), * Formulăm acum ( ), cu măsurată y) a'ștfel incit dstcmulr^^ A- proprietatea dc reglare și o proprietate menționat) anume cea dc stabilitate internă Așadar ' se cere ca ( ) ( ) (R) vezi ( ) ц Cele două tipuri de probleme PR și PRIS reflectă cerințe naturale funcționale, care nu iau in considerare însă și perturbațiile parametrice Acest ultim caz va fi tratat separat în subparagraful SOLUȚIA PROBLEMEI REGLĂRII (PR) în vederea soluționării PR formulăm • Problema primară a reglării (PPR) Dîndu-se sistemul ( ) tind cea de a doua ecuație), adică (excep- z = Dx se cere să se determine o lege de comandă (după starea x e + avem A/C = AF((AF\ Im B> + Ф) = = Af (Af I Im В > + Af( ) c= + ^ + Im В = = ( AF Im В > + CV — ПѴ • * г • % Г \ * și ( ) este demonstrată Din ( ) rezultă că ImBcC ( } în consecință, față de o descompunere de forma vom avea conform cu ui pentru orice F e Din ( ) se vede imediat că pentru orice F' eRfflX/> avem zl-ll T "li ! Z l S unde F’ F^\, rezultînd că d defÂnda de F‘ și deci de (ii) Față de această ultimă constatare, referitoare la A , și ( ) corespunzătoare acestui F, atunci obligator in general* rezultă că dacă ( ) este îndeplinită pentru un anumit F și considerăm structurile ( ) sau ( ) pentru orice F și deci în mod automat ( ) este adevărată pentru orice F I erna este complet demonstrată Cu acest rezultat avem • Teorema PPR ars o soluție ( ) ă«ca și numai dacă Я +(Л) a + ( ) unde = sup {A: B; Ker Z>) * ’'хИл и» /? ( (irt> vțJrJHcft cot^ Т «*— пгт II Dar ( ) și deci Din ( ) și ( ' ) obținem O>+(A -H BF) с ф* Atunci cu atît mai mult, pentru acest F avem ( ) + W°™aPr™ SUcil( ? ») ÎS?‘adevărată reni®“ “"“r ?A + *• fece “( - ”* iîssVu s шto particul“ штм Нй în,m°d din nou și anume S ( - ) pe care le reluăm cu Im В = (V* n Im В ® Im В (X> = ф* ф j ( ) ( ) Im В c J ( ) precum și un Fo structura е * + =^* + = Іт[Г*і ?] ( } unde V* este o matrice bază pentru B) Atunci, cu referire la structura ( ), ( ) se detaliază în ( ) unde BS este matricea de controlabilitate a perechii controlabile (Л , Bsl), fiind deci epică, fapt care a permis scrierea ultimei egalități din ( ) Ultima explicitare în ( ) a matricei bază pentru ф* + | Im By arată, cu referire la structura ( ) că = (A+ BF )j(^* + Яяг Д' ( ) sau ( ) Rezultă atunci din ( ) ■^^ | Bf( ~ •ОЙ ФЯ^,(/) —> cu t —> oo în conformitate cu ( ) Din ( ) rezultă că F introdusă este o soluție a PPR, teorema fiind complet demonstrată ( ) prin ( ) Se poate trece acum la rezolvarea PR Vom distinge două cazma Cazul Perechea (С, A) este detectabilă Avem în acest sens • Teorema Presupunem că perechea (C, A) este detectabilă Atunci PR are o soluție dacă și numai dacă PPR are o soluție, adică dacă si numai dacă ( ) este adevărată Demonstrație Necesitatea Fie (nc, Ac, Bc, Fc, Gc) o soluție a PR Atunci obligatoriu Ж+(Л ) cz Ker ( ) (v ( )) Din ( ) obținem exprimarea ( ) unde și deci ( ) se poate scrie în forma cz ® ( - ) К De asemenea, dacii F„* este o baza a lui (Q* avem zi v* = v: î I В H ' A' t) ’ eJ I eJ c ( ) Punem ад* = o Cu ( ) și ( ), obținem din ( ) ( ) ( ) Punînd atunci СО Im V și COc A Im Vc obținem din ( ) și ( ) că = [СО + ] © ( ) cu СО + с Ж Dar Ж+(Л, ) = «+(Л) © Q C ( ) după cum rezultă din ( ), și în consecință ( ) se poate scrie ] © c ] © a>c w ( ) ceea ce implică în particular că incluziunea ( ), este adevărată, ultima incluziune în ( ) provenind din ( ) Trecem la demonstrarea suficienței condiției ( ) Aceasta fiind îndeplinită, există o soluție ( ') a PPR Întrucît (С, Л) este detectabilă, există L astfel încît | ( ) unde Ol Aq Ker CA' ( ) î^O ' este subspațiul neobservabil al perechii (C, A) Fie acum descompunerea si tem, restr,n (Л, n, ct Л), ch ГЛ -|- nG„C Ли ~| B/? л' - ТЫ Л/-’ /// Л, 'X ur de acum, evident, pwe- ' O ( ) Făcînd în sistemul ( ) inițializarea **■ rezultă cu ( ) că ( ) *( == РФ^( x ( ) Cum ( cu; t - oo, tfcQ g atunci datorită lui ( ), presupunînd că ф г( este nesingulară, ( , ) obținutenrin thscrpiil ' '•Z,f £,’^em(^or netede sau al sistemelor discrete ( ), că dxsciellzarea sistemelor netede, rezultă din ( ), cu ipoteza ( ) ceea ce revine la (vezi ( )) sau^chivalpitf ^eplinită, care în fond este generic adevărată), ( ), sistemul ’ tG îndePhmta' «sternul ( ) este reductibil la Л Г -|~ Jiu Cz c£ (X“ fite detaS? * Ча Ct'm ™ “““t* "» ®«i ₽««- S даrMul w ““e:i ъ ^ pentru sistemul roslrîns are o soluție strarpTSî!? *°Mia(se construiește pentru sistemul restrîns (v demonstram suficienței în teorema ) ' гіеі Ѵе ' тчІ /г' І;\іСГ\ Se obține J£er Q = KerC A Ker CA — Ker cs= I П o o o ~ Se calculează С+(Л) =Ker р+(Л) la forma diagonal canonică, obținem rl unde pi este polinomul minimal Aducînd XZ - чіе nude p (X) = X (X — ) și deci (față de *• -e W Se constată că perechea (Л , Bs) este ■[ : - este o reacție care stabilizează perechea (As, Bs) în sensul că Soluția PPR este atunci Rezolvăm acum a sistemului Alegînd PRfpriu construcția unui estimator stabil de tirul pentru starea > rezultă a(A LC) — { - , t, — , — } Obținem atunci în final I l > Л * ■ Л I Inspectînd compensatorul obținut, se constată că partea sa neobservabilă poate fi eli minată, rezultînd un compensator de ordin mai mic cu o unitate -Fc=[— - ], Gc= unde s-au folosit aceleași notații Funcția de transfer a compensatorului este = Fc (si - * Cititorul este invitat să verifice, aplicînd teorema valorii finale, ca z(t) cu t —> oo & de asemenea să constate că sistemul in circuit închis nu este intern stabil Cititorul poate constata din dezvoltările de mai sus că în rezolvarea efec tivă a PR, ipoteza de persistență ( ) nu a fost expres utilizată SOLUȚIA PROBLEMEI REGLĂRII INTERN STABILE (PRIS) •’ | ' ' *■ T * • • « “ * • — p А Х > unde ipoteza ( ) este adevărată, se cere sâ se determine legea de comandă ( '), care se sene detaliat S** v ’ a t j , •• • • ' ’ ' и = F^, -|- ( S ) astfel încît sistemul rezultant în circuit închis ( ), avînd în acest caz matricele ( ) dc forma particulară (compensatorul ( ) s-a redus la ( )) Лл = Л, + B^, Вя Л + BjF Dr= Dt ( ) (S) f * -e аеш„ se expHc, sau lim z(t) = , Ѵя e £ (obținută din particularizarea lui ( )) /ч Sue ,v SeS’e de Rularea PPR, care s-a făcut la nivelul sistemului (XI ) fără utilizarea expresă a ipotezei ( ), PPRIS utilizează în enun-țul acesteia forma detaliată ( ) necesara explicitărn condiției interne ( b Din ( ) se vede imediat că este necesar ca perechea (A,, BA să fie st abtlizabilă Pentru a simplifica expunerea, vom utiliza o ipoteză mai tare, care dealtfel este insă uzuală, și anume ^perechea (Лх, BT) este controlabilă" ( ) Inspectînd perechea (A, B) din ( , ) se constată imediat că ipotezele ( Л ) și ( ) se pot scrie compact în forma C\l/x(O) Punînd condițiile (S) și ( ?) din PPRIS se rescriu în forma (S) ЯР (А + BF) = (A | Im В > ( ) (Д) Ж+(Л + BF) g: Ker D ( ) Cititorul este invitat bă constate că cele două condiții (R) ale PPR și PPRIS, respectiv ( ') și ( ) sînt identice Așadar avem, într-o formă compacta ©PPRIS (reformulare) Cu ipoteza ( ) adevărata să se determine legea de comandă (după stare) ( ') astfel îneît condițiile (S) fi \R) explicitate prin ( ) și respectiv ( ) să fie adevărate în cele ce urmează ne vom referi numai la această ultimă formulare a PPRIS v v • Teorema PPRIS are o soluție ( ') dacă și numai daca exista ф e {A, B; Ker D) astfel incit r# == (p ф ( ) Demonstrație Condiția, este necesară Evident, avem CB =x (%*(A + BF) ф CB'fA - BF) ( ) \ = CBfA + BF) Ф •unde s-a utilizat ( ) De asemenea ( + BF) O *( + BF) c = F'| ( ) :SL ‘ * ■ > * iț *- - ' • yjy \ z •* , Fj (V = F" | O) ( ) Din ( ) și ( ) rezultă că ( - BF) | =( - BF') | -adică condiția ( ) este îndeplinită Datorită lui ( ), subspațiul — Im , = rB* - (A | Im B> ( ) Corespunzător bazei lui în care este reprezentată structura ( ) avem Im + nV — nZ unde nh = dina *, nz dim SL, nv unitate, precum și Вн АЯВН Bs dim — Im X X ЛВ, i ’i В и ( ) cu X ilgmîmlu se blocurile nesemnificative Alunei din ( ) multă că și ( ) ф* I- Im ( ) Conform ( ), matricea din membrul drept al lui ( ) fie epică și deci Tflllg [ * AgBg “ ^ trebuie să ( ) cu »=» dim , ' adică „perechea (Яя, b\v) este controlabilă" ( ) dacă și numai dacă ( ) este adevărată Conform cu ( ), există un Fs (aleator ales) astfel îneît făcînd modificarea As Bz) corespunde subspațiului f| (Indicație-, se scrie matricea de controlabilitate a perechii (Ap, BG) din ( )) Așadar avem acum descompunerea Ж = бе*ф^П ф ;ф ф T ( , ) X ==: ^ П © S ( ) * л f Л Ч •• * • r ф ( ) rezultă că in (subspațiul SL, & П este un subspațiu Az ilariant și în consecință are sens să ne întrebăm dacă £(\ > descompun: £ relativ la A , adică dacă subspațiul Ъ se poate modifica astfel incit rezultă că în (sub)spațiul SA П este un subspațiu ceea ce revine Az — la faptul că ecuația ( ) ( ) M aibă o soluție (testul Sylvcster) Dacă ( , ) este îndeplinită, atunci utilizînd o transformare de tipul o я o ' o o - ( ) structura ( ) capota forma (finală) Cu dezvoltările de mai sus și cu discuțiile privind subspațiul SA fi (A| Im B) care este Az invariant la nivelul lui în sensul că dacă este adevărată pentru unF (~A(CO*} aceasta continuă să fie adevărată pentru orice F e(J\CO*)> ceea ce rezultă imediat din ( ) Verificarea condițiilor (i) și (ii) este eminamente constructivă, ea revenind la construcția structurii ( ), urmată de verificarea valabilității condițiilor ( ) și ( ) O formulare esențialmente geometrică a condiției (ii) este: în ££/ *, ((£>* П ф СО ( ) Fie e ȘÎ F" е Г^( «F | ( з) F | СО = F* | subspațiul cu subspațiul СО, care conform cu ( ) este AF invariant, rezultă GL relativ la e Atunci se poate se poate alege aleator un ~ unem Atunci, prin procedura într-unul izo- d* u ( ) f cu ^ ( ) că r\ / л I ț tn sensul că dacă ( ’ Ju i ( ) (unde se reamintește că ipoteza ' ) continuă să fie adevărată) Aplicînd atunci teorema cu referire la ( ) și la ( ) rezultă că există (Vf e € (Af; Be; Ker Ds) astfel încît e ( ) Unnînd aceeași procedură ca în teorema , rezultă din ( ) că există СЪ ед{А, В, Ker D) astfel încît ( ) să fie adevărată Suficiența Procedăm similar ca in teorema , implementînd soluția и = Fx — F - F x a PPRIS cu ajutorul unui estimator, ales pentru ușurința expunerii, de stare și de tipul Atunci, demonstrarea condiției (R) se face identic ca în teorema Pentru îndeplinirea condiției (S), cit forul este invitat că probeze că matricea AR a sistemului în circuit închis poate fi adusă la o formă echivalentă de tipul ( ) unde, evident о{Аг + ВЛ) = , а(Л + С)сС sau U^O) Cazul » Perechea (C, A) nu este detectabilă în acest caz, așa cum s-a văzut cu ocazia rezolvării PR, incluziunea ( ) este o condiție necesară pentru îndeplinirea condiției (B) și în consecință continuă să fie o condiție necesară și pentru soluționarea PRIS în cazul PRIS valabilitatea condiției ( ) are unele implicații structurale mai profunde Intr-adevăr, dacă (нс, Л c, Bc, Fc, Go) este o soluție a PRIS, atunci îndeplinirea condiției (S) cere ca această soluție să fie și un compensator stabilizator pentru sistemul (Л,, B,, Cri în consecință, este obligatoriu (v § ) ca „perechea (Q, A j) să fie detectabilă" * ( ) Atunci, dacă stările ncdpcHWb i (von ( )) vom avea, contam си^М^б) ’ ' A) £°™еай « în consecință, se poate face descompunerea Ж — Жі Ф Я ф Ж, Ж = Ж ф « ( ) unde și Q'j smt evident A invariante Atunci într-o bază compatibilă cu ( ) și ținînd cont de condiția ( ) obținem sistemul echivalent ’ k ' ( ) Din ( ) se vede imediat că sistemul redus este ( ) unde perechea (С, Ă) este detectabilă și unde sistemul propriu-zis (Лъ Blt C DA este intact în consecință, ipoteza, privind controlabilitatea perechii (A , Bf) se conservă și la sistermd redus Cu această ultimă constatare, se vede că din punctul de vedere strict procedural, soluționarea efectivă a PRIS revine la teorema Exemplul Să se rezolve PRIS pentru sistemul (de rejecție asimptotică) din fig peicurbația x\ fiind de tip treaptă Ax, DsO Fig hi ! jvJ, ; MttiHceJd А, b, C, D sJnt (exercițiu pentru cititor) г : i - ’ i o F В *= ♦ * ♦ •** Lo oj J c = [ - o : в] D [ jmm ; ] Іідеоіѵ&т mal înfii PPHIS » Construim structura ( ) Să observăm că deoarece q~ , sigur * = Apllcfud procedura din paragraful , obține: Hezolvlnd ecuațiile obținem Cum punem LO de unde o o o- - o o o Alegem obținem un F astfel îneît stabilitatea structurală a proprietăților (S) și ( ?) create de o soluție a PRIS, atunci cînd parametrii sistemului (procesului) și ai compensatorului suferă abateri în jurul valorilor nominale Parametrii unei structuri de reglare sînt elementele matricelor A, A, Л Q, Ca, Dv Z> , precum și Ac, Bc, Fc, Gc ’ De la început vom considera fixe elementele matricelor A , Clt C , Dlt D , adică vom fixa clasa mărimilor exogene și vom admite că precizia’ traduc-toarelor și a elementelor de comparație este ideală în raport cu gradul de cunoaștere al modelului dinamic propriu-zis în consecință, vor fi luate în considerare numai variațiile parametrilor procesului propriu-zis și ale compensatorului, adică variațiile elementelor matricelor Л , Ac, Bc, Fc, Gc In acest context, vom spune că un punct '• * \ \ \ , « л unde V este soluția menționată mai sus Ținînd cont de a doua ecuație în ( ) obținem ț = ARț, Ц ) = î; ( ) z == Dp^ - {ERV - D )x de unde ^( = DR&AR(t) ~ (DRV * D ) Фл (^)^ о* ( ) Cum Флд(£) cu t—* со, conform proprietății (S), iar A este instabilă, (R) are loc dacă și numai dacă DRV - D = , adică (ER) sînt simultan satisfăcute de V, r •? :? , unde cp are drept coordonate o parte din parametri menționați, dacă ER continuă să aibă o soluție într-o vecinătate (parametrică) a lui rp Așadar, dacă ER (v ( )) are o soluție în o, în cazul neted , în cazul discret Cu ( ) ££ devin (Лв —- )Jnn)V = ~ В ( ) ( ) X unde, datorită condiției (S), sigur ,,A — \Inn este nesingulară" ( , ) Detaliat ( ) se scrie Intrucît A intervine liniar în BR, stabilitatea structurală a ER în Л + ДЛ cu ДА arbitrar, ecuația ( ) în creșteri trebuie să aibă q soluție ДУ oricare ar fi ДЛ Această ecuație rezultă imediat din ( ) scriind-o în variații și reamintind că toți parametrii, cu excepția lui Л sînt considerați ficși Obținem atunci ar fi ДА Aceasta se întâmplă + В fi fix ecuație care trebuie să aibă o soluție oricare numai dacă Datorită lui ( ), rangul este n!t = “ uc ținînd cont și obținem atunci din ( ) matricei din membrul sting de structura membrului drept lui lui ( ) ( ) + »«= w, rang - XU ' S •sau încă rang — n Din ( ) avem și ( , ACVlxb L^c^l — "C Comparînd ( ) cu ( ) rezultă în mod obligatoriu că există o matrice Tt convenabil aleasă, astfel îneît [£>i ] = ’c[KC ( - ) Din cea de a doua ecuație ( ) rezultă că lm c Im Dx și în cor' г-•cință ipoteza devine Ipoteza f Matricea Dr este epică Cu ipoteza ' rezultă din ( ) că „Tc este epică' ( ) Din ( ) obținem ( J ) Considerînd acum ultimele două ecuații din ( ) și scăzînd a doua ecu- -u j, în prealabil înmulțiră la stînga cu Tc, din a treia, obținem cu ( ) că ( i) Ținînd cont de ( ), obținem din ( ) și ( ) dependența remarcabilă [Dx D \ = S[C C ] «an sau încă, prin înmulțirea la dreapta ambilor membri ai lui ( ) cu v Vom spune că eroarea de calitate z este deductibilă din eroarea măsurată у dacă există o matrice S (epică), astfel îneît ( ) (sau echivalent( ) sau ( )) este adevărată »« » “„„pictare , a liniilor а/л> Altfel teia astfel încît matricea c , - ; este nesingulară Considerăm compensatorul echivalent prin transformarea ( ) * * o* - o * / C* și obținem cu ( ), ( ) și ( ) ( , ( ) ( ) ( ) Utilizînd aceleași notații rezultă că compensatorul arc structura a " * Ca Х/г cu starea Din ecuațiile generale ( ), obținem cu pensatoi'ului în forma detaliata ( ) și ( ) ecuațiile com- x R v ( ) unde în cea de a doua egalitate a lai ( ) s-a utilizat ( ) Inspcctînd ecuațiile ( ), ( ), se constată, că sistemul compensator încorporează wt subsistem dinamic descris de ecuația ( ), subsistem a cărui dinamică liberă (x«, ■ - ХлА) este impusă dc o copie, multâpldcată de q ori, a valorii proprii a sistemului exogen (Az jBe) (Cc, Ae) sînt stabilizabilă și respectiv detectabilă Demonstrație într-adevăr avem din ( ) rang [Л c — vZ = rang — vZ (X - v)Z + q> Vv e q ( - ) l * '■ \ • вЛ* * д- Л • adică obligator numărul de comenzi nu poate fi inferior numărului canalelor de eroare de calitate Datorită ipotezei de deductibilitate, nu este necesară pentru proiectare decît specificarea tripletului (Ль Blt Cf) Să analizăm cazul cînd ipoteza de lucru ( ) nu mai este adevărată Scriem ecuația modelului intern ( ) în forma Лс ~ x^ ^c CU unde în cazul în care ipoteza de lucru ( ) este adevărată avem ( ) ( ) M = X, N = Dm ( ) se vede că polinomul minimal al lui Л„ ( ) polinomului minimal al lui Аг este arbitrar în consecință; dacă în general, polinomul mi wCLlCcl ^(v) ="Yo + YiV -|- + Yr-iv'- + v ( ) perechea (M, N) din ( ) devine M = ( ) ——Yo—Y •••—Yr ! ] ] ±^;leJ - ^ ( ; ?- rămîl?md intacte (demonstrația acestei ultime rmațn este mai pretențioasa, și nu face obiectul cursului) devin^aemn Я comPensatorului Mi syAOms j wn'Oy( Nțw Y C Bucur, C Popooa Gh Simlnn rit , ™гк, U gică, București, ’ Пт°' ЕШга ^«^ctică șl Pcdag S Căli n, С Вe e a Sisteme automate ailaptivo fi S C ă i n, С В e e a Sisteme automate complexe Editura Tehnică, București S Călin, M Tertișco, I Dnmitrache, C Popcea, D Popc Cu ОЫіті-zări în automatizări industriale Editura Tehnică, București, Г Du mitr ac h e Tehnica reglării automate Editura Didactică și Pedagogică, Bucu-rești, S FI or e а, I Dumitrache, V G ă b u r i ci, FI Mu n t e a nu, S Du mit rin, I C a t a n ă Electronică industrială și automatizări Editura Didactică si Pedagogică București, S F o r e a Elemente de reglaj și automatizare Editura Didactică și Pedagogică, Bucu- rești, A Halanay Teoria calitativă a ecuațiilor diferențiale Editura Academici R S R, A Halanay Ecuații diferențiale Editura Didactică și Pedagogică, București, M Hăngănuț Automatica Editura Didactică și Pedagogică, București, V I o n e s c u, С P o p e e a, Proceduri dc sinteză a sistemelor automate Litografia In- stitutului Politehnic, București, V I o n e s c u Complemente de teoria sistemelor Litografia Institutului Politehnic, Bu- Ж V I o n e s c u Sinteza structurală a sistemelor liniare Editura Academiei R S R , Bu- curești, V I o n e s c u Sisteme liniare Editura Academici R S R , Bucuiești, л , V I o n e s c u Introducere în teoria structurală a sistemelor liniare Editura Academiei R S R , București, V Io ne seu C Pop oca Optimizarea sistemelor Editura Didactică și Pedagogică București, R A К al m an, P L Falb, M A Arbib Teoria sistemelor dinamice Editura lei B ЙЛЕЬ ~ V RSzvan SIMIW ЛЛ* • *** «“ Eu""" R S R , București, H H Rosenbrock State-Space and Multivariable Theory New \oik С I Sa va n t Calcului sistemelor automate Editura Tehnică, București, , D St a no mir O Stănășilă Metode matematice în teoria semnalelor Editura Teh- nică, București, O Stănășilă Analiza matematică Editura Didactică și Pedagogică, București, I G h Șabac Matematici speciale Editura Didactică și Pedagogică, București, J Truxal Automatic Feedback Control System Synthesis Mc Graw Hîll, New Ycrk» j V S Vlad i mir o v Ecuațiile fizicii matematice Editura Științifică și Enciclopedică, București, W M W o n h a n Linear Multivariable Control A Geometric A pprcach Springer-Verlag, W A Wol owich Linear Multivariable Systems Springer-Vcrlag, L A Z a d e h, E Polak Teoria Sistemelor Editura Tehnică, București * * ? ч ш • • V Zoubov- Theorie de la comande Editions MIR, Moscou, , > , - - £ / A BL L л' \ Г Plan editură nr, A Coli de tipar , Pun de tipar Tiparul executat sub comanda nr la întreprinderea Poligrafică „ Decembrie й, str Grlgore Alexancirescu nr - București, Republica Socialistă România - •^ГТх^ - ^L & Ж' * hn%- OTSSfcpiBX^ n'Z к лЯ ? https://neculaifantanaru com/principiile-conducerii html https://neculaifantanaru com/en/leadership-principles html